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ABSTRACT 

The  basic  Mitrovic's  method  for  analysis  and  design  of  linear  feed- 
back control  system  is  studied.   Some  further  properties  of  the  basic 
Mitrovic's  equations  and  curves  are.  investigated  and  applied  to  analysis 
and  design  of  linear  feedback  control  systems.   The  method  developed  in 
this  paper  gives  more  flexibility  to  the  basic  Mitrovic's  method. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance 
and  guidance  given  him  by  Dr.  G.  J.  Thaler  of  the  U.  S.  Naval  Post- 
graduate School  in  this  investigation. 
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1.    Introduction. 

The  Mitrovic's  method  is  an  analysis  and  design  technique  of  linear 
feedback  control  systems  using  algebraic  equations  and  graphical  methods „ 
The  method  is  based  on  conformal  mapping  of  the  S-pl&ne  into  the  real 
domain  through  the  characteristic  equation  of  a  feedback  control  system^ 
and  the  basic  principles  involved  in  the  Mitrovic's  method  are  derived 
by  applying  the  basic  theorem  by  Cauchy  in  the  mapping  procedure „ 

The  fundamental  theories  are  adequately  developed  and  applications 
of  the  method  are  demonstrated  in  references  1  and  2„  The  main  advantages 
of  the  Mitrovic's  method  seem  to  lie  in  the  fact  that  all  operations  are 
carried  out  in  the  real  domain  and  that  it  has  potentiality  of  controlling 
or  adjusting  multiple  variable  parameters . 

The  analysis  of  a  system  in  the  Mitrovicl's  method  is  done  by  analyz- 
ing the  relative  location  of  the  Mitrovic's  working  point  M  with  respect 
to  Mitrovic's  curves s  and  design  of  a  system  is  to  reshape  the  configura- 
tion of  Mitrovic's  curves  and  the  M  point  associated  with  the  system  of 
interest  so  as  to  meet  the  design  specifications .  The  analysis  and 
design  techniques  are  deduced  from  the  fundamental  properties  of  the 
Bo  vs  B][  curves  for  5  =  5l  denoted  by  l-^  with  the  M  point  as  listed 
below, 

1)  If  the  M  point  is  located  at  a^sco^  on  the  Bq  vs  Bi  curve 
for  ^s  Si  then  the  characteristic  equation  associated  with 
the  1 51  curve  and  M  point  has  a  pair  of  complex  conjugate 
roots  r  a  _  ^ co-m  ±Jh^  ox^ 

2)  If  tangent  lines  are  drawn  from  the  M  point  to  the  BQ  vs  Bj_  curve 

for  $=  loO  (the  !/.*  curve)  s  then  the  negatives  of  the  slopes  of 
the  lines  are  the  real  roots  of  the  characteristic  equation,, 
1 


3)  If  the  M  point  is  located  on  a  region  enclosed  by  a  \$±  curve 
encircling  in  the  clockwise  direction,  and  if  the  M  point  is 
in  the  first  quadrant  on  the  Bq  vs  Bj_  piane3  then  all  of  the 
roots  of  the  characteristic  equation  have  ^  greater  than  ~$±. 

In  this  paper  further  properties  of  Mitrovic's  equations  and  curves 
on  the  Bq  vs  B^  plane  are  investigated  to  extend  the  method  to  problems 
which  are  not  treated  in  the  literature s  and  to  give  the  method  a  better 
flexibility. 

In  Section  2  the  V^   curves  are  studied  to  provide  a  procedure  to 
evaluate  frequency  response  of  closed  and  open  loop  systems  on  the 
Bq  vs  B^  plane  so  that  one  can  obtain  information  for  frequency  response 
and  time  response  simultaneously.   Constant  band  width  curves  are  intro- 
duced  on  the  Bq  vs  Bj_  plane  to  design  a  system  to  meet  the  band  width 
specification  precisely  with  other  specifications  such  as  ~$  and  Kv 
specifications . 

In  Section  3  fourth  order  system  curves  are  studied  to  provide 
fourth  order  system  charts  to  fit  any  fourth  order  system,,   The  analysis 
and  design  procedures  developed  in  references  1  and  2  can  be  worked  out 
with  precons true ted  curves  on  the  charts  and  also  a  different  approach  of 
design  technique  is  tried  with  fourth  order  charts . 

In  Section  4  an  analytical  design  technique  is  derived  so  that  the 
labor  required  in  constructing  Mitrovic's  curves  can  be  avoided  or  mini- 
mized. 

Thus  each  section  can  be  regarded  as  independent  from  another; 
however s  the  principles  developed  in  Section  2  may  go  with  any  of  the 
sections  to  give  information  on  frequency  response „ 


2o    Frequency  responses  of  closed  and  open  loop  systems  on  tfee  Mitro- 
vic's Bq  vs  B,  plane. 
Mitrovic's  curves  for  ^  =  09  the  i^  curves,  are  obtained  by  mapping 

the  imaginary  axis  of  the  S-plane  in  to  the  Mitrovic's  plane,, 
For  a  characteristic  equation  of  a  system 

Rs)  =  S^+a^  S"~'+  ■  ■  ■  +  a,  s  +  a0=  o  <2-i) 

Mitrovic's  B  vs  B,  curve  for  0^^,4=1     is  defined  by 


(2-2) 


where  Cdn   and  t,   are  the  natural  frequency  and  damping  rati© 
respectively  of  a  root  of  the  closed  loop  system*  and 
<A(.'S,~)   is  a  function  of  ^  only.   (See  Appendix  I.) 
which  is  a  curve  on  the  Mitrovic's  Bq  vs  B^  plane  independent  of  the  last 
two  coefficients  &q   and  a^  of  the  characteristic  Equation  (2-1), 

Mitrovic's  working  point  M  on  the  Bq  vs  Bj_  plane  is  defined  by  the 
coefficients  a  and  a  such  that 

M  =  M(a1,a0) 

thus,  for  a  given  characteristic  equation s  a  set  of  Bq  vs  B-,  curves  and 
an  M  point  are  uniquely  determined  and  conversely. 

For  a  particular  value  of  damping  coefficient  ^  =  0  the  equation 
defining  the  \q   curve  is,  (using  proper  values  of  ^(0)  and  replacing 
0)7  by  0d   the  frequency  on  imaginary  axis  of  the  S-plane  in  Equation  (2-2) )8 

B.  =  a,co2-arW^-i-a76o6-  •  •  • 

<f2-31) 


Because  of  the  nature  of  the  \o    curve,  the  frequency  dependence  of 
Equation  (2-1)  can  be  represented  by  the  lengths  of  vectors  from  the  M 
point  to  the  I<?  curve.   This  indicates  that  frequency  response  of  a 
closed  loop  system  can  be  evaluated  on  the  Mitrovlc's  plane . 

In  this  section  relations  between  the  frequency  response  (magnitude 
and  phase)  of  a  closed  loop  system  and  the  To   curve  with  the  M  point  on 
the  Mitrovic's  Bq  vs  B|  plane  are  investigatedo   Band  width  ©f  a  closed 
loop  system  is  evaluated  and  the  locus  of  the  M  points  for  a  system  to 
meet  the  specified  band  width  is  cons true ted . 


2-lo  Interpretation  of  equations  (magnitude  of  frequency  response)  on  the 

Bq  vs  B-^  plane . 

Consider  a  closed  loop  transfer  function  of  a  feedback  control 
system 

where  N(S)  is  a  polynomial  of  S  of  order  less  than  n. 
Note  that  F(S)  z  0  is  the  characteristic  equation  from 
which  the  associated  [~^  curve  and  M  point  are  determined „ 

The  magnitude  of  the  frequency  response  of  the  system  represented  by 

Equation  (2-4)  is 


|N(jq))  | 


\kl)l 


Rei^j^' j  and   I/m\Hj  D!]   can  be  represented  in  terms  of  aQS,  BQ  and 

a1 ,  B,  respectively  where  B,  and  Bq  are  the  Mitrovic's  variables  defining 

the  ("^  curve  which  are  functions  of  cO  defined  by  Equation  (2-3). 

Rewriting  F(S)  in  Equation  (2-4)  as  the  sum  of  even  F(S)  and  odd 
F(S), 

F(s>  =  (a&i-  aj  +  a+&*+  •  ■  ')+$(&!  +a3s2+a.5.s^  • .  .) 

then  P(j  CO)   r  j?e{  F  (j  CO)}  4-j  Iw{  F  (j  M)\ 

- 1 1    -  0%  -  •  *  )j  4-j  o)[a ,  -( a3  wa-  b.so)  +  -  -  -)] 

therefore  by  using  Equation  (2-3) 

Re{F(jco)}  =  a0-(<3a6o2-a<to)"+  •  -■•)  =a0-Bd 

l-m{ F  (/«»}  -  w  [a4  - ( a^of-  aj«% -  •  ^  =  a)  (a ,  -  Bi ) 

thus,  Equation  (2-5)  can  be  written  as 

\n(jCo)\ 


Kjco^ 


5 


(2-6) 


"1 

lows.   The  magnitude  A  at  frequency  Od    is  the  ratio  of  the  length  of  two 
vectors  N(jC0  )  and  (aQ  -  BQ)  +>CO(a1  -  Bj). 

Since  a-i  and  a*   are  the  coordinates  of  the  point  M(ajSa0)  on  the 
Bq  vs  B^  plane,  and  B\   and  Bq  at  any  frequency  are  determined  by  the  co- 
ordinate of  a  point  at  CO  on  f7  ,  the  numerical  value  of  the  denominator  of 
Equation  (2-6)  can  be  determined  easily. 

Suppose  the  To    curve  and  M  point  configuration  is  as  shown  in  Fig. 
2-1  and  choose  a  point  P  on  the  I^  curve  at  which  the  frequency  is  CO  8 
then  the  coordinates  of  P  are  (B^co),  Bq(o))).  Locate  a  point  Q  at 
(a,,  Bq(co))  so  that  Q  is  at  the  same  ordinate  as  P  and  is  on  the  line 
**1  ~   al  '    t*ien  t'ie  distances  from  M  to  Q  and  Q  to  P  are 

FlQ«|ao-B0(co)|  >   OP  =  |a  ,-6,(0)')  |     respectively. 
Locate  a  point  R  on  the  line  segment  QP  such  that  Q.R~Ct)'Q.p      and  asso- 
ciate a  frequency  oo  to  the  point  Rs  then  Q.R  =  CO-OP  =  (jO  |a,-8,(C0)| 
and  MR  =  (h&   +  G&   )  —  \(  a o  -  §> ( vX)  +  J  GO  ( a , -  B , (a) Y)  |    which  is  the 
value  of  the  denominator  of  Equation  (2-6)  at  frequency  a)  . 

As  CO  varies,  the  point  P  moves  along  the  17  curve  and  the  point  Q 
moves  along  the  B,  s  a,  line  maintaining  the  same  ordinate  as  P.   The 
point  R  will  describe  a  curve  which  is  represented  by  the  dotted  curve 
and  denoted  by  the  60  (aj_  -  Bi )  plot.   The  &)(a,  -  Bj)  plot  is  a  well  defined 
curve  if  the  Iq  curve  is  well  defined  and  has  the  following  common  proper- 
ties for  all  systems? 

1)  It  has  frequency  on  the  curve  as  a  parameter  such  that  two 

points  constituting  a  pair  of  points  one  on  the  |q  curve  and  the 
other  on  the  00(a|  -Bi)  plot  have  the  same  frequencies  with  the 
same  ordinates. 


CO(a  -B  )  plot 


'  \^  ''Cat  (Magnitude  circle) 


Fig.  2-1  Illustration  of  obtaining  Frequency 
response  for  a-system  with  no- zeros. 


2)  It  starts  (or  the  originating  point  is)  at  (a, ,0)  and  passes 
through  the  intersection  point  of  the  [^  curve  and  the 

B^  =  a-^  line. 

3)  The  points  for  CO  =   1  on  the  VQ   curve  and  0)(a,  -  B|)  plot  are 
identical. 

4)  It  lies  between  the  Q   curve  and  the  B-^  -   a,  line  for  C0<1 
and  dually. 

If  any  circle  centered  at  the  point  M  passes  through  a  point  at 
60  r  60i  on  the  &0(a^  -  B, )  plot  with  radius  of  r,  then  r.  represent  the 
value  of  the  denominator  in  Equation  (2-6)  at  60-  cOp  The  circle  cen- 
tered at  the  M  point  is  called  the  magnitude  circle.  The  numerator  of 
Equation  (2-6)  can  not  be  evaluated  from  the  1$  curve  as  the  denominator. 
It  is  possible  to  draw  separate  Mitrovic6s  curve  for  N(S)  and  obtain 
N(j"o))  but  this  is  labirious  if  N(S)  is  of  high  order. 

The  simplest  case  is  for  N(  jco  )  independent  of  frequency.  This 
is  the  case  for  a  system  that  has  no  zeros  in  its  closed  loop  transfer 
function. 

Suppose  a  closed  loop  system  has  no  zer©8  then 

a        <3q 

M  —   r-r =: ; ; . —  ©r 

\{&0S0)  H- j  CO  (a,-^  | 

The  magnitude  Ai  at  60=  00,  is  a^/r.  and  conversely  if  a  magnitude 
A  =  A,  is  given  then  the  frequency  to  give  A  s  A,  is  found  by  drawing  a 
magnitude  circle  with  radius  equal  to  aQ/A-|,  and  locating  the  point  or 
points  where  the  circle  intersects  with  &>(ai  -  B« )  plot. 

Some  particular  points  in  the  frequency  response  of  a  closed  loop 


system  such  as  the  resonant  pointy  zero  db  point  and  half  power  point 
(band  width)  can  be  found  very  simply. 

Drawing  a  magnitude  circle  tangent  t©  the  0)(a,  -  Bi)  plot  gives 
the  resonance  point.  Magnitude  circles  with  radii  of  a^   and  fl  &q   give 
the  0  db  point  and  band  width  respectively. 

For  a  system  with  N(S)  in  Equation  (2-4)  of  order  other  than  zero, 
the  procedure  is  slightly  modified  but  still  the  basic  idea  is  the  same. 

This  is  illustrated  with  Fig.  2-2.  Modification  must  be  made  to 
take  care  of  the  frequency  dependence  of  the  numerator  in  Equation  (2-6). 
Since  the  denominator  in  the  equation  defining  the  magnitude  of  fre- 
quency response  bears  the  same  form  regardless  of  the  order  of  N(S)S  the 
construction  of  the  (x)(a     -  B. )  plot  is  the  same  for  all  systems. 

Rewrite  Equation  (2-6)  as 

|(Vty+j^(-eo  I-tH^j^I  (2-7) 

and  interpret  as  follows : 

The  quantity  at  any  60  on  the  left  side  of  Equation1  (2-7)  is  deter- 
mined by  drawing  a  magnitude  circle,  and  it  is  represented  by  the  length 
of  radius  of  the  magnitude  circle.   On  the  right  side  of  Equation  (2-7), 
|  NfjcoOj   has  some  finite  value  for  every  a)   .  Then  a  number  which 
multiplies  to  ||\)(j'o))|   to  give  the  equality  in  Equation  (2-7)  is  the 
reciprocal  of  the  magnitude. 

In  Fig.  2-2,  the  value  of  |fS/(jCo)j  is  represented  on  a  line  drawn 
from  M  point  in  an  arbitrary  direction.  This  line  is  denoted  by  |n(j'0o)| 
and  at  each  point  on  the  line,  there  associate  a)  such  that  if  the  point 
N  is  at  a)   on  the  line  then 

MA/=|wCjW)|  and       A=MTJ/r 
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)  plot. 


Fig.  2-2  Illustration  for  obtaining  frequency  response 
for  a  system  with  zeros. 


Conversely,   0)  =  6),  for  given  magnitude  A  -   A.    is  determined  by  search- 
ing such  0)    that  for  two  magnitude  circles  with  radii  of  r  and  Ar  cross 
the  0)(a  -  B1)  plot  and  the  |f\Kja>)|  line  respectively  at  the  same  fre- 
quency.  This  can  be  done  by  a  few  tries  and  errors, 

The  arguments  developed  above  are  not  restricted  by  the  order  of 
system  equation.  The  only  restriction  is  to  have  V0    curves  plotted  on 
the  Bq  vs  B^  plane  with  the  same  scales  of  B^  and  Bq  axes  to  provide 
simplest  graphical  solutions „   Bq  vs  B  curves  can  always  be  plotted  on 
such  axes  scales  by  using  suitable  frequency  scaling  technique „ 

As  mentioned  previously,  the  higher  the  order  of  the  numerator 
polynomial  of  a  system,  the  more  complicated  the  method  is.   The  problem 
is  to  provide  a  convenient  method  to  represent  the  frequency  dependence 
of  the  numerator  of  a  closed  loop  system  equation. 

In  the  following  paragraphs,  such  methods  are  presented  for  the 
cases  of  first,  second  and  third  order  numerators. 


2-2.  Representation  of  frequency  dependence  of  N(S)  on  the  Bq  vs  B][  plane. 

In  the  previous  section,  it  is  seen  that  the  frequency  dependence 
of  the  numerator  polynomial  N(S)  in  a  closed  loop  system  equation  must 
be  represented  as  a  distance  from  the  M  point. 

In  Fig.  2-3  it  is  illustrated  for  N(S)  of  first  order.   When  a 
closed  loop  system  has  one  zero  then  N(S)  is  of  the  first  order  as 

M/r\   I  III  I 

The  real  term  of  N(j(W  )  is  independent  of  frequency.   By  construct- 
ing a  rectangle  with  a  fixed  side  of  length  equal  to  a~  and  another  side 
which  is  in  right  angle  with  the  fixed  side  varying  linearly  with  (x)   , 
N(joJ)    can  be  represented  by  the  length  of  the  diagonal. 

In  Fig.  2-3,  a  point  T2  is  located  on  the  Bq  s  a~  line  such  that 
MT£  =  b, .   If  the  line  segment  MT2  is  calibrated  linearly  with  Cd   such 
that  63  z   0  at  M  and  Cd  z   1.0  at  T2s  and  if  the  same  calibration  is  ex- 
tended beyond  the  point  T^s  then  at  eaeh  point  on  the  line  BQ  =  aQ  to 
the  right  of  M,  there  associate  CO  linearly  with  the  distance  from  M  to 
each  point.   For  any  point  T^  at   0)=  <^,  MT,  =  b^  6)^  .   Noting  that 
TT^  =  aQ,  the  length  of  diagonal  of  the  rectangle  MT,Ta,  is  equal  to 

MT=|a0+jbia>i|-  |n(jo>.)'| 

For  a  system  with  two  zeros  in  its  closed  loop  system  equation^  N(S)  is 
of  second  order  polynomial  as 

N(S)=kSVfc^+a0   and   |N(/'w)|=|(a0-l5Wl)-+-jbia)  | 
The  imaginary  term  can  be  represented  in  exactly  the  same  way  as  in  case 
of  first  order  N(S). 

The  real  term  is  now  also  frequency  dependent ,  and  can  be  repre- 
sented on  the  B,  =  a,  line  with  quadratic  calibration. 
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fcft— j 

4- 


|N(^)|=|a0+^lffi]J 


Fig.   2-3     Representation  cf   I  N(j €■*)-« |a0+;jbCo| 


*^- 


b  <i> 


Fig.    2-4     Representation  of    JN ( jco)  |  =|( aQ-b  CO  )+jb  coj 
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In  Fig.  2-4,  b]4)  is  represented  as  in  Fig.  2-3 „   To  represent  the 

2 
term  (a  -  b  cO  )  the  line  b^  =  a^^  is  calibrated  as  follows.   The  point 

PQ(als0)  is  labeled  0)z   0.   From  the  point  PQ  up  each  point  on  the 

B^  =  a^  line  is  calibrated  with  a)  such  that  if  a  point  P,  on  the  B,  s  a. 

2 

line  reads   a)  =  ti>±      then  the  distance  P^  r  bjw^  .   Noting  that 

SP"0  =  aQ  , 

By  forming  a  rectangle  MT,TP-,  ,  the  length  of  the  diagonal  MT  is 

For  a  system  with  N(S)  of  third  order  or  higher,,    )\J(i  co)   can  be  repre- 
sented in  a  similar  way  as  in  case  of  second  order  N(S)  with  some  modifi- 
cations as  the  case  of  first  order  N(S)  modified  to  account  the  second 
order  N(S);  however,  it  requires  laborious  calculations.   For  example, 
for  third  order  N(S) 

N(S)   =bJ53-f  b2SZ  +  b,S  4  a0  (2-8) 

N(  VcO  )   =  (a0-  fc>2 cxT  )  +  J(tj  GO-  b3  a/) 

which  indicates  that  the  Bq  =  aQ  line  must  be  calibrated  to  account 

3 

(b,G0  -  b3 GO  )  which  requires  cubic  calculations ,  and  quadratic  calcula- 
tions for  fourth  order  N(S)  and  so  on. 

Those  calculations  can  be  avoided  by  applying  Mitrovic's  basic 
equations  for  the  Yq    curve  to  Equation  (2-8)  and  by  proceeding  the  same 
way  as  was  done  to  represent  F(  jco  )  by  means  of  0)  (a^  -  B,  )  plot. 
This  is  illustrated  with  Fig.  2-5. 

Since  N(S)  is  of  third  order,  the  V0    curve  obtained  from  Equation 
(2-8)  is  a  straight  line  and  this  is  denoted  by  T^H    .   Let  Mitrovic's 
variables  associated  with  Equation  (2-8)  be  B^N  and  BqN.   Then 
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CO  (aj-I^L-T 


T2r- 


-0>(b1-Blv)    plot. 


/ 


/ 


/> 


x^COCb^Bj   plot  translated. 


// 


// 


^(j^) 


Fig.  2-5  Representation  of  |n( ja>)  =  K^-B^  )+ja>(b,  -BlH  ) 
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1^  is  defined  by 

z  f 2-9) 

The  point  M^b^  ,  slq)  which  is  analogous  to  the  point  M(a,  s  aQ)  for 
the  characteristic  equation  F(S)  r  0,  is  determined  by  the  coefficients 
b^  and  aQ  of  Equation  (2-8) o 

If  Bjh  and  BQN  axes  are  chosen  identically  with  B*  and  Bq  axe.ss  then 
has  a  slope  of  bj/b^  and  It-  is  located  on  a  horizontal  line  passing 
through  M  on  the  Bq  vs  Bj^  plane . 

If  0d(bi   -  Bjjg)  plot  is  constructed  from  £N  and  the  line  B^  =  F,N  = 
b^  ,  then  by  the  same  reasoning  as  for  G3(a^  -  Bj_)  plot,  for  any  fre- 
quency C0X  the  quantity   |  N  (j  a\)  |  =  \(&0-  E^)+-  j  fi>j  (  bj.  -  BJN)  | 
where  E^N  and  BQN  are  coordinates  of  a  point  at   6)=  CO,   on  [^  ,  is 
equal  to  the  distance  from  MN  to  T2  where  the  point  T~  is  at   6}=  GO.  on 
the  0)(b1  -  B1N)  plot. 

Since  N(j<i)  )   and   (a«  -  Bq)  +  j  CO  (a^  -  Bj_)   must  be  compared  to 
get  A  the  magnitude,  and  since  N(jco  )   is  measured  from  Mjj  whereas  the 
latter  is  measured  from  M,  it  is  necessary  either  the  C0(b|  -  B,J  plot 
is  translated  horizontally  to  the  right  or  the  C0(a|  -  B^)  plot  to  the 
left  so  that  MN  and  M  coincide .   In  Fig.  2-5 8  the  6d(b^   -  BlN)  plot  is 
translated  to  the  right  by  amount  of  a,  -  b,  ,   then 

|N(jco^|=  MJi  =TlT3 

The  method  representing  |n(.'C0  )|  for  N(S)  of  up  to  third  order  is 
presented.   For  high  order  N(S),  |N(j0O  )|  can  be  represented  by  an  analo- 
gous method  as  in  case  of  third  order  N(S);  however,  it  is  inevitable  to 
be  laborious  to  represent  N(jC0  )  as  the  order  of  N(S)  increases. 
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On  the  Bq  vs  B,  plane  for  feedback  control  systems  with  adequate 
relative  stability s   the  M  point  is  usually  located  in  the  region  near  to 
the  origin  as  compared  with  whole  region  covered  by  the  lo  curve.  This 
is  due  to  the  fact  that  I  -5  curves  for  higher  X,    are  biased  towards  the 
origin.   Therefore  in  obtaining  frequency  response  usually  a  first  part 
of  the  (^  curve  is  required. 

It  must  be  noted  that,  since  the  basic  Mitrovic's  equations  are  not 
restricted  to  the  characteristic  equation  of  a  feedback  control  system 
but  apply  to  any  polynomial  with  positive  coefficients s  the  principles 
developed  above  apply  to  any  transfer  function  which  is  expressed  in  the 
form  of  a  quotient  of  two  polynomials  provided  the  Bq  vs  Bj_  curve  for 
~$  =  Q  and  the  M  point  associated  with  the  transfer  function  are  defined,, 
Therefore  it  is  obvious  that  the  magnitude  in  the  frequency  response  of 
an  open  loop  system  can  be  evaluated  in  a  similar  way  on  the  same 
Bq  vs  B;l  plane  where  the  magnitude  of  the  closed  loop  system  is  evalu- 
ated. 

In  general  the  open  loop  frequency  response  can  be  obtained  by  con- 
structing the  lo  curve  and  by  locating  the  M  point  associated  with  the 
open  loop  transfer  function  and  then  by  proceeding  exactly  the  same  way 
as  for  the  closed  loop  system.   In  Section  2-7 9   some  examples  for  evalu- 
ating the  phase  of  the  frequency  response  for  both  open  and  closed  loop 
systems  are  shown,  which  will  well  illustrate  the  above  arguments. 


2-3,  Determination  of  band  width  of  a  closed  loop  system  on  the 

Bq  vs  B,  plane. 

For  any  given  system,  band  width  defined  at  magnitude  of  l//z     or 
half  power  point  can  be  determined  precisely  by  the  method  described  in 
Section  2-2. 

For  a  system  with  no  zeros,  band  width  is  found  at  the  intersection 
point  of  CO (a,  -  B,)  plot  and  a  magnitude  circle  with  radius  equal  to 
/2aQ. 

For  a  system  with  zeros,  band  width  can  be  found  by  trying  a  set  of 
2  concentric  magnitude  circles  with  radii  of  r  and  fl  r      such  that  the 
two  circles  cross  the  6)(ai  -  Bj_)  plot  and   JN{jcO  )|   respectively  at  a 
same  value  of  frequency. 

Frequently  in  the  study  of  feedback  control  system  the  point  of  band 
width  alone  of  the  frequency  response  is  of  interest  and  also  sometimes 
upper  bound  and/or  lower  bound  of  band  width  is  of  interest. 

By  the  property  of  the  equation  defining  band  width  and  associated 
geometry  on  the  Bq  vs  B^  plane  it  is  possible  to  set  up  upper  and/or 
lower  bound  of  band  width  for  a  given  system,,   Among  the  various  points 
in  a  frequency  response  curve  of  a  given  system,,  if  band  width  only  or 
approximate  value  of  band  width  is  required 8  then  by  setting  up  the  lower 
and/or  upper  bound  of  band  width,  the  effort  in  plotting  £>(ai  ~  Bx)  can 
be  economized  or  even  without  plotting  any  portion  of  60(ai  "  Bi)  trae 
approximate  value  of  band  width  can  be  determined. 

In  this  section  the  three  cases  for  a  system  with  (1)  no  zeros, 
(2)  with  one  zero,  (3)  with  two  and  three  zeros  are  considered. 

For  a  system  represented  by  Equation  (2-4) 8  and  using  Equation  (2-6), 


the  band  width  defined  at  half  power  point  is  the  frequency  to  satisfy 
the  equation 

|crf.Nl      1n(jq>)1        j_ 
|(a0-B^+jw(a,-B,)|  =/a  |n(jco)|  C2-io) 

All  that  has  to  be  done  is  to  interpret  Equation  (2-10)  on  the  B«  vs  B-. 
plane.  The  left  side  of  Equation  (2-10)  can  be  interpreted  by  the 
60  (a,  -  B,)  plot  as  shown  in  Section  2-1 ,  and  for  the  right  side  is 
needed  a  modification  to  the  method  described  in  Section  2-2  to  account 
for  the  factor  /T. 

Consider  a  closed  loop  transfer  function  with  no  zero. 

Let  0^  be  the  band  width s  then  from  Equation  (2=10) 

|^-^ra)^Vj^b(a.-Bj(^)|=/2^  (2-11) 

must  hold9  where  B^(  60^)  and  BQ(  6)^)  are  the  coordinates  of  a  point  ©n 
the  \0   curve  at  C0z    0)b  . 

Equation  (2-11)  implies  that  a  magnitude  circle  with  radius  of 
f~l   aQ  intersect  with  the  C0(aj_  -  Bj_)  plot  at  Cd^     which  is  described  in 
Section  2-2.   Since  the  \0    curve  tends  to  infinity  in  some  direction  as 
Cd   tends  to  infinity s  G)(a]_  -  B]_)  plot  also  tends  to  infinity.  Noting  the 
finite  radius  of  a  magnitude  circle  and  that  the  U)(ai  -  Bi)  plot  starts 
from  a  point  inside  the  magnitude  circle  with  radius  fl   a^,  it  is  seen 
that  there  exists  at  least  one  intersection  point  between  the  circle  and 
the  6)  (aj_  -  B^)  plot.  Depending  upon  the  relative  location  of  the  M 
point  with  respect  to  the  V6   curve  and  also  depending  upon  the  geometric 
shape  of  the  Vq   curve  there  may  be  more  than  one  such  intersection  point 
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implying  more  than  one  band  width  to  exist. 

On  the  BQ  vs  B,  plane  with  the  \0   curve ,  by  drawing  a  magnitude 
circle  of  radius  fl   aQ  and  by  inspecting  the  values  of  6t)  at  all  inter- 
section points  of  the  circle  with  the  \$   curve  if  there  are  any9  and  by 
estimating  the  shape  of  the  CO  {a.-,    -  B-.)  plot.,  information  about  the  band 
width  can  be  obtained.   Also  the  point  at  CO  s   I  on    \q   and  the  points 
where  the  1^  curve  intersects  with  Bi  =  a, 9  Bq  s   a^  and  Bq  i  (1  +  /~2)a0 
lines  provide  additional  information. 

Suppose  a  system  is  represented  by  the  l^  curve  and  M  point  as  in 
Fig.  2-6,  A  magnitude  circle  C.b  with  radius  of  (l  ag  is  drawn.  The 
points  P^j,  ?2i  **3»  ?4»  ant*  ?5  represent  the  points  where  the  1^  curve 
intersects  with  C-  ,  BQ  s  aQ  ,  BQ  =  (1  -r  /2)a  or  B^  -.  a^  . 

Suppose  at  each  point  of  those  intersections,,  frequencies  are  read 
to  be  60l  ,  CO  ,  Q l2  60  and  &)_  as  shown. 

The  latter  part  of  the  \$   curve  is  neglected  since  it  is  assumed 
that  for  that  part  of  \Z   9    the  CO  (a-,  -  B*)  plot  stays  outside  of  C.. 
In  an  actual  situation  this  can  be  judged  by  looking  at  the  quantity 
C06(a,  -  BjCfiOg))  and  GOyCa^  -  B  (&)y))  where  <5)&  and  CO-,   are  the  fre- 
quencies at  the  points  where  the  latter  part  of  1^  intersects  with 
Bq  =  (  1  +  /~2)aQ  and  BQ  =  a^  respectively. 

The  C0(a,  -  B,)  plot  resulting  from  the  portion  of  the  1^  curve 
above  the  line  Bq  =  (  l+/2)aQ  can  never  cross  CAb  ,  therefore,  then 
there  are  four  possibilities  for  band  width  as 

(1)  6^60^^  (3)   0)^(301,^6)2 

(2)  C0^  60b^  C05  (4)        0db^0dt 
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2.414a, 


Q4                                                              B0=(l+/2)a 

^ — 
p2>4        ■    _* 

0 
3            ^x     Mag.   circle 

VAb 
M                            1                          Bn=an 

/ 

r=/za0      J 

Fig.  2-6  Determining  approximate  "band  width 
for  a  'system  with  no  zeros. 
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Each  of  the  four  possibilities  above  can  be  evaluated  by  looking  at  the 
point  at  £0  =  1  on  the  1^  curve. 

Suppose  cDg^"!  then  the  fi)(a^  -  Bi)  plot  for  0)~  GV  stays  outside 
of  the  Ifl  curve  in  turn  outside  of  Cajj  *       indicating  63b  <  6)3  . 

Suppose  6)3 <  1  then  Q3P3  >Wj(a1  -  B  (A>3))  and  the  fc>  (a,  -  B,) 
plot  must  cross  C  ,  somewhere  between  63  ~  and  Cl),    therefore 
C^<G)b<ca«<l      Or     6)3  <6>b<i<a>* 

Suppose   6)=  1  is  located  between  ^  an^  ^3  or  tdy<  \  <  C0*>    3 
then,  since  the  i~^  curves  for  0)2  <  03  <  Co  3  are  inside  of  CAb  ,   the 
CO  (a,  -  B,)  plot  must  cross  C.^  between  4)  r  1  and  6)^  ,       therefore 
1<  6),<  0)3°   In  this  case  there  is  still  a  possibility  that  u3b  <  0)3 
but  this  can  be  checked  by  looking  at  the  quantity  03{a,  -  B«)  for  some 
(£   less  than  &>«  „ 

Suppose  OX  <  1  <  <D2  ,   then  the  a)  (a^  -  B^)  plot  stays  outside  of 
CAb  for  1^C0<C02  and   0)^0)3  •   The  portion  of  the  43(aj_  -  B^)  plot  for 
0d2<  0^<  Od-x   maY  be  inside  of  CA2  .   If  this  is  the  case  then  6J>2<  ^b*  ^3 
again  the  possibility  that  GJb<  C02  must  be  checked.   If  0)(a1  -  Bj_)  for 
(x>2<  (J^<  Oh   nas  no  P°i-nt  inside  C.,   then  0)(a^  -  Bj_)  for  0)£  1  is  out- 
side CAb  and   ^^1  .   If  the  points  P,  and  P2  are  very  close,  then 

Suppose  03  >  1 ,  then  there  is  no  possibility  that  CO,    is  in  the  in- 
tervals (1,  C02)   and  0)>  0)3  .    Since  the  OX®!  -  B1)  plot  at  03  =  0  is 
at  the  point  (a^t   0)  and  the  63(a,  -  B,)  plot  at  63  ^  1  is  outside  of 
CAb,  band  width  is  definitely   Q<  63  <  |  <  00}   and  possibly 
0O2  <  OJ^  Cc>3- 

If  either  one  of  the  frequencies  CO     and  633  is  equal  to  ls  then 
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00b  =  1  and  other  possibilities  can  be  investigated  as  above . 

Thus  the  four  possibilities  listed  above  are  analyzed „   The  geometry 
involved  in  Fig.  2-6  is  so  simple  that  the  relative  locations  of 
6i(a1  -  B^)  plot  for  C0*0)4,   6)34oKg)4,  (O^COSO^,  CO  <  6)  $  CO 
and   0*  C0S&)   with  respect  to  C^b  are  determined  almost  by  inspection 
using  the  point  C0=  1  as  a  guide,  resulting  in  determination  of  upper 
and /or  lower  bound  of  band  width . 

It  must  be  noted  that,  depending  upon  the  system  of  interest^  there 
will  result  various  configurations  of  the  ["J  curve  and  the  working  point 
M  (a, ,3q);  however,  the  basic  principle  of  reasoning  is  the  same  as 
above . 

For  a  system  with  one  zero,  the  band  width  defining  Equation  (2-10) 
is  to  have  frequency  dependent  terms  in  both  sides  of  the  equation  and 
requires  some  modification  in  the  method  applied  for  a  system  with  no 
zeros.   If  a  system  has  one  zero  in  its  closed  loop  transfer  function,, 
then  the  equation  defining  band  width  is 

|(a0-BD(a)^Vj£«>b^a.-B,(a)^|=/2  |aD+jb,a)J      (2-12) 

where  b^  is  the  coefficient  of  the  S  term  in  the  numerator  of 
the  given  closed  loop  transfer  function   BqCoO^)  an^  biC^k) 
are  the  same  as  in  Equation  (2-11). 
Suppose  the  [J  curve  and  M  point  configuration  is  as  shown  in  Fig» 
2-7.  To  account  for  the  frequency  dependent  term  fl   bj&)  on  the  right 
side  of  Equation  (2-12)  the  line  Bq  =  (  1  +  J~2)  a0  (or  BQ  ■  aQ  line) 
is  calibrated  with  CO  linearly  such  that  the  points  Q  and  T,  correspond 
to  6)r  0  and  GO  =  1  respective ly9  and  the  length  QTj_  is  equal  to 
f~2   b,  .   Then  the  distance  from  Q  to  any  point  at   C0=(0j   en  the  line 
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Fig.  2-7  Determination  of  band  width  for  a  system 
with  one  zero. 


B0  =    (  l+/2)a0     will  represent    /Tb^j      and  it  is  obvious   that 

MW  s     /Tl^+jb,  CO1  j     which  is   the  quantity  on  the  right  side  of  Equation 

(2-12)   for      0)=  CO 

Suppose  the  (d(a     -  B,)  plot  is  as  shown  by  the  dotted  curve  in 
Fig.  2-7.   Noting  fehe  nature  of  the  iO(&i   -  Bj)  plot  that  it  always 
starts  at  the  point  (a  ,  0  ),   crosses  V0    at  a)  =  1  and  at  the  point 
where  (^  intersects  with  Bj_  s  a^  lines  it  is  easy  to  guess  or  sketch  the 
Cd  (a,  -  B-.)  plot. 

The  band  width  CO,  is  some  0)  such  that  Equation  (2-12)  holds. 

This  can  be  interpreted  in  Fig.  2-7  as  follows  s  Suppose  a  magnitude 

circle  C^     is  drawn  such  that  it  pass  through  the  point  W  at   6)sa)' 
on  the  B6  r  (1  +V~2)  a.0  line  and  suppose  it  reads  6i~    &>" 

on  the  co^  -  B^)  plot.   If   oo'ro)1'   then  Equation  (2-12)  holds  at 

frequency  equal  to  oS   ,  therefore  oifc  =  a)'  , 

To  determine  band  width  one  must  find  a  C.  such  that  it  crosses 
both  A)(a^  -  B^)  plot  and  Bq  -  (  1  ■+  fl   )aQ  at  the  same  frequency. 

For  a  Qi    if  0)V  w  then  the  two  frequencies  are  either  oJ  >o) 
or  0)'<  CO '  .    By  drawing  another  C^1   until  the  inequality  between 
£0   and  cO   is  reversed  from  the  previous  set  of  CO   and  (>0     the  bounds 
of  band  width  can  be  set  up. 

For  example,  for  the  configuration  of  curves  as  shown  in  Fig.  2-7 
1^2  "^  1   is  assumed. 

For  the  two  magnitude  circles  Caq  and  C*, 

Frequency  on  Bq  =  (1  +  /~2)ao  line      Frequency  on  <u)(a^  -  Bi)  plot 

cA0 coz  o < o)'  -  a)0 

CAi  d)'  ■  1  >  a)"  -  4>, 

in  which  the  inequality  sign  is  reversedj,  therefore  o<  c0.<  Cd^<  Cd,  <   I 
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It  is  obvious  in  Fig.  2-7  that  there  is  no  possible  band  width  for 
60  <  0)q   .   There  may  be  another  band  width  for  00  >  a)«  „  This  can  be  in- 
vestigated by  trying  some  magnitude  circles. 

C   in  Fig.  2-7  is  identical  to  CAb  in  Fig«  2-6.   If  h±   =  0  in  Fig0 
2-7,  then  C**    stays  on  C.q  for  all  00   9  which  is  the  case  for  a  system 
with  no  zeros  and  agrees  with  Fig,  2-6.  It  is  readily  seen  from  Figo 
2-7  that  the  effect  of  zero  and  forward  gain  of  a  system  on  band  width 
can  be  analyzed. 

For  two  systems,  one  for  no  zero  and  the  other  for  one  zero,  if  the 
characteristic  equations  are  the  samej,  then  band  width  of  the  first  system 
is  less  than  that  of  the  second.  This  is  confirmed  in  Fig.  2-7  as  for 
the  first  system  £0^  »  Odn     and  for  the  second  system   0),  >CL)q  - 

The  effect  of  varying  forward  gain  and  the  zero  on  band  width  for  a 
system  with  one  zero  is  seen  from  Fig.  2-7  by  writing  the  numerator  of  the 
closed  loop  transfer  function  as 

N(S)  s  bjS+  aQ  -  KS  -f  KZ 

Increasing  either  K  or  Z  is  to  expand  both  CAQ  and  C**  which  is  to 
increase  OOq  and  00^   resulting  in  increased  band  width.   Since  K  and  Z 
do  not  enter  into  any  coefficient  other  than  a,  or  a^  of  the  character- 
istic equation9  the  «7  curve  remains  unaffected  by  either  K  or  Z. 

Suppose  K  is  fixed  and  Z  is  variable.  Since  a^  and  b|  are  independ- 
ent of  Z,  the  Cd(ai   -  Bj_)  plot  and  the  lines  Bj  s  a-^   s  B|  :  a^  4  /I  bj_ 
remain  unchanged.  Varying  Z  moves  the  M  point  and  the  line  Bq  ■  aQ  ver- 
tically, thus  increasing  Z  expands  Caq  and  CA|  resulting  in  increased 
band  width. 
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Suppose  Z  is  fixed  and  K  is  variable g  since  a^  ,  aQ  and  b^  are 
linearly  varying  with  K,  all  of  the  lines  in  Fig.  2-7  and  the  0)(a]_-  Bj_) 
plot  are  moved.   The  M  point  is  moved  along  a  line  with  some  positive 
slope  and  Caq  and  CA]_  are  expanded  with  increasing  K  resulting  in  in- 
creased bandwidth. 

In  a  system  which  has  positive  coefficients  in  both  numerator  and 
denominator  polynomials  of  its  open  loop  transfer  functiona  b^  is  always 
less  than  a.  assuming  unity  feedback.   This  nutans  that  the  bandwidth  of 
a  system  with  one  zero  is  always  less  than  that  of  a  system  which  results 
from  replacing  b^  with  a, .   Therefore  for  a  system  whose  V^   and  M  point 
configuration  as  shown  in  Fig»  2-7  bandwidth  must  be 

600<  Cdb   <  0d2 
where  0)2  is  the  frequency  on  the  60 ( a,  -  Bj)  plot  at  which  the  6t)(a^  -  Bj) 
plot  and  magnitude  circle  C»2  intersect.   CA2  &s  drawn  such  that  it  passes 
through  T2  where  T2  is  located  on  the  Ba  :  (  1  +  /~2~  )a  line  at  the  dis- 
tance of  fl   a,  measured  from  Q. 

Above  inequality  is  useful  to  set  up  the  bounds  of  6),  when  b  is  a 
variable  and  when  it  is  required  to  locate  M  point  to  meet  some  required 
bandwidth  specification.   It  must  be  noted  that  the  bounds  determined  by 
the  above  inequality  become  poorer  as  the  difference  between  bj_  and  a, 
gets  larger. 

It  is  interesting  to  note  that  if  bj  =  a^   then  the  open  loop 
transfer  function  is  of  the  form,  assuming  unity  feedbacks 

&(s)  = a,s+a, 

^J      *"+a«.,s'-,+  ---  +a;s2 

which  is  a  type  2  system,  and  this  system  has  greater  bandwidth  than  a 


system  whose  open  loop  transfer  function  is 


££)=■ 


s*  ■+  a^.,s'n '  +  -  -  -  +  %  s*  +  («a  -  b, )  s 

For  a  system  with  two  zeros s  bandwidth  can  be  evaluated  in  a  simi- 
lar way  as  illustrated  in  Fig.  2-7.   One  modification  is  required  in 

Fig.  2-7  to  account  for  the  frequency  dependence  of  /"2~(aQ-  b£iO  ), 

2 
where  b£  is  the  coefficient  of  the  S   term  in  the  numerator  of  the 

closed  loop  transfer  function.   This  is  done  by  calibrating  the  line 

B-i  =  a-,  as  in  Fig.  2-4  to  represent  (&q  -  ^2^1   )  ^y  the  distance  from  M 

to  the  point  at  cOz    CO,  on  the  B*  =  a,  line.   To  account  for  the  factor 

J~~2   8  the  origin  of  calibration  or  the  Cd~   0  point  is  chosen  at  a  point 

on  the  B,  =  a,  line  below  the  B.  axis  so  that  it  is  apart  from  the  M 

point  by  a  distance  of  [2   aQ.   Then  the  line  B^  2  a^  Is  calibrated  with 

—     Z 
60  such  that  any  point  at  CO  z    CO,,  on  the  line  will  measure  /2  b2^>i 

by  the  distance  from  the   00=  0  point  (originating  point)  t©  the  point 

where   00  =  ul  . 

__  I         2  2  I 

For  a  system  with  3  zeross  the  quantity  7  2  (aQ-b2&)  ) -i-j'COCbi^bjfjO  ) 

which  is  the  right  side  of  the  bandwidth  defining  equation  can  be  repre- 
sented in  a  similar  way  as  in  Fig.  2-5.   The  only  difference  is  to  ex- 
pand the  translated  d)(b^  -  B,J  plot  radially  from  the  M  point  by  factor 

of  fi  . 

Then  for  both  cases  of  the  above „  bandwidth  can  be  evaluated  by 
applying  the  same  analysis  as  for  the  system  with  one  zero  in  Fig.  2-7. 
The  following  numerical  examples  illustrate  the  arguments  developed  so 
far. 
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Example  2.1 

Obtain  frequency  response  (magnitude)  of  the  closed  loop  system 

i®=  sw^s+fr-  (2'13> 

The  system  represented  by  Equation  (2-13)  is  a  third  order  system. 

Since  it  has  no  zeros  the  bandwidth  will  be  less  than  the  frequency  at 

the  point  where  the  f~0   curve  and  Bq  =  (1  +  J~2   )aQ  line  intersect.   If 

the  frequency  response  is  required  up  to  OJ  r  oO,  s  the  portion  of  the 

b 

I  o     curve  needed  is  below  the  line  B_  =  (It/  2  )a~  =  9.65. 

In  Fig.  2-8  the  frequency  response  of  the  system  (2-13)  is  repre- 
sented on  the  Bq  vs  B,  plane. 

From  the  characteristic  equation  of  the  system  (2-13)  the  \0   curve 
is  defined  by 

2.  (2-14) 

which  is  a  straight  line  with  slope  of  3  passing  through  the  origin,, 
The  magnitude  of  the  frequency  response  is  defined  by 

a6       _      n 

where  B  and  Bq  are  defined  by  Equation  (2-14). 
The  point  M  (2,  4)  is  located,  and  the  lines  B^  =  a^  =  28  Bq  =  &q  -   4 

and  Bq  =  (1  +  f2   )aQ  =  9.65  are  drawn  as  shown.  The  0)  (a^  -  B^  plot  is 

constructed  as  illustrated  in  Section  2-1. 

The  &}(aj.  -  B^)  plot  starts  from  the  point  (2,  0)  and  crosses  the 

f|j  curve  at  0)  =  1  and  the  point  (2S  6)  at  which  the  ^  curve  and  the 
line  B^  =  2  intersect. 

For  0  <  00  <  1  and  for  W>1.42  the  fl)(aj_  -  Bj_)  plot  lies  to  the 
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right  of  the  (^  curve  and  for  1.0  <  CO  <  1.42  t©  the  left.  The  point  at 

0)  s    .5  on  the  C0(al   -  B^  plot  is  located  on  a  horizontal  line  passing 
through  the  point  at  6).-  .5  on  the  ^J  curve  and  the  distance  from  the 
line  B1  =  2  to  the  point  at  0)  =  .5  on  the  0)(a1  -  B1)  plot  is  one-half 
of  the  distance  from  the  line  B,  =  2  to  the  point  at  G)s  .5  on  the  1~J 
curve  so  that  the  distance  from  the  Bx  -  2  line  to  the  point  at  60s  .5 
on  the  a)(ai  -  BL)  plot  is  equal  to  CO  (2  -  B1(C0  ))  |        which  is  the 
value  of  imaginary  term  at  60  >  .5  in  the  denominator  of  the  equation 
defining  the  magnitude  of  the  frequency  response „ 

Some  points  on  the  0>(ai  -  Bi)  plot  for  some  pertinent  values  of 
frequencies  are  determined  likewise  and  those  points  are  connected  by  a 
smooth  curve  to  give  the  C0(a^  -  B|)  plot  as  shown.   Thus  each  pair  of 
two  points  on  a  horizontal  line,  one  on  the  1^  curve  and  the  other  on  the 
60  (a,  -  Bj_)  plot  has  one  value  of  frequency  and  the  distance  from  the 
point  M  (2,  4)  to  any  point  at  frequency  ft)  on  the  C0(&i   -  Bi)  plot  is 
the  value  of  the  denominator  of  Equation  (2-15). 

Magnitude  circles  are  constructed  as  shown.  Noting  that  the  numer- 
ator of  Equation  (2-15)  is  constant  at  aQ  -  U.,  it  is  seen  that  the  ratio 
of  aQ  s  4  to  radius  r  of  a  magnitude  circle  which  passes  through  a  point 
at  frequency  of  ft)  on  the  0)(aj  ~   Bj_)  plot  is  the  magnitude  A  (ft)). 

At  CO  =  03  the  radius  of  magnitude  circle  is  4  which  gives  A  =  1. 
As  the  frequency  increases  from  zero,  the  radius  of  the  magnitude  circle 
decreases^  in  turn  increasing  A.   At  CO  *   1.18  magnitude  circle  is 

at  to  the  a)(ax  -  Bt)  plot  and  has  radius  r  s  . 75  which  implies 
peak   Mpt0  =  4/.  75  =  5.33 
ft)  =  1.18 
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For  frequencies  from  CO  =  1.18  on,  the  radii  of  magnitude  circles  in- 
crease s  in  turn  decreasing  A. 

A  magnitude  circle  with  radius  of  /2  afi  :  5.65  crosses  the 
CO  (ajL  -  B^  plot  at  Cd  -   lo76  which  implies  the  bandwidth  is 

Cd^z   1.76 
The  frequency  response  up  to  Cd  z   1.8  is  thus  represented  on  the  B«  vs 
B,  plane. 
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Example  2.2 

Given  an  open  loop  transfer  function  of  a  unity  feedback  servo  as 

obtain  the  magnitude  vs  frequency  plot  of  the  frequency  response. 
The  closed  loop  transfer  function  is 

Since  Equation  (2-16)  deals  with  large  numbers ,  it  is  not  convenient  to 
work  on  the  Bq  vs  B^  plane  associated  with  Equation  (2-16). 
By  using  a  transformation  or  frequency  scaling 
S  =  100s 
Equation  (2-16)  transforms  to 

_C_(    v   /?./£     -f  4-»g ___ 

K^S)~   S3+/6.'7/5a  +  /2.z^7s  +  ^.^  (2-17) 

Then  the  equation  defining  frequency  response  is 
_C_,.       v   J  g.gg  +jJ2.t  C6t\ 

and  B~  are  Mitrovic" s  variabl 
for  %  =   0   associated  with  the  system  (2-17)  and  £0   is  re- 
lated with  <x),  the  frequency  associated  with  the  system  (2-16) 

by 

£0=  100  Cdt 

In  Fig.  2-9  the  \0  curve,  0\(a,.  -  B,  )  plot  and  magnitude  circles 
representing  the  numerator  of  Equation  (2-18)  are  constructed. 

The  frequency  dependence  of  the  imaginary  term  in  the  numerator  of 
Equation  (2-18)  is  represented  on  the  BQ 
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B^t  :  a  line.  (Note  that  the  frequency  dependence  of  the  imaginary 
term  in  the  numerator  of  Equation  (2-18)  is  represented  on  a  line  dif- 
ferent from  that  illustrated  in  Section  2-1.) 

As  illustrated  in  Section  2-2,  the  distance  from  the  point 
M  (12.267,   4.48)  to  a  point  at  frequency  C0±  on  the  BQt  s   2aQt  s  8.96 
line  is  equal  to   |4. 48  +j  12.267  CO  A 

The  magnitude  circles  are  drawn  so  that  at  any  frequency  the  numer- 
ator and  denominator  of  Equation  (2-18)  can  be  compared. 

Since  the  M  point  is  located  at  the  point  where  the  coordinates  are 
not  convenient  numbers,  reading  off  the  radii  of  magnitude  circles  is 
not  convenient. 

This  can  be  simplified  by  providing  a  scale  which  is  calibrated 
with  the  same  scale  of  coordinate  axes  so  that  the  radius  of  any  magni- 
tude circle  is  read  off  easily.   This  is  shown  in  Fig.  2-9,  for  example , 
at  C0t  s  .5,  the  value  of  |N(jcO  )|  v       is  read  to  be  7.5  on  the 
scale  provided  on  the  line  B-.t  =  12.267  by  drawing  a  magnitude  circle 
passing  through  the  point  at  0)fc  =  .5  on  the  line  Bqc  s  2a~t  and  simi- 
larly for  the  denominator   (aQt  -  Bqc)  -Hi |  (^(a,  -  B.  )  I  \Q)_     ^    s  6.0S  thus 
the  following  table  is  constructed; 
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Example  2.3 

Suppose  for  the  system  in  Example  2.2  only  the  bandwidth  is  of  in- 
terest. The  equation 

|^^-Bot(^t^  +}&*  02'a7'^(^it)\'/z\^^V/;ij6)t\       (2-19) 
defines  the  bandwidth  ft).   of  the  system  in  Equation  (2-17). 

By  drawing  a  magnitude  circle  with  radius  equal  to  1. 414(4. 48+j 12.1) 
which  is  the  quantity  on  the  right  side  of  the  equation  above  at  a)   =  1, 
and  noting  the  frequencies  in  the  neighborhood  of  the  point  on  the   P0 
curve  where  the  circle  intersects  with  the  \0   curve,  and  also  noting  the 
relative  location  of  the  ^t(ait  *  B-. ,.)  plot  with  respect  to  the  V0 
curve,  some  information  can  be  obtained  about  the  bandwidth.   This  is 
illustrated  with  Fig.  2-11. 

In  Fig.  2-11,  the  h   curve  is  drawn  for  the  system  (2-17),  the  work- 
ing point  M  (12.27,  4.48)  is  located,  and  B1(.  =  au  -  12.27  is  laid  off 
the  same  way  as  in  Fig.  2-9.  The  line  Egt  =  (1 f  1.414)a0t  =  10.8  is 
drawn. 

From  the  intersection  point  of  BQ.  =  10.8  and  B^fc  «  12,27  lines  to 
the  left,  a  distance  of  1.414  b^t  -  17.1  is  measured  off  on  the  Bgt  = 
10.8  line  and  labeled  Od      -   l>0.        The  radius  of  a  magnitude  circle  pass- 
ing through  the  U)     =  1.0  point  on  the  BQt  r  10.8  line  is  equal  to 

1.414(4.48  +j  12.1)  I  r  1.414  N(j  1)  I.   This  circle  denoted  by  CAl  crosses 
the  \0    curve  between  the  points  at  G).    =   1.0  and  i^t   =  1.1. 

It  is  easy  to  predict  that  0)      <  1  is  not  possible.   If  it  is  sus- 
bt 

pected  that  CO,  <  1  is  impossible,  it  can  be  checked  as  follows? 

The  check  is  based  on  the  fact  that,  for  the  bandwidth,,  the  magni- 
tude circle  must  cross  the  ti>t(.&it  "   Bqc)  plot  and  1.414  b^t  ^t  scale 
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represented  on  the  line  BQt  =  (l-t,/~2  )a,Qt   z   10.8  at  the  same  frequency, 
and  the  fact  that  if  any  magnitude  circle  is  drawn  and  if  it  crosses  the 
&)t(alt  -  B0t)  plot  and  1.414blt^scale  (BQt  -  10.8  line)  at  o)t  s  o)t" 
and  L0t  =  Oj'   respectively  with  either  OK  00"  or  oj>  a)  ,  then  for 
existence  of  a  bandwidth  frequency,  there  must  be  another  magnitude  circle 
to  give  a  new  set  of  CO   and  ck>      with  inversed  inequality  from  the  previous 
set  of  0)    and  0)"  . 
magnitude 
0)£  >  J    and    60t'  =  i         then      to  '  <  u)'t 
Since  C.^  is  drawn  with  cu'r  1  and  it  crosses 
lx)t   ■  1.0  and  0)  =1.1  and  1.0<O\<l.l.   This  is  reasoned  without 
drawing  any  portion  of  the  ^t(ait  -  B,  )  plot« 

Draw  a  magnitude  circle  with  CO  z   0.   This  circle  is  of  radius  equal 
to  [2   aQt  tangent  to  the  line  Bqj.  =  (1  +  fl   )aQt  z   10.8  and  represents 
fT|N(jO)|  =  [~2   aQt  by  its  radius.   Since  60t(au  -  Bu)  plot  ap- 
proaches to  the  point  at  a)  =  1  on  the  1^  curve  from  the  point  (12.27,  0), 
the  circle  must  cross  the  Od    (a*.    "  Blfc)  plot  at  some  frequency  greater 
than  zero  to  give  an  inequality 

0)^=0  }      COt  >  0  then     o)t'  <  0)^' 

The  inequality  sign  is  not  reversed  from  the  previous  set  of  u>t  and  a)^', 
and  this  implies  there  might  be  no  frequency  such  that,  for  0  <  W.  < 1  , 
a  magnitude  circle  crosses  the  (jO+(a,  t  -  B,  .)  plot  and  the  fl   b,  fc  Qt 
scale  at  the  same  frequencies;  in  other  words,  60  -  00,.=  uOht  may  not 
be  possible  for  the  frequency  range  from  zero  to  1.0.   However,  there 
may  be  some  frequency  between  0  and  1.0  at  which  the  inequality  is  re- 
versed since  the  interval  of  6)+   from  0  to  1.0  is  equivalent  to  a  rela- 
tively large  interval  of  original  system  frequency  00  -    100  &)  , 
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Locate  the  point  at  0)     .-  .5  on  the  ^  (a.   -  B.  )  plot  which  is 
half  way  from  the  line  Blt  =   a^t  =   12.27  to  the  point  at  u)     =  .  5  on  the 
Yq     curve.   This  point  is  located  inside  the  magnitude  circle  correspond- 
ing to  (jSz   0,  therefore,  no  bandwidth  is  possible  for  Od .   ^  .5. 

For  the  magnitude  circles  corresponding  to  „  7  <  a) '  <:  1,  the  inequal- 
ity 03  <£0  is  not  changed.   This  can  be  checked  by  comparing  the  dis- 
tances from  the  M  point  to  the  'j  curve  and  the  [2   bjt  cdt   scale  and 
noting  that  the  <^t(alt  _  B1(.)  plot  lies  to  the  right  of  the  \0   curve  for 
.7  <  C0t<  1.0. 

For  example,  a  magnitude  circle  for  0)t  z    . 7  will  cross  the  ^  curve 

at  a  point  where  cd     is  greater  than  .8  to  give  id       .8  and  o)t  <  tD^ 

Thus  for  the  bandwidth,  .  7  <-.  0J,  <  1.0  is  not  possible, 
bt 

The  only  part  left  to  check  is  the  interval  .5  <60t<  ,7.   Knowing 

that 

for  CO^z'   .5  to'  <  0)'      and 

for  d£=  .7  ti)'<(0" 

it  is  quite  sure  that  for  the  interval  .5  <  id   <  .7  the  same  inequality 
will  be  maintained.   If  a  magnitude  circle  for  fl)t  r  „6  which  is  the  mid 
point  of  the  interval  of  the  frequency  from  .5  to  .7  is  drawn,  this  cir- 
cle will  pass  the  \n   curve  at  the  point  where  0)      =  .6  implying  the 
point  at  C0t   =  .6  on  the  ^t(ait  -  Blt)  plot  will  be  inside  the  circle 
in  turn  Od't<(Jdt   which  is  the  same  inequality  as  the  cases  of  (rt=    .5  and 
60^=  .7.   Thus  the  same  inequality  is  maintained  for  Cdt   from  0  to  1.0  and 
t£>t>t<  1  is  not  possible. 

The  tests  made  above  are  carried  out  without  plotting  any  portion  of 
the  Ot(ait  -  Bit)  plot  except  the  point  at  o)t  =  .5  on  the  0dA&it-  B^t) 
plot  is  located.   The  tests  can  be  carried  out  with  a  sketch  of  the 
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60t(alt  -  B1(.)  plot  and  the  sketch  can  be  made  simply  with  reasonable 
accuracy  by  noting  that  the  horizontal  distance  from  the  line  Bjt  r  alt 
=  12.27  to  the  COt(alt  -  Blt)  plot  is  a  fraction  of  that  to  the  ^  curve 
and  the  fraction  is  the  value  of  0)  . 

For  an  interval  of  frequencies  (CO,,  ^) ,  if  the  inequality 
0)'  <  O)"  (err      60 '  >  CL)") 
holds  at  both  ends  ft),  and  Clk  then  the  inequality  is  maintained  through- 
out the  interval  unless  the  >0   curve  or  the  k>(a,   -  Bit^  plot  is  oddly 
shaped  in  the  interval,  or  in  other  words,  unless  the  rates  of  change  of 
the  distance  from  the  M  point  to  the  a)t(ait  -  Bj_t)  plot  is  not  fluctuat- 
ing in  the  interval. 

In  this  example,  if  the  COfc(a.   -  B,  )  plot  is  sketched t    then  one 
can  see  that  the  rate  of  change  of  distance  from  the  M  point  to  the 
LOt(a^t  -  B^t)  plot  is  steadily  increasing  except  in  the  region  of  small 
frequencies. 

Knowing  that  0),     >1,0,  the  next  thing  to  do  is  to  find  another 
magnitude  circle  to  insure  an  upper  bound  of  t)   as  illustrated  in  Sec- 
tion 2-3. 

Noting  that  C»i  crosses  the  Iq  curve  at  the  point  where  the  fre- 
quency is  slightly  less  than  1.1  and  also  noting  that  the  ^t(ait  -  Bit) 
plot  lies  to  the  left  of  the  \q   curve  for  CO    >  !?>  it  is  predicted  that 
C0bt  will  not  be  very  far  from  CO  z   1.1.   Therefore  it  is  natural  to 
try  a  magnitude  circle  which  passes  through  CiL  =  1.1  or  ft).  =  1.2  on 
the  f2   b^tCt)t  scale.   The  magnitude  circle  which  passes  the  point  at 
(X)t  =  1.2  on  the  line  Bot  r  (1  +  /2  )aot  is  tried  and  denoted  by  CA2° 
The  circle  CA2  represents  the  right  side  of  Equation  (2-19)  at  0)t-1.2o 
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It  passes  the  \0   curve  at  a  point  where  the  frequency  is  less  than  1.2 
which  implies  it  will  pass  a  point  on  the  &)t(a^t  -  B]_t)  plot  where  fre- 
quency is  less  than  1.2,  thus  id      =  1.2  clearly  indicates  an  upper  bound 
of  bandwidth  and 

To  determine  OX   precisely,  the  ^t^lt  '  Blt^  Plot  for  1.0<4).< 
1.2  is  needed. 

A  magnitude  circle  C.u  which  passes  through  the  point  at  o)      :  1.1 
on  the  <f2   b^t  Ci)t  scale  crosses  the  ^t(a^t  -  B^t)  plot  at  cOt   =  1.1  and 
bandwidth  is  determined  to  be 

Ht= 1-1 

This  agrees  with  Fig.  2-10. 

Since  Fig.  2-10  or  2.11  is  for  the  system  represented  by  Equation 
(2-17)  which  is  related  to  the  original  system  by  S  =  100s  ,  the  original 
system  bandwidth  OX  is 

(A    =  100  60bt  =  110 

By  substituting  s  =j  6_)  =   1.1  and  S 

(2-17)  and  (2-16)  respectively,  the  magnitudes  calculate  to 

1 1(0 1  .     1=  -707 


&L.\- 
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2-4.  Locus  of  M  points  on  the  B^  vs  B-.  plane  for  constant  bandwidth. 

In  Section  2-3,  bandwidth  of  a  given  system  is  determined.   For  a 

given  system  with  fixed  M  point  which  is  determined  by  the  coefficients 

a,  and  a^,  bandwidth  is  determined  by  varying  a)   until  some  frequency 

(or  frequencies)  Cx)  =  <j  is  found  to  satisfy  Equation  (2-10). 
b 

Suppose,  conversely,  CO   in  Equation  (2-10)  is  fixed  at  bandwidth 
60b,  then  there  must  be  some  a-^  and  ag  to  satisfy  Equation  (2-10),  and 
the  pair  (a  ,  ag)  defines  the  M  point  on  the  Bq  vs  B\   plane  for  a  system 
to  have  bandwidth  equal  to  CO..      If  all  of  such  pairs  (a,  ,  ag)  satisfying 
Equation  (2-10)  are  plotted  on  the  Bq  vs  Bi  plane,  then  a  curve  which 
represents  the  locus  of  M  points  on  the  Bq  vs  B^  plane  for  fixed  band- 
width 0)l  is  obtained.   This  curve  is  called  the  constant  bandwidth  curve 
for  bandwidth  equal  to  £0,  or  merely  constant  0)^   curve. 

If  a  constant  CO,  curve  is  constructed  on  the  Bq  vs  Bj^  plane,  and  if 
the  M  point  is  chosen  at  any  point  on  the  curve,  then  the  resulting  sys- 
tem will  have  bandwidth  equal  to  COv. 

Since  the  synthesis  of  a  feedback  control  system  by  Mitrovic's 
method  is  to  choose  a  suitable  M  point  to  guarantee  the  required  system 
performance  specifications,  the  constant  C0^   curve  will  provide  informa- 
tion in  choosing  the  M  point  to  meet  specified  bandwidth  requirement. 

By  fixing  6J=  0),    in  Equation  (2-10),  the  constant  (0,    curve  is 
defined  by 

|(V^^V+j^(S-B,«0^  |  -ft  |  N  (j  C0t)|  (2-20) 

where  B,  (60b)   and  B~  ( CO.)   are  the  coordinates  of  the 
points  on  "  curve  at  CO  =  £0^. 

On  the  left  side  of  Equation  (2-20),  Bi(6)b)  and  Bq(OL)  are  fixed  and 
those  values  can  be  read  off  on  the  l7  curve. 
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On  the  right  side  of  Equation  (2-20),  noting  that  N(S)  is  the  numer- 
ator polynomial  of  the  closed  loop  transfer  function,   |n(  /a)b)j  contains 
aQ  which  is  treated  as  a  variable  on  the  constant  6)^  curve.  Moreover, 
if  N(S)  is  of  order  other  than  zero,  it  may  contain  b,  which  is  the  co- 
efficient of  the  S   term  in  N(S),  and  b,  is  not  independent  of  a,  which 
is  the  other  variable  on  the  constant  03,  curve.   Thus,  even  though 
Equation  (2-20)  is  a  perfectly  general  expression  for  defining  constant 
&)|j  curves,  it  does  not  clearly  indicate  the  geometrical  shape  of  the 
constant  U),  curve  on  the  B_  vs  EL  plane. 

In  the  following,  the  cases  of  zero  order,  first  order  and  second 
order  of  N(S)  are  considered. 
N(S)  of  zero  order. 

When  a  closed  loop  transfer  function  has  no  zeros,  N(S)  r  a~  and 
Equation  (2-20)  becomes 

|  ( VBoK))  +j  *>b  {'  a-  b,  wi)  |  =  £  do 

by  removing  the  absolute  sign,  and  manipulating 

( aD+ B0(oV>y-  <  (4 ,  -  B,(*vf  -  *  BX^  (2-2i) 

Since  6X,  BQ(0)b),  B,  (  6.),)  are  fixed,  Equation  (2-21)  indicates  a  hyper- 
bola with  following  geometrical  properties  on  the  Bq  vs  B-^  planes 

the  center  of  symmetry  at    (  B/a)^  ,  -Bd0%)) 

the  vertices  at    (  B^t) ,  (-l±Jz)B0(titS) 

major  axis  is    j31  =  J3,  (o)^ 
The  two  branches  of  hyperbola  (2-21)  lie  symmetrically  with  respect  to 
the  Bq  =  "Bo( ^k)  line.   The  upper  branch  lies  above  the  B^  axis  with 
minimum  point  at  (B^(a)b),  .414  Bo(OJt>))  which  is  one  of  the  vertices. 


For  a  stable  system,  only  the  portion  of  constant  fr)b  curve  which  lies 
in  the  stable  region  is  of  interest. 

It  is  interesting  to  note  that  the  locus  of  the  vertex  of  the  branch 
of  the  hyperbola  which  lie  in  the  first  quadrant  is  similar  to  the   [7 
curve,  and  the  locus  of  the  center  of  the  hyperbola  is  the  mirror  image 
of  the  \0   curve. 

Those  can  be  seen  as  follows.  As  shown  previously  the  vertex  of  the 
upper  branch  of  the  hyperbola  is  at  (B^CO^),  ,414B0(O)b))  and  the  center 
point  is  at  (B^a^),  -BQ(  ft)fe)). 

Considering  00,    as  a  continuous  positive  variable,  then  the  locus 
of  the  vertex  point  is  defined  by 

X  =  B±(0)) 

Y  =  .414  B0(ti) 

and  in  like  manner  the  locus  of  the  center  point  is  defined  by 
X  :  B^U)) 

Y  =  -BQ(C0) 

where  X  and  Y  axes  are  identical  to  B^  and  Bq  axes  respectively. 
The  two  sets  of  equations  can  be  compared  with  the  [7  curve  defined  by 
Bl   z   BjCft)) 
BQ  =  B0(CO) 
Consider  the  consistency  of  principles  involved  in  constant  0).  curve 
and  in  determining  bandwidth  described  in  Section  2-3. 

In  Fig.  2-12,  suppose  a  constant  Q),    curve  is  plotted  from  Equation 
(2-21)  assuming  the  \^   curve  as  given.   The  point  P  is  at  OJ  z    0},    on 
the  I  a  curve.   The  constant  &),    curve  is  denoted  by  C   .   The  vertical 
axis  denoted  by  Y  is  on  the  B,  =  Bi(^b^  ^ne   which  passes  through  the 
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Fig.  2-12  Showing  consistency  of  Section  2-3  and 
Section  2-4  about  bandwidth. 
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point  P  and  V  is  the  vertex.   The  horizontal  axis  and  the  lower  branch  of 
the  hyperbola  which  lie  below  the  B^  axis  are  not  shown. 

The  portion  of  the  constant  CO.  curve  which  is  of  interest  for  a 
stable  system  is  the  portion  of  C^.,  under  the  ^  curve. 

Suppose  an  M  point  is  chosen  arbitrarily  on  C  cQh   as  shown.   Then  a 
Od  (a,  -  B, )  plot  results  as  described  in  the  previous  sections  and  a 
magnitude  circle  C»  with  radius  equal  to  fl   a  must  pass  through  the 
point  Q  at  which  (X>-    (x),    on  the  &  (a^  -  B^)  plot  to  have  a  consistency 
for  the  principles  involved  in  the  constant  CO.    curve  and  determining 
bandwidth  as  described  in  Section  2-3.   This  is  to  show  that 

Locating  a  point  R  at  the  intersection  point  of  the  B^  =  a,  line  and 
B0  =  B0(0)b)  line, 

1/2. 

But  (a,  ,  a«)  satisfies  Equation  (2-21)  and  it  rearranges  to 

and         MQ=/^a0 

N(S)  of  first  order. 

When  a  closed  loop  system  has  one  zero 
N(S)  =b,S+  20 
For  unity  feedback  control  system,  b-^  enters  into  a^  and  a^  and  b-^  are 
not  independent  of  each  other. 

Suppose  an  open  loop  transfer  function  of  a  unity  feedback  control 
system  as 
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where  b^  ,  a^  and  a.  are  variable  and   cLs  are  fixed  constants, 

few-- 


or  r<s>=  g%a,g^...ta,stat  (2"22) 

where  a     ,    -   d„    ,    ,      ...      a0   =   d0     and     a,  :d,  +  bn 
n-1         n-1  2d,  111 

N(S)   can  be  written  in   the   form  of 

N(S)    =    (aL   -   dL)   S+aQ 

With  this  N(S),  Equation  (2-20)  is  manipulated  to 

[  VBa((0^f  +  «b[a,+(^iKV^il]  =  P(<M,)  (2-23) 

where   D(0)b,  d.)  =  2[p£(0)b)  +  0)f  (B,(fl)^-^fJ  -+-df^l 

Equation  (2-23)  defines  the  constant  oV  curve  for  any  system  whose  closed 

D 

loop  transfer  function  has  one  zero. 

Noting  that  for  D(<Db,d1)>0,  Equation  (2-23)  is  an  ellipse  centered 
at   (2d^  -  B^(0>b)  ,  -B_(  <^>,) )  with  semi-axes  of  lengths  equal  to 


^V'>    and  V^ 


(^fi^i)   •    (Horizontal  and  vertical  axes  respectively. 

For  each  d-,  ,  as  CO.    varies  there  results  a  family  of  ellipses. 
•L        b 

It  is  seen  in  Section  2-3  that  for  a  system  represented  by  Equation 
(2-22),  larger  b^  results  in  larger  bandwidth. 

For  any  unity  feedback  system  assuming  all  positive  coefficients  in 
its  open  loop  transfer  function,  b-^  is  always  less  than  a-,  and  its  upper 
limit  is  ai  which  is  the  case  of  a  type  2  system,  whereas  its  lower  limit 
is  zero  which  is  the  case  of  a  system  with  no  zeros. 

For  the  upper  limit  of  b.  =  a,  corresponds  to  d,  =  0,  and  letting 
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6^   -   0  in  Equation  (2-23),  the  constant 

(a0+  B0(o>of+ &)£(*!  +B,C^f-^(B0>l,H<3lH^     (2„24) 

Equation  (2-24)  can  be  applied  for  a  type  2  system  with  one  zero  in  its 
closed  loop  transfer  function.   For  any  system  with  b, <  a,   or  d^>  0  , 
if  the  constant  0),  curve  is  plotted  by  Equation  (2-24),  and  if  the  M 
point  is  located  on  the  curve,  then  the  system  will  have  bandwidth  less 
than  C0b,  thus  Equation  (2-24)  defines  a  locus  of  M  points  on  the 
Bq  vs  B^  plane  for  some  upper  bound  of  bandwidth  for  a  system  whose 
closed  loop  transfer  function  is  represented  in  the  form  of  Equation 
(2-22). 

This  will  be  illustrated  with  a  specific  system 

with  associated  constant  (d,    curves  as  shown  in  Fig.  2-13. 

For  the  system  (2-25)  the  l0  curve  is  defined  by 
BL  =  CO2" 
BQ  =  CO  -  0) 
and  this  curve  is  plotted  in  Fig.  2-13  as  shown. 

If  the  system  had  no  zero,  then  constant  Cd,    curves  would  be  defined 
by  Equation  (2-21)  with  B,  and  Bq  defined  as  above.   For  this  case,  con- 
stant Cl),    curves  are  constructed  for  (A    :  .1  to  6)^  =  .5  as  shown  in 
Fig.  2-13.   If  M  point  is  chosen  at  any  point  on  C^  which  is  a  constant 
6)b  curve  for  60  r  .5  then  the  system 


will  have  bandwidth  equal  to  .5. 
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(2-26) 


The  dotted  curve  which  is  a  parabola  is  the  locus  of  the  vertices 
of  the  constant  (Db  curves  for  the  system  (2-26). 

A  constant  &)b  curve  of  Cd^   =  i„o  for  a  system 

which  is  obtained  by  letting  b,  =  a  in  Equation  (2-25)  is  plotted  by 
Equation  (2-24)  and  denoted  by  C.  0  in  Fig.  2-13.   If  the  axes  scales 
are  made  the  same,  the  C^  q  would  be  an  arc  of  a  circle,  since,  for 
03b  =  1,  BQ(1)  ;  0,  BL(1)  =  1   then  Equation  (2-24)  is 

a*  +  (al  +  if  «z 

which  is  a  circle  centered  at  (-1,  0)  with  radius  equal  to J2 . 

If  the  M  point  is  chosen  at  any  point  on  C,  Q  then  the  system 
(2-27)  will  have  bandwidth  equal  to  1.0. 

Suppose  the  point  M(.407,  .089)  is  chosen  for  both  systems  (2-26) 
and  (2-27).   Since  the  point  M  (.407,  .089)  is  on  both  C>5  and  Cl   0  » 
the  systems  (2-26)  and  (2-27)  have  bandwidth  0>b  =  .5  and  6J>b  =  1.0  re- 
spectively. 

For  system  (2-26) 


m** 


.ofl 


S    +Sd+S*+,*>7£  +  ,029  -Sb.s 


=  .J07 


For  system  (2-27) 


Kji.»)  = 


.uo7S  +  .0*9 


S*  +  S3  +  £*  +  •  *b7£ '  +  •  &#9 


S*jl* 


•1°1 


Then  the  point  M  (.407,  .089)  for  the  system  (2-25)  guarantees  that 
0.5  <    CO,  <   1.0 


for 


<  bl  <  al 
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N(S)  of  second  order  and  up. 

When  a  closed  loop  system  has  2  or  more  zeros 

N(s>  z  aD+  b.s'-f  t?zgz+  -  - ' 

If  the  coefficients  b^S   for  k^2  are  fixed  (in  order  for  a  sys- 
tem to  be  analyzed  or  synthesized  by  Mitrovic's  method  on  the  Bq  vs  B^ 
plane,  the  coefficients  a^S   for  k  ■>  2  must  be  determined  and  likewise 
for  b^S  ),  then  by  proceeding  in  the  same  uay  as  the  case  of  N(S)  =  b^S 
■+-  aQ,   the  equation  defining  the  constant  ,;   curve,  which  is  similar  to 
Equation  (2-23)  with  due  modification  to  account  for  the  coefficients  b2, 
b  ,  .  .  .  will  be  obtained  and  likewise  for  the  equation  defining  constant 
CO^   curve  of  certain  upper  bound  for  bandwidth. 

It  must  be  noted  that  for  a  high-velocity  constant  system  b,  is  very 
close  to  a,  since  the  velocity  constant  Ky  is  related  to  a*,  aQ  and  b^  by 

K'~    4,-b, 

Therefore,  for  a  high  K  system,  the  constant  60-^   curve  obtained  by  let- 
ting b^  r  a^  will  be  a  good  approximation  to  the  actual  bandwidth. 

Also  in  lag  compensation,   a-^  =  b,   and  as  shown  in  references  1 
and  2,  lag  compensation  problems  are  solved  by  letting  a,  =   b,  ,  conse- 
quently, the  constant  60,  curve  obtained  from  a^  =  b-^  will  be  good  in 
lag  compensation  problems. 
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2-5.  Bandwidth  of  third  order  systems. 

It  is  shown  in  references  1  and  2  that  the  equations  of  third  order 
systems  can  be  normalized  to  give  a  universal  Mitrovic's  chart  for  third 
order  systems,  and  this  chart  is  constructed  in  references  1  and  2. 

In  the  normalized  chart,  \Q    is  a  straight  line  extending  from  the 
origin  to  infinity  with  slope  of  unity  on  the  BQn  vs  B^n  plane,  and  the 
frequency  CO   in  the  normalized  system  is  related  to  the  frequency  CO  in 
the  unnormalized  system  by 
U>  =  k  cotY1 

where  fc    is  a  constant,  the  normalizing  factor,  and  sub- 
script n  denote  normalized  quantity. 
Assume  a  third  of der  system  with  no  zeros  in  its  closed  loop  trans- 
fer function: 


it(s)=^r- 


It  normalizes  to,  by  transformation  S  z   a£S 


(2-28) 


where  dIT)  =  <2,  /<2Z 

?  (2-29) 

%>/?  -  2o/d* 

Equation  (2-29)  defines  the  relation  between  coefficients  of  the 
normalized  and  unnormalized  systems,  consequently  the  Mitrovic's  variables 
are  related  by  Equation  (2-29)  between  normalized  and  unnormalized  sys- 
tems, 

namely      B|71=B,/<2^ 

(2-30) 
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Applying  Equation  (2-2i)  to  the  system  represented  by  Equation 
(2-28),  the  constant  0),  curve  is  defined  by 

(a^  +iyfiWf  +  %n  fai*  -B|7.(ft)^f  =  a  B^  (<*W  (2-31) 

where   U)   is  the  bandwidth  of  normalized  system  and  is  re- 
lated to  the  unnormalized  bandwidth  00,  by  a),    :  a„(J,  . 

b  J        b    2  bn 

Bin  and  BQn  are  Mitrovic's  variables  for  the  ^7  curve  in  the 

normalized  third  order  chart,  and  are  related  with  unnormalized 

variables  B^  and  Bq  by  Equation  (2-30). 

Equation  (2-31)  defines  a  family  of  constant  to,   curves  on  the 

bn 

BQn  vs  B^n  plane  if  a  set  of  60^n  is  supplied. 

If  this  family  of  constant  GO,   curves  are  superposed  on  a  third 
order  chart,  then  it  will  provide  additional  information  about  the  third 
order  system. 

Fig.  2-14  shows  the  family  of  constant  bandwidth  curves  for  third 
order  systems  whose  closed  loop  transfer  function  is  represented  by  Equa- 
tion (2-28)  on  the  normalized  Bq  vs  B^  plane.   Note  that  the  portions  of 
constant  bandwidth  curves  defined  by  Equation  (2-31)  which  lie  in  the 
unstable  region  are  omitted. 

Fig.  2-14,  the  constant  bandwidth  chart  for  third  order  systems,  can 
be  superposed  on  Fig.  2-15  which  is  the  third  order  Mitrovic's  chart  con- 
structed in  references  1  and  2  to  provide  additional  information  concern- 
ing system  bandwidth  in  both  analysis  and  synthesis  of  third  order  sys- 
tems. 

To  determine  system  bandwidth  is  to  locate  the  M  point  of  given 
system  on  the  constant  bandwidth  chart  and  consult  with  the  constant 
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r*~U/ 


;Pig»  2-14     Constant  bandwidth  curves 

(3rd  order  systems-normalized, for  b..  =0) 


,'.!■  ■  •■■!>-  ±L 


....,—  *  i tt 
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bandwidth  curves.   Also  it  is  an  easy  matter  to  choose  an  M  point  such 
that  it  guarantees  a  given  bandwidth  specification. 

From  Fig.  2-14 ,  the  effect  of  moving  the  M  point  (in  the  stable 
region)  on  system  bandwidth  can  be  analyzed  as  follows s 

Suppose  a  point  M  (a^n.«aQ  )  is  located  in  the  stal  For 

fixed  a1  ,  increasing  a   results  in  increasing  the  bandwidth  as  is 
easily  seen  from  Fig.  2-14.  This  agrees  with  common  eaqper.- 

For  fixed  a.  ,  the  effect  of  changing  aT   ©n  band*  ! 
upon  the  value  of  a~  .  This  can  be  analyzed  by  noting  the  slope  of  con- 
stant Ct),   curves. 

It  is  seen  in  Section  2-4  that  the  locus  of  the  vertex  point  of  the 
family  of  constant  bandwidth  curves 9  the  hyperbolas,  is 

Y  =  0.414  367)((x)t7>) 
where  X  and  Y  axes  are  identical  to  the  B«   and  B^*  axes  respectively 
and  B«  (&>tn)  amJ  B0  ^  ^tn^  are  the  e«iuati'Q)ns  which  define  the  ^  curve. 
Noting  that  the  ■  $  curve  for  the  normalized  system  is  a  straight  line 
emanating  from  the  origin  with  unity  slope  in  the  first  quadrant 
the  locus  of  the  vertex  points  of  the  constant.  ft)fe  curves  is  a  straight 
line  with  slope  of  .414  defined  by 

B0n  =  °414Bln  C2-32) 

on  the  Bq  vs  B,   plane. 

Since  a  vertex  point  of  a  hyperbola  is  the  point  at  which  the  slope 
of  the  hyperbola  changes  its  sign,,  every  constant  <D,   curve  lying  in 
the  stable  region  is  such  that  the  portion  of  the  curve  which  is  above 
(to  the  left  of)  the  line  defined  by  Equation  (2=32)  has  negative  slope 
and  dually. 
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Dividing  the  stable  region  into  two  regions  by  the  line  defined  by 
Equation  (2-32),  it  is  easy  to  see  that  if  a-  is  fixed,  then  by  increas- 
ing a-  ,  the  bandwidth  is  increased  as  far  as  the  M  point  is  in  the 
region  above  the  line  and  dually  for  the  region  below  the  line. 

If  a  closed  loop  system  has  a  zero  in  its  closed  loop  transfer 
function,  then  the  normalized  system  equation  is 
,  bmxi  ■+  2-a  •/? 

fo-^+xf +»„,.*■,  a..  :-"> 

where  b,  is  normalized  b^  ,  the  coefficient  of  S  term  in  the 
numerator  polynomial  of  the  unnormalized  system,  and  the  rela- 
tion between  b^  and  b.   is 

The  arguments  developed  in  Section  2-4  for  a  system  whose  closed  loop 
transfer  function  has  one  zero  applies  to  the  third  order  system  repre- 
sented by  Equation  (2-33). 

By  letting  bin  =  a^n  in  Equation  (2-33),  then  the  resulting  constant 
bandwidth  curves  which  are  defined  by  Equation  (2-24)  in  Section  2-4 
(with  due  account  of  normalized  third  order  F^  )  ares 


(aw+^4Va)^(am+3J^ 


(2-34) 


where  B,   and  B.  are  Mitrovic's  variables  for  the  17  curve 
In      On  ° 

on  the  normalized  third  order  charts. 
The  portions  of  the  ellipses  lying  in  the  stable  region  are  plotted 
in  Fig.  2-16  which  can  also  be  superposed  on  Fig.  2-15  to  provide  infor- 
mation about  bandwidth  for  any  third  order  system  with  one  zero  in  its 
closed  loop  system  transfer  function, 
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Fig  2-16  Constant  bandwidth  curves 

(Normalized  3rd  order  system 


H 


:!:■   ,  ;:r 


Im:] .;;; 


f°r  vS! 

1       .,    ...i..        .         i                  j   ..    ;           ..      i:    .'    ■     ■ 

-^T.-i      :,:.;!;■:■::     1  :  :  I-.  '  ■  :  ■  k  \  A    [ 

1 

As  illustrated  in  Section  2-4,  Fig.  2-16  constitutes  a  bandwidth 
chart  for  a  system  represented  by 

and  a  bandwidth  chart  for  the  upper  bound  of  bandwidth  for  a  system  rep- 
resented by 

kJS\  , ^^ 

K^>        £3-t2Lz£*+2.,S   +  a.Q 

or  equivalently  for  the  normalized  system  in  Equation  (2-33), 

The  effect  of  moving  the  M  point  in  the  stable  region  on  bandwidth. 
can  be  analyzed  in  Fig.  2-16  analogously  as  in  Fig.  2-14. 

Note  that  the  ellipse  defined  by  Equation  (2-34)  has  a  center  point 
moving  along  the  line  Bqd  =  B-^n  in  the  third  quadrant  as  6),   varies , 
and  also  note  that  the  ellipse  has  negative  slope  in  the  first  quadrant. 
Thus  it  is  easy  to  see  that  by  increasing  either  a.   or  a~  ,  the  bandwidth 
is  increased. 

In  Fig.  2-17,  2-18,  2-20  and  2-21,  the  constant  bandwidth  charts, 
Fig.  2-14  and  2-16  are  plotted  on  expanded  scales.   The  third  order  chart 
Fig.  2-15  is  correspondingly  expanded  in  Fig.  2-19  and  2-22  which  go  with 
Fig.  2-17,  2-18  and  2-20,  2-21  respectively. 

The  following  examples  illustrate  determining  bandwidths  of  third 
order  closed  loop  systems  by  using  the  bandwidths  charts  developed  in 
this  section,  and  also  illustrate  a  design  of  system  to  meet  exact  band- 
width  by  using  the  third  order  system  bandwidth  chart. 
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Fig.  2-19  3rd  order  chart (Normalized  B  vs  B  curves) 


Fig.  2-21.  Constant  bandwidth  curves 


Fig.  2-22  3rd  order  chart. (Normalized  B  ys  B.  curves) 
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Example  2.4 

For  the  system  in  Example  2.1 


few- 


the  bandwidth  is  found  to  be  CO,  =  1.76.   Now  one  can  find  the  bandwidth 
on  the  third  order  bandwidth  chart. 

The  given  system  normalizes  to,  by  using  the  transformation  S  -  3ss 
JLuy.   _  -J^l 

and  the  normalized  bandwidth  60,   is  related  with  the  unnormalized  band- 
bn 

width  CO  by  O)  =  3  &)   .   Since  the  system  has  no  zeros,  one  enters  into 
d      b      bn 

the  chart  in  Fig.  2-14  or  the  expanded  chart  in  Fig.  2-17  or  Fig.  2-20. 
The  M  point  of  the  normalized  system  seems  the  chart  in  Fig.  2-17  is 
suitable.   By  locating  the  point  M(.222,  .148)  on  the  chart  in  Fig.  2-17, 
6Jbn  is  found  to  be  .59  by  interpolating  0),       z  .55  and  oO,   =  .6  curves. 
Then  the  bandwidth  of  the  original  system  is 

which  agrees  with  the  result  in  Example  2.1. 
For  the  system  in  Example  2.2 

$S  *~    S3  +  / i>i /./£*+  /22  67/S  +  v-v-8000 

the  bandwidth  is   found   to   be    C0fc   =    110   in  Example  2.3. 

This  system  is  third  order  with  one  zero  in  its  closed  loop  transfer 

function.   If  the  bandwidth  is  determined  from  the  third  order  bandwidth 

chart  for  a,   =  b,    then  the  bandwidth  so  determined  must  be  not  smaller 
In    In 

than  110.   Noting  that  aL  =  122671  and  hl  z   121000  are  not  greatly  dif- 
ferent, the  bandwidth  determined  from  the  chart  for  b^   =  a^  will  be 
very  close  to  the  actual  bandwidth  110. 
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The  system  normalizes  to,  through  the  transformation  S  =  1671.1s, 

C_ ,   v bind   +  ■  £>0t>9l  __ 

R^    '         J3 -f  J2 -r  .oV38A  +  .ooofl 

Entering  into  the  bandwidth  chart  in  Fig.  2-21  with  the  normalized 
M(.0438,  .00096)  point,  the  normalized  bandwidth  is  found  to  be 

^>  fetfj  =  •  o  7 
by  interpolating  between  the  two  curves  for   bn  r  "^  an(^   b  =  «075. 
Then  the  unnormalized  bandwidth  is 

which  is  very  close  to  the  actual  bandwidth  d)     =   110  and  agrees  with 
the  arguments  developed  in  Sections  2-3,  2-4  and  2-6. 
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Example  2.5 

Given  a  system  block  diagram  below 


f-% 

-\ 

K 

c 

\ 

SCS+OU+f) 

JQS*  +  KtS 

adjust  K,Ka  and  Kt  with  the  following  performance  specifications; 

1)  ^   r   .3 

2)  Kv  =   10 

3)  0)b  --  20 

For  the  system  without  feeding  back  derivative  signals  K  must  be 
set  to  90  in  order  to  meet  the  Kv  specification,  and  the  system  is  at 
stability  limit  which  can  be  seen  from  the  third  order  chart. 

The  closed  loop  transfer  function  with  the  derivative  signals  feeded 
back  is 


K 


and  normalizes  to,  by  using  the  transformation  S  =  a^s  =  (10-+-KKa)ss 


^<\_  . 


where ;   aln  -  (  9  +  AT  Kt )  /0  0  +  K  Kd  f 

a0n  =  F/('o  +KK3.i 
and     a3  =  10 f-  KKa  =  normalizing  factor. 
The  velocity  constant  and  bandwidth  of  the  unnormalized  and  normal- 
ized systems  are  related  by 
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Kv  =  (l0+KKa^v7, 

COb=({0-r  K  Ka)  tijbn 

Kv/6)b=  K^  /  <^>u»  =  '/«- 

where  the  subscript  "n"  denotes  the  normalized  system  quantity. 
Thus  the  ratio  of  the  velocity  constant  Kvn   and  the  bandwidth  6)bn  of 
the  normalized  system  must  be  maintained  at  1/2. 

In  Fig.  2-22a  the  'o  and  1,3  curves  are  duplicated  from  the  third 
order  chart  (Fig.  2-15)  and  then  some  constant  bandwidth  curves  are 
superposed  from  the  third  order  system  bandwidth  chart  (Fig.  2-14). 

In  order  for  the  system  to  have  "$=    .3  the  M  point  must  be  located 
below  the  I  .3  curve  in  the  first  quadrant  on  the  Bgn  vs  B^n  plane.   For 
the  normalized  system  Kvn  =  a0ny'aln  indicates  that  for  an  M  point  located 
anywhere  on  the  normalized  Bq  vs  B^  plane  if  a  line  is  drawn  from  the 
origin  to  the  M  point,  then  the  slope  of  the  line  is  equal  to  the  normal- 
ized velocity  constant  Kvn.   Therefore  for  this  problem  the  M  point  must 
be  at  the  intersection  point  of  a  constant  CO^  =  ^b'n  curve  and  a  line 
drawn  from  the  origin  such  that  the  slope  of  the  line  is  (l/2)£0bn  and 
moreover  that  M  point  must  be  such  that  it  insures   ^=  .3.   For  any 
chosen  constant  &)bn  =  ^bn  curve  tne  corresponding  straight  line  with 
slope  1/2  CdL*   can  be  drawn  and  then  the  line  almost  always  intersects 
with  the  constant  (j)un  =    &)b'  curve.   ./hen  intersection  occurs,  location 
of  the  M  point  at  the  intersection  guarantees  that  both  Kv, bandwidth 
specifications  have  been  satisfied,  but  may  not  guarantee  the  "£  specifi- 
cation. In  order  that  all  three  specifications  be  satisfied,  the  inter- 
section must  lie  in  the  area  enclosed  by  the   .3  curve  and  the  B^n  axis 
in  the  first  quadrant, 

71 


iBcn 


To 


'^  ^,71  =1.0 


.  ^r~  Constant  Q^n  curve 


■&** -9 


"Kvn=.35  line 
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C0brj=.6 


n. 


^(T^sSCM)' 


•8         .9         1.0   B_ 


Fig.  2-22a  Solution  of  Example  2.5  using 

the  3rd  order  system  band  width  chart. 
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If  one  chooses  the  constant  bandwidth  curve  for  uX   =0.7  and 
draws  the  line  with  a  slope  of  0.35  from  the  origin  as  shown  in  Fig. 
2-22a,  then  the  intersection  point  is  at  (.575,  .204)  and  this  point  in- 
sures  ^=  .3  which  is  evident  from  the   1.3  curve  configuration,  thus 
one  chooses  the  M  point  at  M(a^n,  a«  )  =  M(.575,  .204)  to  give 


9  i-  KKt 
K 


a0n- 


.5-75" 


—  =  .244 


Since  6J>bn  :  .7  is  chosen  the  normalizing  factor  a^  =  10-f-KKa  is  de- 
termined from  the  relation 

tOj,  =(lO  +  KKa)cdh71^\C) 
from  which  10+*KKa  =  28.6  is  obtained.   Substituting  this  value  of  the 
normalizing  factor  into  the  denominators  of  the  two  equations  for  the 
coordinate  of  the  selected  M  point  above,   K,Kt  and  then  Ka  are  deter- 
mined and  the  solutions  are 

K  r  4675 

Kfc  =  0.0987 

Ka  =  0.00388 


2-6.  Frequency  response  of  the  fourth  order  systems. 

The  lo  curve  of  a  fourth  order  system  can  be  normalized  so  that  a 
single  curve  can  represent  the  l^  curve  of  any  fourth  order  system. 

Suppose  a  fourth  order  characteristic  equation  is; 

F(S)  r  S4+  a.,?T+  a^S^-  4,£  ■+  &0  =0 

then  the    \0  curve  is  defined  by 

(3|Sb  a,aT 

a       ¥  (2-34a) 

Equation  (2-34a)  represents  a  concave  down  parabola  whose  axis  of  sym- 
metry is  Bi  =  aoao/2  with  vertex  at  (a2a3/2,  3233/4),  Bq  =  0  occurs  at 
CO  a  0  and  CO-   a|  . 

h 
Assume  a  normalized  frequency  ft)   such  that   a)s  a  4)    then 

Equation  (2- 34a)  becomes 

3,  =  3.*  a3  a>£,  (2-35) 

2 

Let  B       =   B,/a2a~     and     B       =   *Wa2        t*ien  E<luation   (2-35)   can  be  writ- 
ten as 

(3,^  =  0)^  (2-36) 

Equation  (2-36)  is  exactly  the  same  equation  for  the  Iq  curve  as  is  asso- 
ciated with  a  characteristic  equation  whose  coefficients  of  the  S  ,  S 

2 
and  S   terms  are  unity. 

If  a  curve  defined  by  Equation  (2-36)  is  constructed,  then  with 

proper  scaling  factors  of  coordinates  and  frequencies,  it  will  represent 

exactly  the  same  curve  represented  by  Equation  (2-34a)  and  thus   f^ 

curves  of  fourth  order  systems  are  normalized. 
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The  normalizing  property  of  the  \0   curve  for  fourth  order  systems 
can  be  used  to  obtain  the  frequency  response  of  any  fourth  order  system. 
Suppose  a  fourth  order  system  closed  loop  transfer  function 

which  has  a  second  order  polynomial  in  the  numerator.  For  the  case  of  a 
fourth  order  system  with  third  order  numerator,  add  b3S  to  the  numerator 
of  Equation  (2-37).   By  using  a  transformation  or  frequency  scaling 

Equation  (2-37)  can  be  written  as 


af^^-f  a5  a/^r-r-  &?*$*+  a,  3/^  +  a0 


The  systems  represented  by  Equation  (2-37)  and  (2-38)  are  equivalent 

since  the  transformation  S  =  a|s  is  one  to  one.  Note  in  Equation  (2-38) 

that  the  coefficients  of  the  S^  and  S^  terms  in  the  denominator  are  made 

unity,  and  also,  the  coefficients  in  the  numerator  and  denominator  of 

Equation  (2-38),  by  comparing  those  of  same  power  terms,  are  scaled  by 

-3/2 
the  same  amount.  For  instance,  b^  and  a^  both  are  multiplied  by  a2         , 

and  likewise  for  a 2  and  b2   >  etc. 

Let  a0n  =  a./^ 

aln  —  <3,/<2z2.5 

a2n  =a2/^    =1 

a3n   =  a5/d3  =1  (2-39) 

bln   =  b,/4,a3 

b2n    =  b*  /  &*. 

a2t   =a/7a3 


then  Equation  (2-38)  is 


T«- 


b>27)  st*  a2t  kn4+2* 


4*+  *£?  *3  + /s~+  2^pat7iy<i+  a.o7l 


(2-38a) 


,.l* 


By  substituting  s  -  j  ^  ,  and  noting  that  Bln  z   a)  tn»  BQn  s  a)  £n  -  (0t„ 
which  are  the  equations  defining  the  normalized  ^  curve  for  fourth  order 
systems,  the  magnitude  of  the  frequency  response  is  given  by 


^i^fn) 


1  (doy-ky^tn  )  +j  4zt  bin  C0tn  \ 


(2-40) 


thus  the  frequency  response  of  fourth  order  systems  is  represented  in 
terms  of  normalized  quantities. 

The  methods  described  in  Sections  2-1,  2-2,  2-3  and  2-4  can  be  ap- 
plied to  get  the  frequency  response  of  the  fourth  order  system.  The 
frequency  response  of  the  original  system  represented  by  Equation  (2-37) 
can  be  obtained  from 


T\J 


((CO) 


(a0-kiox)+Jb,£0 


(2-41) 


which  results  from  applying  Equation  (2-6)  to  the  system  (2-37).   The 
advantage  of  using  normalized  Equation  (2-40)  over  the  original  Equation 
(2-41)  is  that,  since  the  normalized  |0  curve  for  any  fourth  order  sys- 
tem is  a  fixed  curve,  it  is  not  required  to  plot  the  \Q    curve  for  each 
system. 

Knowing  how  to  interpret  Equation  (2-41)  on  the  Bq  vs  B,  plane  which 
is  associated  with  the  original  system,  Equation  (2-40)  can  be  inter- 
preted on  the  normalized  B^  vs  B^n  plane  without  difficulty. 
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The  quantities  or  coordinates  in  the  normalized  and  unnormalized 
systems  or  Mitrovic's  plane  are  related  by  Equation  (2-39)  plus  the  fre- 
quency relation 

CO  =  3~2  cOtn 
The  M  point  is  located  on  the  normalized  Bqu  vs  B^n  plane  by  the  normal- 
ized coordinate  (a.   ,  a_  ). 

The  imaginary  term  in  the  denominator  of  Equation  (2-40)  can  be  rep= 

resented  in  the  same  way  as  for  Equation  (2-41),  namely  by  CO  (a^  -  B.) 

-%  -k 

plot;  however,  it  must  be  multiplied  by  a«t  to  get  the  a,^  ^Vn^lrT8!  ^ 

plot,  and  likewise  for  the  imaginary  term  in  the  numerator  of  Equation 
(2-41)  if  a  system  of  interest  has  b,  other  than  zero. 

In  Fig.  2-23,  the  normalized  VQ    curve  for  fourth  order  systems  is 
plotted  on  the  Bqr  vs  Bi   plane  with  the  same  scales  for  both  axes.  Any 
point  M„(a,  ,aQ  )  located  on  the  Bqq  vs  B^n  plane  is  equivalent  to  the 
point  M(a,,a_)  on  the  Bq  vs  Bi  plane  where  a,  ,  a^  and  a  ,  a  are  re- 
lated by  Equation  (2-39).   Any  frequency  COtn   on  the  normalized  V0   curve 

i 
is  equal  to  CO   divided  by  a~  where  CO   is  the  frequency  on  the 

unnormalized  V0   curve  at  the  point  which  is  equivalent  to  the  6i)   point 

on  the  normalized  \0   curve. 

The  following  table  is  a  comparison  of  the  V0   curves  for  normalized 
and  unnormalized  systems; 

Normalized  \0   curve   Unnormalized  l0  curve 

Axis  of  symmetry  Bin  r  0.5  B^  =  .5  8l2&3 

Coordinates  of  vertex         (.5,  .25)  (.5  a2a3  ,  .25  a|  ) 

Frequency  at  vertex  C0tn  =  .707         CO  =      .707/^2 

Freq.  at  B^  intersect  COtn  =   0,  C0tn  =  1     COs   0S  CO-  J~*2, 

C°indCe"::rt  "  Bl  ».  •>•  d.  °>     (0,  0).  (O,  a2a3) 
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The  equation  defining  constant  bandwidth  curves  of  fourth  order  sys- 
tems on  the  normalized  BQ  vs  B,  plane  contains  a  parameter  Sot5*  therefore 
it  is  not  a  practical  way  to  try  to  provide  a  universal  chart  of  constant 
bandwidth  curves  of  fourth  order  systems ,  because  for  each  &2t   there  re- 
sults a  family  of  constant  bandwidth  curves,  Howevery  it  would  be  worth 
while  to  investigate  the  effect  of  varying  a**.  on   bandwidth ,  since  if 
the  effect  of  a£t  on  bandwidth  is  known,  then  if  a  constant  bandwidth 
curve  is  drawn  or  bandwidth  is  determined  for  some  a2t  then  it  may  be 
possible  to  get  some  information  about  bandwidth  for  another  &2f 

Assume  two  fourth  order  systems  such  that  when  they  are  normalized 
they  result  in  identical  system  equations  except  they  have  different 
a2t's,  namely  a2*t  =  k^  for  one  of  the  two  systems  and  a^t  =  k2  for  the 
other  where  k^  and  k2  are  constants,  k^  k2  . 

The  two  systems  have  a  common  M  point  on  the  normalized  Bq  vs  B^ 
plane.   Suppose  this  M  point  is  located  as  shown  in  Fig.  2-23.   Since 
the  normalized  \0   curve  and  M  point  are  common  to  both  systems s  the 
&)tn(a.  -  B|  )  plot  will  also  be  common,  but  each  system  will  have  a  dif- 
ferent  -^3-  (aln-  Bln)  plot  as  shown  in  Fig.  2-23.   (kj_  z   1  and  k2  s  2 

are  assumed  for  a  simple  illustration.) 

Since  the  bandwidth  of  a  system  is  obtained  at  the  point  where  some 
magnitude  circle  intersects  the    -=^  (aln  -  B-  )  plot9  and  noting  that 
the  smaller  a2t  causes  the     "Tip  (aln  "  Bln)  Plot  to  be  displaced  farther 
from  the  B,   =  a1   line  than  the  larger  a2t  doess  any  circle  centered  at 
the  M  point  intersects  the    7== (aln  "  Bln^  plot  at  lower  frequency  for 
the  smaller  a.,  .   For  example  in  Fig.  2-23  s  if  the  system  has  no  zeros 
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in  its  closed  loop  transfer  function,  then  the  magnitude  circle  CAb  to 
determine  bandwidth  has  radius  equal  to  fl   a~  s  .0707  and  it  intersects 
with  the  two    — =(aln-  Bln)  plots  at  Pj  and  P2  to  give  the  bandwidths 

03bn  =  .342     for     (Tlt     ;   1 
<^bn  =  «35      for     /~^Cc  -  2   as  showilu 
It  is  obvious  the  above  reasoning  holds  for  systems  whose  closed 
loop  transfer  function  has  any  number  of  zeros s  thus  a  conclusion  Is  ob- 
tained as  follows: 

The  effect  of  a„  on  bandwidth  of  a  normalized  fourth  order  system 
is  to  increase  the  normalized  bandwidth  with  increasing  Ant' 
It  must  be  remembered  that  the  normalized  bandwidth  &)bn  is  related  with 
unnormalized  bandwidth  60.  by 

^-/^^ 
As  a  simple  application  of  the  above  conclusion9  refer  t©  Fig.  2-13,   In 
Fig.  2-13,  the  hyperbolas,  the  constant  bandwidth  curves  are  for  a  system 

^  a" 

This  system  can  be  regarded  as  some  system  which  normalizes  so  that  the 
coefficient  of  the  S   term  is  unity  with  a„  -   1.   In  other  words g   if  for 
a  normalized  system 

if  constant  CO.      curves  are  plotted  with  a„  r  1,  then  there  will  result 
bn  2t 

exactly  the  same  configuration  of  curves  as  in  Fig.  2-13.   For  the 
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M(a,  ,  a«  )  point  as  indicated  on  the  curve  C  c  ,  O       is  known  to  be  .5, 
If  there  is  another  normalized  system  with  the  same  M  point  but  with  a«t 
slightly  different  from  1,  then  CO.   of  this  system  will  not  be  far  off 
from  .5  in  the  well-known  direction  of  deviation  from  .5.   Same  reasoning 
can  be  made  if  the  M  point  is  located  anywhere  on  the  BQ  vs  B1  plane. 

The  following  numerical  examples  are  to  illustrate  using  the  normal- 
ized \0   curve  in  obtaining  frequency  response  of  fourth  order  systems. 
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Example  2.6 

Given  a  closed  loop  system  transfer  function 

T(S)=     S*+22S3+,rzSSL+ IMS +  1*00  (2"42) 

obtain  the  frequency  response  (magnitude). 

The  problem  can  be  worked  out  without  normalizing  the  systems ;  how- 
ever, it  will  be  shown  on  the  normalized  BQn  vs  B^n  plane. 
For  the  given  system,  the  coefficients  are 
aQ  =  1600 
ax   r  680 
a,  =  172 


22 


*3 
By  using  relations  between  normalized  and  unnormalized  quantities  in 

Equation  (2-39),  the  coefficients  are  normalized  to 

£0n=  a,/  aj  =.oS4-2 

and  fet=fax/a3~.&S 

The  frequencies  of  the  normalized  and  unnormalized  system  are  related 

From  Equation  (2-40),  the  frequency  response  (magnitude)  of  the  given 
system  is,  in  terms  of  normalized  quantities, 


¥^\ 


.o£-*Z 


^z-E^+j-^.tf-B,^ 


In  Fig.  2-24  the  point  M(aln ,  aQ  )  =  M  (.18,  .0542)  is  located  and 
the  lines  Bln  =  aln  =  .18,  B0n  :  a0n  =  .0542,  B0n  s  (1  +  f!  )aQn  z    .131 
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are  drawn  on  the  normalized  Bq  vs  B*  plane. 

Since  the  line  Brt  z  (1  -+-/~2)a.  =  .131  crosses  the  '  o   curve  at  the 
un  un 

point  where  the  frequency  is  slightly  less  than  0)       -   .49  and  since  the 
system  has  no  zeros,  the  bandwidth  <D,   must  be  less  than  .4.   From  the 
|fl  curve  and  M  point  configuration  on  the  normalized  Bq  vs  B^  plane  as 
shown  in  Fig.  2-24,  it  is  seen  that  the  first  one-half  part  of  the 
curve  is  sufficient  to  see  the  frequency  response  of  the  system. 

In  Fig.  2-25,  the  \0  curve,  M  point  and  the  B.  -  a.  ,  BQn  z  &Qn 
and  B0  =  (1  +  J~2  )a~  lines  are  drawn  to  expanded  coordinate  scales. 
Frequency  response  is  represented  in  the  same  way  as  in  Fig.  2-8 ,  the 
only  different  thing  from  Fig.  2-8  is  the    -j=^-  {&\n   -  B|n)  plot. 

As  illustrated  previously,  the    -=  (a-   -  B,  )  plot  is  con- 
vert 

structed  such  that  the  horizontal  distance  from  the  line  B,   =  a,   =  *  18 

to  any  point  on  the    -7=-  (a,   -  B,  )  plot  is  equal  to  the  horizontal 
J82+      In    in 

distance  from  the  line  Bin  =  a^n  =  .18  to  the  1^  curve  multiplied  by 
J^L    m  J^n   m      The   %  , 

!0  if  the  ordinates  are  the  same  as  in  the  &>(a,  -  B.)  plot. 

Some  magnitude  circles  are  drawn  to  show  the  magnitude  A  z  LOtU^ 

at  some  frequencies.  The  magnitude  circle  tangent  to  the  ~z=  (ain~Biin) 
plot  which  is  for  resonant  point  touches  to  the   r=r  (a-,  =  B,  )  plot  at 

C0tn   =  .24  with  radius  equal  to  .05  to  give 

the  resonance  frequency   M   =  .24  or   0)r  *  (13.1)&)rn  =  3.14 
the  resonance  peak       fl^  =  .0542/  .05  =  1.08 

The  magnitude  circle  with  radius  equal  to  /2  &0  :  .765  crosses  the 

Ok 
-—^-{^    -Bln)  plot  at  0)tn  =  .388  to  give  the  bandwidth 
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60bn  =  .388   or   60b  r  (13.1)  (.388)  =  5.09 
Substituting  S  =j&)   in  Equation  (2-42),  the  magnitude  for  Cd  ~    a)     = 
3.14  and   0)=  (^  =   5.09  calculate  to 

— -(j3.i¥)\z   1.1  for  resonance  peak  MpW 

and  Ad'^o^l:  •  71  for  bandwidth 

Magnitude  circles  are  drawn  in  Fig.  2-25  for  several  frequencies  and  mag- 
nitudes are  indicates,  thus  the  frequency  response  is  obtained. 
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Example  2.7 

Given  an  open  loop  transfer  function  for  a  unity  feedback  control 
system  as 


3J.6(s+\)(S+5) 


determine  the  bandwidth  of  the  closed  loop  system 
The  closed  loop  system  transfer  function  is 


£fe)- 


(2-43) 


which  has  a  second  order  numerator. 

The  quantities  to  be  normalized  to  work  on  the  normalized  BQ  vs  B, 

plane  are,  using  the  relations  in  Equation  (2-39), 

a-  =  158/55.6  r  .0512 
On 

aln  =  196.66/(12.333)  (55.6)  =  .286 
bln  =  190/(12.333)  (55.6)  =  .277 
b2n  =  31.6/55.6  =  .568 
/ajt  =755.6/12.333^  =  .603 
{&i  z   7.44 
6)  =  7.44  0)tn 

The  frequency  response  of  the  system  in  terms  of  normalized  quantities 
is  defined  by,  using  Equation  (2-40), 


f<j*W 


(-05l2-.5&M&)+j^£0*'7'7) 


(.0512-0^)^  $*(.2*i-BrJ 


For  the  bandwidth 


|(«2-^Vj|i^^W|-^|(.«ia-W^)+j||(.aT^ 


(2-44) 
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The  left  side  of  Equation  (2-44)  can  be  represented  by  the  usual 
-g^-(.286  -  Bln)  plot.   On  the  right  side  of  Equation  (2-44),  both  real 
and  imaginary  terms  are  frequency  dependent.   As  illustrated  in  Section 
2-3,  the  real  term  with  the  factor  fl   is  represented  on  the  B,n  =  a^n  = 
.286  line  and  the  imaginary  term  with  the  factor  fl   is  represented  on 
the  BQn  =  a^  =  .0512  line.  Those  are  shown  in  Fig.  2-26. 

In  Fig.  2-26,  the  normalized  1^  curve  with  the  point  M(.286,  .0512) 
is  shown.  On  Bln  =  .286  line,  fl  ^20.^^.  ~    ,805  ^tn  is  calibrat«d  as 
follows.   The  point  corresponding  to  60t_  =  0  which  is  located  below  B, 
axis  is  at  the  distance  of  fl  aQ  =  .0723  from  the  M  point.  From  this 
point  up,  the  Bln  r  .286  line  is  labeled  with  6)tn  such  that  the  distance 

from  the  CO   =0  point  to  a  UJL  =  CO.,    point  is  equal  to 
tn  tn     1 1 

2         2 
fl   bln  0dtl  z   .805  (K>tl     so  that  at  00 

M  point  to  the  0)  fc  z    0)t]_  point  on  the  B^n  =  .286  line  is  equal  to 

|yz(X.-^tt)H-o7*3--*°^l 

For  example  at  C0tn  =   .9  ,    fl   (aQn  ~  bln  ^tn^    is  e<*ual  to  the  dis~ 
tance  from  the  M  point  to  the  6)tn  r  .9  point  on  the  Bln  r  .286  line, 
which  is  equal  to  .73. 

The  imaginary  term  in  the  right  side  of  Equation  (2-44)  is  repre- 
sented on  the  BQn  r  aQn  r  .0512  line  as  follows.  On  the  B0n  =  aQn  -z  .0512 
line,  each  point  to  the  right  of  the  M  point  bears  frequency  a)fc  such 
that  a  point  corresponding  to  Cdtn  z     0)ti     on  the  BQn  =  .0512  line  is 
located  at  a  distance  of  fl   bin^tl/  /  a2t  =  *^^ti   .   For  example, 

at  6)_  :  ,9  on  the  BA  =  .0512  line,  the  distance  from  the  M  point  is 
tn  Un 
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(.65) (.9)  r  .585,  thus  the  line  BQn  =  .0512  has  linearly  displaced  a>tn 
whereas  the  line  B,   r  .286  has  quadratic  o)   . 

With  those  two  lines,  the  quantity  on  the  right  side  of  Equation 
(2-44)  can  be  read  off  at  any  frequency.   For  example s  at  CO   s  .9  the 

quantity  fl  |(.512  -  .568  0)^)  +j-^§(.277)   is  the  length  of  the  diagonal 

of  the  rectangle  formed  by  the  point  M,  Q)        z    .9  points  on  B1n  r  .286 

tn  in 

line  and  Bq  r  .0512  line,  which  is  the  radius  of  the  magnituce  circle 

CA2. 

Then  the  bandwidth  is  determined  in  the  same  way  as  in  Example  2.3« 
In  this  example,  as  seen  in  Fig.  2-26, 

<°bn  •  -9 
Noting  that  frequencies  are  scaled  by  Cd=   7.44  fi)   ,  bandwidth  of  the 

original  system  is 

"b  = 6-7 

By  substituting  S  r  j  6.7  in  Equation  (2-43),  the  magnitude  calculates 


|i-Cj6.7)|=.7»7 


1TJ' 
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2-7.  Evaluation  of  phase  angle. 

In  the  previous  sections  the  magnitudes  of  the  frequency  response 
of  a  system  are  evaluated  on  the  Bq  vs  B^  plane.   Since  the  numerator 
and  the  denominator  of  the  equation  defining  frequency  response  can  be 
represented  as  vectors  on  the  Bq  vs  B^  plane,  it  is  obvious  that  the 
phase  angles  in  the  frequency  response  of  a  system  can  also  be  evaluated 
on  the  Bq  vs  Bi  plane. 

In  this  section  the  method  for  evaluating  the  phase  angles  in  the 
frequency  response  of  a  system  is  illustrated  with  some  of  the  numerical 
examples  which  are  considered  in  the  previous  sections  of  the  magnitude 
evaluation. 

Again  writing  the  general  equation  for  the  frequency  response  of  a 
system  as 


the  phase  of  the  frequency  response  is 

2^£^=Z.tUj^  -/.((ao-Bo^-tW^-^S)         (2-45) 

K  0 

Since  the  vectors  N(  j  0)  )  and  (a  -  B  )  +.'t0  (a^  -  B^)  are  measured 
from  the  M  point  and  the  real  and  imaginary  parts  of  the  vectors  are 
represented  on  the  BQ  vs  B,  plane  in  such  orientations  that  they  are 
parallel  to  the  Bq  and  B^  axes  respectively,  a  proper  reference  in  meas- 
uring the  angles  must  be  set  up. 

Consider  the  angle  ^K(aQ  -  BQ)  +;M(a.     -   bq))  with  the  aid  of 
Fig.  2-27.  The  4)  (a-^  -  B,)  plot  is  constructed  as  illustrated  in  the 
previous  sections  with  the  [^  curve  and  the  M  point  configuration  as 
shown.  Select  a  point  such  as  P  in  Fig.  2-27  at  frequency  of  u3j_  and 
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foUj-  Bt)  plot 


270°    Bo  =  a„  line 


Fig.  2-27  Evaluation  of  the  phase .angles 
in  the  frequency  response. 
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draw  the  vector  PM.  The  components  of  this  vector  define  the  real  and 
imaginary  terms  (aQ  -  Bq)  and  6\(aj_  -  B^)  in  magnitude  and  sign.   The 
horizontal  component  6).(a^  -  B^)  is  the  imaginary  part  and  positive  when 
directed  to  the  right.  The  vertical  component  (aQ  -  BQ)  is  the  real  part 
and  positive  when  directed  upward.   Considering  the  signs  of  the  compon- 
ent vectors,  the  angle  of  the  vector  PM  =  (a.Q   -  Bq)  ■+■;  CO^a^  -  B^)  is  the 
angle  B-,   which  is  measured  from  the  vector  (aQ  -  Bq)  to  PM,  and  this  is 
equal  to  the  angle  measured  from  the  line  B^  r  a^  below  the  M  point  to 
MP  and  positive  when  measured  to  the  clockwise  direction.   The  angle  of 
the  vector  N( j  "> )  is  naturally  62  which  is  measured  from  the  line  B^  =  a^ 
above  the  M  point  to  the  vector  N(jA>  )  with  the  positive  direction  as  in- 
dicated, thus  the  phase  angle  of  the  system  in  the  frequency  response  at 
frequency  CO.  is 


^  ■»■<]■«£-  e*-&4 


It  must  be  noted  that  the  reference  line  of  the  angle  of  N(  led )  de- 
pends upon  how  the  vector  N( ;co  )  is  represented  and  it  is  seen  that  the 
vector  N(j60  )  can  be  represented  either  to  the  right  or  left  of  the  line 
Bi  =  a,  and  either  above  or  below  the  line  Bq  -  aQ  ;  however,  for  the 
angle  of  the  vector  (a0  -  B0)  HCO  (a,  -  B,),  the  reference  line  is  always 
the  line  B-,  =  a,  below  the  M  point  and  the  positive  angle  is  measured 
from  the  reference  line  to  the  vector  MP  in  the  clockwise  direction,  thus 
the  0°,  90°,  180°  and  270°  positions  are  as  indicated. 

Evaluation  of  phase  angles  in  the  frequency  response  of  an  open  loop 
system  is  also  possible  as  the  case  of  obtaining  magnitude  in  the  fre- 
quency response.  The  following  examples  illustrate  evaluation  of  phase 
angles  in  the  frequency  response  for  closed  loop  systems  as  well  as  for 
open  loop  systems. 
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Example   2.8 

For   the   system  equation  in  Example   2.1  which  is   rewritten  below 


G(s)=- 


5(StiXS  +  2) 

obtain  the  phase  characteristics  in  the  frequency  responses  for  both 
closed  loop  and  open  loop  systems. 
Closed  loop  system 

The  closed  loop  system  has  the  transfer  function 


&®s      S3  +  5S*  +2S  + 


S3   +  3S*   +  2.£  +  4- 
in  which  the  numerator  polynomial  is  of  order  zero  and  the  angles 
Z  N(j&>)  =  0  for  all  frequencies,  therefore  only  the  angles 
^  ((aQ  -  BQ)  4-J60(a  -  B  ))   are  frequency  dependent. 

In  Fig.  2-28,  the  ^c,   curve  and  the  CO  (a,  -  B,)  plot  are  duplicated 
from  Fig.  2-8  which  is  for  obtaining  the  magnitude  of  the  closed  loop 
frequency  response  of  the  same  system. 

The  reference  angles  0  ,  90  ,  180  and  270  are  indicated.  The 
radial  vectors  directed  from  the  M  point  to  the  <D(a,  -  Bj)  plot  give 
the  angles  of  the  denominator  of  the  closed  loop  frequency  response  equa- 
tion which  is  the  vector  (4  -  BQ)  -4-jO>(2-B..)  .  Noting  that  the  angles 

o 
measured  from  the  0  reference  line  to  the  vectors  are  of  the  denominator 

and  also  noting  that  the  angles  of  the  numerator  are  zero  degrees  for  all 

frequencies,  the  phase  of  the  closed  loop  system  at  each  frequency  is  the 

negative  of  the  angle  of  the  vector  so  measured,  thus  at  a):  0  the  phase 

is  zero  degrees  and  as  0)    increases  the  phase  increases  in  the  negative 

direction.  For  example,  at  CO r  1.0  the  angle  of  the  vector  measures 

45  and  the  phase  of  the  system  is  -45  .  The  phase  assumes  -180  at 
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'    (a,-  B,  )  plot 
/ 

P« -l.  a)^=1.76 


Gain  cross-over  at  =1.05 
phase(open  loop)=-l68° 


Bo=4  line 


270° 
The  mag.  "circle  for  the  gain 
cross-over. (Centered  at  Mo,  r=4 ) 

Fig.   2-28  Evaluation  of  phase  angles. 
(Example  2.8) 
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CO  z  1.4  and  at  this  point  the  '0  curve ,  CO  (a^  -  B^)  plot  and  the  line 
Bl  =  al  =  ^  intersect.  If  a  protractor  is  overlapped  on  the  figure  such 
that  the  center  point  of  the  protractor  coincides  with  the  M  point,  then 
one  can  read  off  the  phase  at  each  CO   conveniently. 
Open  loop  system 

Rewriting  the  open  loop  transfer  function  as 

one  can  construct  the  \Q   curve  associated  with  the  denominator  polynomial 
of  the  open  loop  transfer  function.   By  comparing  the  equation  with  the 
closed  loop  transfer  function,  it  is  obvious  that  the  \0   curve  for  the 
open  loop  system  is  identical  with  that  of  the  closed  loop  system.  The 

Mn  point  defined  analogously  as  the  M  point  is,  by  noting  that  the  co- 

o 
efficient  of  the  S   term  in  the  denominator  polynomial  of  the  open  loop 

system  is  zero,  at  B^  =  2  on  the  B,  axis  and 

is  located  in  Fig.  2-28  as  shown.   Since  the  Mn  and  M  points  are  on  the 
same  line  B.  =  2,  and  the  17  curves  for  the  two  transfer  functions  (open 
and  closed  loop)  are  identical,  the  d)(a^  -  B^)  plot  is  not  altered. 

By  applying  the  same  principles  to  the  open  loop  system  the  phase 
can  be  evaluated  with  the  vectors  directed  from  the  Mq  point  to  the 
CD  (a1  -  Bj_)  plot.   Naturally  the  90°---270  line  must  be  shifted  down  to 
the  B^  axis.  Noting  that  as  CO  tends  to  zero  the  ^(a^  -  B^)  plot  ap- 
proaches the  Mq  point  from  the  direction  of  the  origin  the  phase  at 
CO  =  0  is  -90  .   As  CO  increases  the  phase  increases  in  the  negative 
direction.   Some  vectors  representing  the  phase  of  the  open  loop  system 
are  drawn  with  dotted  lines,  for  example  at  tf)  =  .5  the  angle  of  the 
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vector  measures  130  to  give  a  -130  phase. 

At  CO  =  1.42  which  is  the  point  where  the  phase  of  the  closed  loop 
system  is  -180°,  the  open  loop  system  also  assumes  -180°  phase  and  it  is 
interesting  to  note  that  this  occurs  for  any  system  which  has  a  frequency 
independent  numerator  polynomial. 

As  CO  tends  to  infinity,  from  the  behavior  of  the  i0  curve  and  the 
CO  (a.  -  B,)  plot  it  is  seen  that  the  phase  approaches  -270°  which  is  the 
expected  phenomenon. 

In  the  Bode  diagram  manipulation,  the  gain  crossover  frequency  and 
the  phase  at  that  point  are  used  for  the  stability  check  and  as  a  figure 
of  merit  in  a  design  problem.  The  gain  crossover  point  which  is  the  zero 
db  or  unity  magnitude  point  can  be  found  simply  by  drawing  a  magnitude 
circle  with  radius  equal  to  aQ  =  4  centered  at  the  Mq  point,  which  is  to 

account  the  magnitude  of  the  open  loop  system,  thus  the  gain  crossover 

o 
frequency  is  found  to  be  1.05  and  the  phase  at  this  point  is  -168  to  give 

+12  phase  margin.  This  phase  margin  indicates  a  stable  but  too  oscilla- 
tory system  which  checks  with  the  Iq  curve  and  the  M  point  configuration. 
The  gain  margin  can  be  found  by  drawing  another  magnitude  circle  centered 
at  the  Mfl  point  passing  through  the  -180  phase  point. 
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Example  2.9 

For  the  system  in  Example  2.2 

\2.\000(S  +  4-7) 


6(5)  = 


the  transformed  closed  loop  system  equation  is,  by  using  the  transforma- 
tion S  =  100s 


fM  = 


In  order  to  obtain  the  phase  characteristic  of  the  closed  loop  system 
the  Vb  curve  and  the  6)t(a  -  B,  )  plot  of  Fig.  2-9  are  reproduced  in 
Fig.  2-29. 

The  angles  of  the  denominator  vectors  &      are  measured  in  a  similar 
way  as  in  the  previous  example.   Noting  that  the  system  has  the  first 
order  numerator  and  the  vector  N(  j'4>t)   at  each  0d     is  represented  by  a 
vector  directed  from  the  M  point  to  the  line  Bqc  s   2a   (the  two  vectors 
for  CO  z    .2  and  ft),  z    .3  are  shown  in  the  figure  with  dotted  lines)  the 
angle  Z_N(  |0)  )  at  each  CO     is  measured  from  the  line  B.   =  a.   above  the 
M  point  to  the  vector  N(j'A)t)  at  each  C0t.  For  example  at 
COt  r  .2^Zn(j  .2)  r  29  .   Denoting  the  angles  of  the  denominator  and 
the  numerator  vectors  at  the  same  frequency  by  0,  and  5L  respectively, 
then  the  phase  of  the  closed  loop  system  L  ■—  (  j'&)t)  =  ^2  "  ^1  »  ^or 
example  at  0)fc  r  .2,  Z.A(j.2)  =  29°  -  31.5°  «  -2.5°. 

For  the  open  loop  system  the  transformed  system  equation  is 

The  ]~l   curve  for  the  open  loop  system  is  not  changed  from  that  of 
the  closed  loop  system  but  the  Mq  point  is  moved  to  the  point  (.167,  0) 

L.IIC   S     <tllU   " 

98 


1    B°: 

I13 

B,t=alt  =12.267 

u« 

r 

\        32 

1 

\ 
u>t(alt-B,t)   plot 

• 

\ 

1     28 

\ 

\ 

J/.3 

» 

\ 

V 

x/.2 

I 
\ 
\ 

V 

» 

\ 

1 

20 

\ 
\ 
\ 

\ 

\ 

0 

1 

16 

-      J\ 

\ 
\ 

V 

12 

u 

\ 
\ 
\ 
\ 

\ 

•       • 

/.5       A*      M       1.2       /•/ 

1 

\ 
9        -f      \;7       .« 

.5         .*        .3          Z        -l^r 

Bbt=2a.* 

8 

■  I'7 

I"6 

> 

T-T(l2. 267, 4.48) 

4- 

p5 

■i\                        y} 

f"3 

5° 

•2 

*^.. 

16   Bit. 


Fig.  2-29  Evaluation  of  phase  angles. 
(Example  2.9). 
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denominator  of  the  open  loop  transfer  function.   This  implies  that  a  new 
(i)  (a,   -  B,fc)  plot  for  the  open  loop  system  must  be  constructed  with 
a.   =  .167  in  order  to  obtain  the  frequency  response  of  the  open  loop 
system  and  the  procedure  to  construct  this  plot  is  straightforward. 
Since  the  MQ  point  is  so  close  to  the  origin,  the  0)  (.167  -  B.  )  plot 
will  not  give  good  resolutions  of  the  angles  of  the  denominator  for 
small  0d     if  the  plot  is  constructed  in  Fig.  2-29.   This  can  be  overcome 
by  plotting  the  curves  for  small  6)t  on  an  expanded  scale. 


100 


3.   The  fourth  order  system  charts. 

In  references  1  and  2,  a  universal  third  order  system  chart  is  con- 
structed on  the  Bq  vs  B^  plane  so  that  a  single  chart  will  provide  solu- 
tions to  any  third  order  system  problem. 

This  is  done  by  using  a  suitable  transformation.   Any  third  order 

system  is  transformed  into  a  normalized  system  so  that  it  has  unity  co- 

3 
efficients  of  the  terms  S  and 

a  set  of  fixed  Bq  vs  B,  curves. 

In  a  fourth  order  system, 

term  in  its  characteristic  equation,  Bq  vs  B^  curves  are  changed  by  chang- 

3 
ing  either  a-^   the  coefficient  of  the  S   term  or  a£  the  coefficient  of  the 

S2  term. 

Since  for  any  polynomial  it  is  possible  to  make  any  two  of  the  coef- 
ficients unity  (or  the  same  numbers),  by  a  suitable  linear  transformation, 
in  a  fourth  order  system  characteristic  equation  any  two  of  the  coeffi- 
cients of  S   ,  S  and  S   terms  which  are  the  coefficients  affecting  the 
Bq  vs  B^  curves  can  be  made  unity  so  that  all  Bq  vs  B^  curves  of  any 
fourth  order  system  depend  upon  only  one  coefficient. 

For  each  value  of  the  variable  coefficient,  there  results  a  different 
set  of  the  Mitrovic's  curves.   If  the  sets  of  the  curves  for  some  suitable 
pertinent  values  of  the  variable  coefficient  are  constructed  and  if  the 
behavior  of  the  set  of  the  curves  for  the  variable  coefficient  changing 
from  one  pertinent  value  to  others  is  understood,  then  the  sets  of  the 
curves  may  serve  as  fourth  order  system  charts  which  can  be  applied  to 
any  fourth  order  systems. 
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3-1.  Transformation  of  the  fourth  order  system  equation. 

Suppose  a  fourth  order  system  characteristic  equation  is: 

F(S)  =  S^+  &3^  +  a^+  £,£  +  a0  (3-1) 

and  associated  B~  vs  B.  curves  are  defined  by 

Bi  =  ^(pjs^oo^  +•  agfecsycow  +  <bcs)ctx£ 
b0  =  5^  -2/^W  -^(S)co.f 

Using  a  transformation  or  frequency  scaling 
Equation  (3-1)  can  be  written  as 
or  by  dividing  through  by   <3L» 

f{4)-4*+  4*+*+/^  a%+.zj*Ia  +  ^/4*° 

or        i(A)**AU+  xS3  +  3^t4Z  +   a.|t^  4-  a6t  =6  (3-3) 

where      3.^  =  2.2  /  <2.5 

Thus  Equation  (3-1)  is  transformed  into  Equation  (3-3)  in  which  the 
two  coefficients  of  the  S  and  S   terms  are  unity. 

Equations  (3-1)  and  (3-3)  are  equivalent  since  the  transformation 
S  s  a3s  is  linear  and  one-to-one. 

The  Bq  vs  B^  curves  associated  with  Equation  (3-3)  are  defined  by 

(3-4) 
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Equations  (3-2)  and  (3-4)  are  equivalent  to  each  other  with  the  fol- 
lowing relatione. 

Comparing  Equations  (3-2)  and  (3-4)  which  are  Mitrovic's  equations 
for  a  system  represented  on  different  coordinate  scales,  it  is  seen  that 
studying  a  fourth  order  system  with  the  second  equation  is  more  convenient 
than  with  the  first,  because  in  the  second,  only  one  variable  coefficient 
a2t  appears  whereas  in  the  first  the  two  coefficients  a2  and  a3  appear. 

A  fourth  order  system  closed  loop  transfer  function  represented  as 

transforms  into 

to  which  the  Mitrovic's  curves  defined  by  Equation  (3-4)  associate. 

The  coefficients  and  frequencies  of  the  transformed  system  are  related 
with  those  of  the  original  system  by 

plus  the  relations  defined  as  in  Elation  (3-3)  for  the  coefficients  of 
the  denominators. 

The  velocity  constants,  B^  and  K^  of  the  original  and  transformed 
systems  respectively  are  related  by 
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a* 
since      KV=a7K 

Thus  a  fourth  order  system  is  transformed  into  a  system  from  which 
resulting  BQ  vs  B,  curves  are  dependent  upon  only  one  coefficient  a2t. 
Comparing  Equations  (3-6)  and  (2-38a),  it  is  noticed  that  the  characteris- 
tic equations  are  of  different  forms  which  indicates  that  the  normalized 
\  q     curve  in  Section  2-6  is  different  from  the  I o   curve  resulting  from 
Equation  (3-4)  and  any  I  ^  curve  defined  by  Equation  (3-4)  cannot  be 
plotted  on  the  normalized  Bq  vs  B^  plane  unless  the  system  Equation  (3-6) 
is  transformed  into  the  form  of  Equation  (2-38a)  and  Equation  (3-4)  is 
forced  into  a  form  to  give  the  normalized  >0   curve.   If  one  wishes  to  work 
on  the  normalized  Bq  vs  B^  plane  with  l-g  curves,  then  the  system  Equa- 
tion (3-5)  or  (3-6)  is  transformed  into  the  form  oT-Equation  (2-38a)  to 
give  the  equation  defining  Mitrovic's  curves  as 


■On  =  0\t   -  X  tpJSIcdn     -  %(S*)Cdn 
/3at 


(3-7) 


Equation  (3-7)  gives  the  normalized  \q      curve,  and  the  !■$  curves 
can  be  plotted  on  the  normalized  Bq  vs  B^  plane;  however,  /a^t  aPP©ars 
in  too  many  terms  and  more  complicated  behavior  of  the  curves  is  expected. 
Thus  the  normalized  Bq  vs  B^  plane  is  not  suitable  to  study  fourth  order 
systems.  However  the  normalized  curve  can  do  an  excellent  job  in  determin- 
ing absolute  stability  of  a  fourth  order  system. 

The  coefficient  a2t  appearing  in  Equation  (3-4)  indicates  that  Bq  vs 
B^  curves  for  a  fourth  order  system  depend  upon  the  value  of  a2t  , 
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therefore  it  seems  worth  while  to  investigate  the  possible  range  cf  a^. 

The  coefficients  a2  and  a«  of  Equation  (3-1),  the  characteristic 

equation  of  the  original  system,  can  be  written  as 

a3  s  r,+  r^-*-  r5+  r^ 

(3-8) 
a2  z  r,  v  r,  r3+  rir«  +  rxr34'^r<t+-r?r^ 

where  -r^  ,  -r2  ,  -r3  and  -r^  are  the  roots  of  Equation  (3-1). 
The  roots  r^  ,  r2  ,  r3  and  r^  fit  into  one  of  the  following  three  cases  : 

1)  All  4  roots  are  real. 

2)  Two  real  and  one  pair  of  complex  conjugate  roots. 

3)  Two  pairs  of  complex  conjugate  roots. 

2 
Noting  that  a~t  s  &2> ^a3  »  by  using  a2  and  a3  expressed  in  terms  of 

roots,  it  is  possible  to  determine  the  rang 

the  fourth  order  systems  for  Case  1  and  Case  2. 

Case  1 

a  r'  r*  +  n  ^  *  r'  r<*  +  rz  ^  "*"  r2  P*  "*"  ^  r* 

2t  '  ( r,  +  r*  +  r3  +  IV)2- 

r,r2  -»•  r, r3  +  rrr<,  +  r2n  •>- ra r»  +  r3  r» 

noting  that        r,  V  r*  +  r3*  +  V*   >0 

o  <  a2t  <  .  5 

It  can  also  be  seen  that  for  a  fourth  order  system3  if  its  open  loop 
transfer  function  has  a  numerator  polynomial  of  order  less  than  twos  and 
if  it  has  all  real  poles,  then 

a2t  < .  5" 

In  such  a  system,  the  coefficients  a2  and  a3  of  its  characteristic 
equation  can  be  expressed  in  terms  of  the  poles  of  the  open  loop  transfer 

105 


function.  Those  poles  of  the  open  loop  transfer  function  replace  the 
roots  of  the  characteristic  equation  in  Equation  (3-8)  to  get  the  same 
result. 

Case  2 

Let  the  two  real  roots  be  r,  and  r^  ,  and  the  complex  conjugate 
roots  r^  and  r,  . 

r3    =  3  co^  4- j  yr^-i  «>n 

r4     =     3*4,  -j/^*^ 


and 


*\    -  ( r,  +  rz  +  %  s  co^  f=.(rt  +rzf+  o-  5<4,  (r,  +r*0  -»-  4-sX? 
a2   =  r,ru  +2*^(1-, +&)  +  *£ 

to  compare  a^  and  a« 

<£-  aa  =  r,a+  nrj*r*V2^«j,(r,+r^  +(*^-  t  K*  0-9) 

Assuming  a  stable  system,  Equation  (3-9)  is  examined  for  several 
situations. 

If   rL  }  r2  >  S  uXn  then 

<£-  22  >  (l  1  3*-  0  &£  implies  that 

if    5>.3o2      then  2.*-   2a  >  Q  and   ^t  <  LO 

It  must  be  noted  that  the  inequality  above  may  have  considerable  tol- 
erance which  implies  that  even  for  %    some  value  less  than  .302  a2t  may 
be  less  than  1. 

Let  r2  >  r^  and  by  neglecting  r^   in  Equation  (3-9) 

4'^>  £  +  2lU>*rx+   (4*2-0<a  (3=10) 

For  a  well  behaving  system,  r2  is  usually  greater  than  ^oX^    the 
real  part  of  the  complex  conjugate  roots. 
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Let  r2  s  k^oO^   where  k  is  a  positive  real  number  greater  than  1. 
Substituting  this  r2  into  inequality  relation  (3-10) 

Inequality  relation  (3-11)  indicates  that  for  a  reasonable  number  k,  if 
^  is  not  too  small,  then  a|  -  a2  is  greater  than  0  and  a2t<1.0. 

For  example  if  k=2.5  then  a^  -  a    >  0  for   £>.25  and  if  k  =  3 
then  a|  -  a2  >  0  for   S  > . 229 

Thus  for  a  fourth  order  system ,  if  it  has  no  more  than  one  pair  of 
complex  conjugate  roots  of  the  characteristic  equation  and  if  it  is  ade- 
quately damped,  then  the  coefficient  a?  is  less  than  1. 

For  Case  3,  which  is  the  case  of  two  pairs  of  complex  conjugate  roots 
present  in  the  characteristic  equation,  it  is  hard  to  set  up  an  upper 
bound  of  a2  as  in  Case  1  and  2.  This  can  be  seen  as  follows  % 

For  Case  3 
Let  fl"f*    )      k-^V 

and  rj  -  ^COm  -j  /i=5»flXni       ,      lj  -  % Ol„a  - j fi-g  ^ 

then  4^(2:5,^  4  2§i%) 

&z  =  CUaf  +  U.  %%  tO* to**  -+-  ^a. 

Let  0^,  <  0)^,2.  a?)d  £=  ^  /^, 

a*  i  -*-<^  ?*£  +  £* 

then  —j-  ae    5 5 — - — 
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Thus  the  ratio   2/^2  depends  upon  %,   ,    S>  and  k  .   If  both  $, 
3  , 

and  5o  are  small,  then  2/  2  assumes  a  large  number.   However s  if  a  well 
a3 

damped  system  is  assumed,  then  Case  3  can  be  approximated  to  Case  2  and 

&2t   "^  1  is  a  reasonable  upper  bound.  For  example  if  $^   =  .3,  ^  :  .7 

and   <%,=<0»a    then  ?2  =  2,84  =  #  n 
aZ     4.0 

3 

In  this  section,  in  studying  fourth  order  system  curves ,  a2t  <  1  is 
assumed.   It  must  be  noted  that  a«t  <  1  is  not  the  sufficient  condition 
for  a  fourth  order  system  to  be  well  damped.   Even  though  a2t  is  very 
small,  a  fourth  order  system  may  be  unstable  depending  upon  the  location 
of  the  M  point.   It  turns  out  that  one  can  find  a  necessary  condition  on 
a2t  ^or  a  system  to  nave  a  relative  stability  greater  than  ^.   a  given 
system  damping  coefficient.  This  will  be  discussed  later. 
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3-2.  Three  dimensional  representation  of  fourth  order  curves. 

The  Bq  vs  B^  curves  for  fourth  order  systems  defined  by  Equation  (3=4) 

2 
depend  upon  a2t  which  is  the  coefficient  of  the  s  term  in  the  character- 

istic  equation  of  the  transformed  system  equation,  a.  can  assume  any 

value  within  its  allowable  range  and  independent  of  frequency  on  the  B0  vs 

Bi  plane. 

If  a  new  variable  B2t  =  a.  is  defined  and  if  B2fc  varies  on  B2t  axis 

which  is  chosen  to  be  perpendicular  to  both  B.   and  B.  axes  to  form  a 

three-dimensional  cartesian  coordinate  system,  the  equations 

ftzt  =  d-zt 

B,t=^t&W^  +  ^)<^"  +  tfW-S)*^  (3-12) 

Bot  =  a^  -  &(S)c4J!  -  %(*)">»* 

define  the  Mitrovic's  surfaces  on  the  B2t  -  B.   -  BQt  coordinate  system. 

For  fourth  order  systems.  Equation  (3-12)  defines  fixed  surfaces 
which  can  be  applied  for  any  fourth  order  system  as  the  case  of  the  third 
order  chart  for  third  order  systems. 

Since  the  variable  B2t  is  continuous  and  BQt  vs  B,..  curves  are  con- 
tinuous, the  Mitrovic's  surfaces  are  also  continuous.  Noting  that  any 
plane  B2t  z  a2't  =  constant  which  cuts  the  surfaces  gives  the  set  of  Bqc  vs 
Blt  curves  f°r  various  values  of  <£  ,  by  applying  the  continuous  property 
of  the  Mitrovic's  surfaces,  all  of  the  basic  property  of  the  BQ  vs  B, 
curves  and  the  analysis  and  design  principles  of  feedback  control  systems 
on  the  Bq  vs  Bj^  plane  can  be  extended  to  the  three-dimensional  Mitrovic's 
surfaces. 

For  example,  the  stable  region  for  the  M  point  on  the  Bq  vs  B]_  plane 
which  is  the  region  in  the  first  quadrant  enclosed  by  the  1^  curve  can  be 
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extended  to  the  three-dimensional  Mitrovic's  surface  as  the  region  (volume) 
in  the  first  octant  enclosed  by  the  V0   surface  which  is  the  surface  made 
up  by  the  \  0   curves  for  all  values  of  a2t  ^  0«   Locating  the  point 
M(a  ,  aQt)  on  B0(.  vs  Blt  plane  is  to  locate  the  point  M(a2  ,  a,t,  aQt)  in 
the  Mitrovic's  volume. 

Fig.  3-1  is  a  schematic  sketch  of  the  l0  surface.  B2t  axis  is  chosen 
so  that  B2t,  B^t  and  BQt  axes  form  a  right-handed  coordinate  system.  The 
equations  defining  the  F^  surface  are 

&lt=  Cd£  (3=13) 

which  are  the  equations  for  a  parabolic  surface.  The  surface  forms  a 
cave  with  closed  end  at  the  origin  on  the  B«t  ■  0  plane.  Applying  the 
stability  criterion  to  each  BQt  vs  B^t  curve  which  is  formed  by  cutting 
the  surface  with  each  B2t  =  &ot  Plane»  tne  stable  region  is  the  whole  vol- 
ume inside  the  cave. 

If  the  surface  is  cut  by  a  plane  Bq.  a  a,*.  =  constant ,  then  there  re- 
sults a  B, t  vs  B2t  (or  B  vs  B  )  curve  fo  •  £  =  0  on  the  plane  and  like- 
wise for  the  BQt  vs  B2t  curve. 

It  is  interesting  to  note  that,  Bgt  vs  B£t  curves  for  $  =  0  on  BQ. 
vs  Bot  planes  for  fourth  order  systems  are  straight  lines  whose  slopes  and 
axis  intersections  are  determined  by  a-^t  only  and  an  extremely  simple 
method  to  test  the  absolute  stability  of  a  fourth  order  system  can  be 
derived. 

By  eliminating  Od  .  and  a2  in  Equation  (3-13),  the  \Q  surface  of  a 
fourth  order  system  is 

Bot^e.t-Brt  C3-14) 
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,  /  B  =constant  plane 


Fig.  3-3   The  I   surfece  for  4th  order  systems. 


By  cutting  the  surface  represented  by  Equation  (3-14)  wttb  a 
Blt  =  alt  Plane>  tne  B0t  vs  B2t  curve  f°r  $   -  0  Is  obtained  as 

30t=  a.ti^t-a,*  (3-i5) 

which  is  an  equation  defining  a  family  of  straight  lines  on  the  BQt  vs 
B2t  plane  (or  B^t  r  alt  =  constant  plane). 

For  each  value  of  a,  ,  there  results  a  different  line.   By  reason- 
ing with  the  stable  region  in  the  three-dimensional  volume,  the  region 
for  the  point  M  (a2t  *  aQt)  on  BQt  vs  B2t  plane  to  guarantee  a  stable 
system  is  the  area  in  the  first  quadrant  below  the  line  defined  by  Equa- 
tion (3-15)  or  mathematically 

o<  aot  <art.  axx  -2-,* 

which  agrees  with  the  Routh  criterion. 

The  r©  curves  for  some  values  of  a,   are  plotted  in  Fig.  3-2.   As 
the  Equation  (3-15)  indicates,  a  f^7  line  for  a   =  k  cuts  the  B2t  axis 
at  B2t  =  k  and  it  has  a  slope  of  k.  The  frequency  on  these  1"^  lines  is 
constant  at  all  points  on  a  line  which  is  a  peculiar  nature  of  the  f7 
curve  on  a  BQt  vs  B2t  plane.  Thus  the  f^  line  itself  is  a  constant  0)nt 
curve  on  the  Bqj.  vs  B2t  plane. 

For  instance,  along  the  \q   line  for  a^fc  s  .2,  0)   =  .446  which  is 
equal  to  /SL.x       as  tne  second  of  Equation  (3-13)  implies.   Thus  a  point 
M  (aQt  ,  a2t)  located  anywhere  on  a  single  \0   linSj  say  for  &,.  ~  a  '  , 
represents  a  closed  loop  system  which  has  a  pair  of  pure  imaginary  roots 
r  =   ±  (Ct)    -   ±  i'  -/aiV  which  is  readily  seen  from  the  basic  nature  of  the 
1^  curve. 

As  an  example,  suppose  a  characteristic  equation  of  a  fourth  order  system 
is 
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Since  a,   =  .2  the  1^  line  is  the  line  labeled  as  a.   r  .23 
COt   =  .446  in  Fig.  3-2.   The  Mitrovic's  working  point  M  (a0t  ,  a2t)  lo- 
cated anywhere  above  the  i^  line  represents  an  unstable  system^  and 
dually  for  below  the  line.   If  the  M  point  is  followed  along  the  line 
then  resulting  every  characteristic  equation  has  a  pair  of  imaginary 
roots  r  =  ±j  .446.   For  instance,  if  M  (.7,  .1)  is  chosen,  then  the 
characteristic  equation  is 

and      -f,4.><^)=0 

The  point  M  (.7,  .1)  is  also  on  the  \Q   line  for  a,   -  .5,  therefore 
another  characteristic  equation 

has  the  roots  r  :  ±i'/5    and  it  calculates  to 

&#/*>«  0 

From  the  Mitrovic's  basic  equation  and  the  way  of  defining  a  three- 
dimensional  plot,  the  concept  of  a  Mitrovic's  surface  can  be  generalized 
to  any  order  system  and  to  a  set  of  any  three  coefficients  of  the  charac- 
teristic equation.  Then  fixing  one  of  the  three  axes  to  a  constants 
there  results  a  well-defined  curve  which  is  the  one  defined  by  the  gener- 
alized Mitrovic's  equation  in  Reference  3,   though  the  generalized  Mitro- 
vic's equation  turns  out  to  be  undefined  and  uninterpretable  at  some 
particular  values  of  ^  for  some  particular  set  of  coefficients,   (For 
example,  the  f^  curve  for  the  BQt  vs  B2t  plot.) 

The  most  significant  feature  in  the  concept  of  the  Mitrovic's  sur- 
face is  the  possibility  of  controlling  three  coefficients  simultaneously 
since  the  surface  is  independent  of  the  three  coefficients  which  are  to 
be  controlled. 


3-3.  Study  of  fourth  order  curves. 

As  pointed  out  in  Section  3-1,  a  fourth  order  system  can  be  trans- 
formed into  an  equivalent  system  so  that  there  remains  only  one  param- 
eter (coefficient)  in  the  characteristic  equation  on  which  the  fourth 
order  system  curves  depend.   This  feature  greatly  simplifies  studying 
the  behavior  of  the  Mitrovic's  curves  for  fourth  order  systems. 

In  Section  3-2,  it  was  shown  that  if  a  fourth  order  system  charac- 
teristic equation  is  transformed  such  that  the  coefficients  of  the  S 
and  S   terms  are  unity,  then  the  coefficient  of  the  S2  term  a2t  which  is 
the  only  one  affecting  Bqj.  vs  B^t  curves  is  less  than  1.0  for  most  prac- 
tical fourth  order  systems. 

Since  the  Bqj.  vs  B^  curves  are  continuous  for  all  values  of  Cd    .  s 
^  ,  and  the  coefficients  of  characteristic  equation,  if  the  behavior  of 
the  curves  is  well  investigated  and  if  the  curves  for  some  pertinent 
values  of  parameters  are  constructed,  then  those  precons true  ted  curves 
may  serve  as  a  guide  or  charts  for  both  analysis  and  design  of  fourth 
order  systems 

In  this  section,  the  effect  of  ayt   ©n  B,>t  vs  B*t  curves  is  studied 
to  provide  the  fourth  order  charts. 

In  Equation  (3-4)  which  is  the  equation  defining  B~  vs  B,   curves 
for  a  fourth  order  system,  there  are  two  parameters  a2t  and  5  •   F°r  one 
value  of  aoj.  ,  each  ~£  gives  a  different  curve  I  ^  to  form  a  family  of 
curves  corresponding  to  the  particular  value  of  a^.   This  family  is  de- 
noted by   Ifeh3*t    »  or  simply  \{i}        .   Then,  for  each  a2t  there  re- 
sults a  different  family. 

Fig.  3-3  to  Fig.  3-22  show  such  families.   On  each  sheet  the  curves 
for  five  values  of  *>  are  shown.   For  one  value  of  a2ts,  there  are 
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Fig.  3-3  BQt  vs  Blt  curves(4th  order  systems) 


Fig.  3-9  B_,  vs  B-,  curves(4th  order  systems) 
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Fig.  3-10  BQt  vs  Blt  curves(4th  order  systems) 
for  a2t=.4  (^=.5, .6, .7, .8, .9,1.0) 
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Fig«  3-18  BQt  vs  Blt  curves (4th  order  ayntema) 
for  a2t=.8  ($=.5#.6,.7,.8,.9,1.0) 


for  a2t=.9  (S=0.,.l,.2,.3,.4,.5,.6,.7,.8,.9,1.0) 
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Fig.  3-20  B   vs  B   curves (4th  order  systems) 
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associated  two  sheets  of  graphs,  one  for  0  ^  "5  ^  «5  and  the  other  for 
.5  ^  3  ^  1.0  with  interval  of  <,    .  1.  The  values  of  a2t  are  from  .1  to 
1.0,  thus  for  every  two  sheets  of  graphs,  a2fc  increases  by  interval  of 
.1.   (Note  that  for  .4  £  a2t—  1.0,  the  first  graph  for  each  value  of  a2t 
contains  the  family   l$$]       for  0  ~  $  ^  1*0  and  the  second  graph  con- 
tains the  family    l^  ;   for  -5  £  £  £  1.0.)  The  constant  CO   curves 

are  also  shown  for  some  pertinent  values  of  0)      . 

nt 

It  is  observed  that  as  a2t  increases,  the  curves  are  in  general  ex- 
panded in  both  abscissa  and  orxinate.   This  is  well  expected  from  the 
equations  defining  Bgt  vs  B^t  curves.   The  equations  are  rewritten  below 

B,t  =a2t03rs)cowt  +  &<^)^*t  +<zw$)^t 

a  y  (3-4) 

The  terms  with  which  a2t  contributes  to  B^t  and  BQt  are  non-negative 
for  0  ^^>  —  1.0,  thus  increasing  a2t  increases  B^t  and  Bgt. 

For  small  values  of  a2  each  member  of  the  family   li-^i   is  so 
shaped  that  the  whole  region  under  the  curve  in  the  first  quadrant  is 
the  region  for  the  M  point  where  the  relative  stability  is  guaranteed  to 
be  greater  than  the  value  of  £  of  the  curve.   Such  a  region  is  denoted 
by  R^  .   For  example  in  Fig.  3-3,  the  whole  area  enclosed  by  the  I  5 
curve  and  B^t  axis  in  the  first  quadrant  is  the  region  R5  on  which  if 
the  point  M  (a,  ,  a_  )  is  located,  then  the  system  has  "p  greater  than  .5. 

As  a2fc  increases,  the  manner  in  which  the  R^  region  is  formed  is 
changed  in  such  a  way  that  the  curves  for  high  ^  do  not  form  R^  . 

The  behavior  of  the  curves  depends  upon  the  values  of  a2t  and  ^  . 
Examination  of  Equation  (3-4)  reveals  that  for  a  particular  value  of 
a?  ,  there  is  a  certain  critical  value  of  ^   such  that  if  ^  is  increased 


beyond  the  critical  value,  then  the  '  ^  curve  exhibits  counter-clockwise 
encirclement  on  the  Bgt  vs  B-^t  plane  so  that  there  is  formed  a  small  R^ 
near  to  the  origin  or  no  R^  at  all.   Equivalently  for  a  particular 
value  of  ^  there  exists  a  critical  value  of  a2f 

Investigation  of  the  critical  value  of  ^>  or  a2t  is  important  since 
in  synthesis  of  a  system  one  must  choose  a  point  M  (a^ts  aQt)  so  that  the 
specified  system  performances  are  met  and  particularly  for  the  relative 
stability  of  the  system,  the  M  point  must  be  located  on  a  R^  where  ^ 
is  some  suitable  value  to  meet  the  specified  relative  stability.   For 
example  if  the  system  relative  stability  is  specified  by  ^  s  .5  then 
one  should  not  choose  the  M  point  anywhere  on  the  I  ^5  curve  for 
a2t  =  1.0.   (See  Fig.  3-21  and  3-22.) 

In  synthesis  of  a  system,  usually  ■£  is  chosen  at  some  pertinent 
value  according  to  the  required  relative  stability.   Therefore  it  is  con- 
venient to  investigate  the  critical  value  of  a2t  for  given  values  of 
-£>  ,   thus  a  study  of  the  behavior  of  the  curves  with  changing  a2t  for  a 
fixed  £  is  needed.   This  requires  that  the  family  of  curves  must  be 
represented  in  a  different  fashion  from  those  of  Fig.  3-3  through  Fig. 
3-22  in  which  the  families  of  Bgt  vs  B^t  curves  with  changing  5  f°r 
fixed  a   (the  family      |  f$fLa2t)  are  represented.   That  iss  the  family 
of  curves  must  be  represented  with  changing  a2t  for  fixed  ^  .   Those 
families  are  denoted  by        £{a2t}  * 

At  this  point,  in  connection  with  the  study  of  the  families 
I  £j{a2t\  »  a  constant  £  plane  is  defined  as  follows; 

A  constant  ^  plane  is  a  Mitrovic's  plane  on  which  Mitrovic's  curves 
for  fixed  5  are  plotted,  and  denoted  by  the  ^  =  k  plane  where  k  is  a 
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constant  and  0  ^  k  ^  1.0.   With  this  definition,  for  the  fourth  c 

systems,  the  family  of  Bqc  vs  B,.  curves  with  various  values  of  £2*-  ^! 

""$  fixed  to  zero  consists  of  all  of  the  Bqc  vs  B^t  curves 

the  "$=0  plane,  and  likewise  for  any  other  values  of  0  ^  ^  £ 

The  constant  6)n  lines  on  the  constant  •$  plane 

In  general,  if  the  Mitrovics  curves  are  plotted  on  a  constant  ^ 
plane,  and  if  any  one  of  the  coefficients  which  appears  i 
defining  the  Mitrovic's  curves  is  varied,  then  all  of  the 
the  frequencies  are  the  same  are  on  a  single  straight  li 
the  coordinate  axis  intersections  are  determined  by  <?  and  U) r>. 
line  is  called  the  constant  0>n  line.   (Similarly  the  constant  Ct) 
for  a  transformed  system.) 

For  Bq  vs  Bi  curves  defined  in  general  form  as  in  Equat 
by  fixing  £,  z    "£  1  ,  a  constant,  the  equation  of  Bq  vs  B^  curves 
Sj_  plane  are: 

If  the  coefficients  aiS  ,  k=j^  i  are  fixed  and  if  a,  is  varied, 
then  on  the  S 1   plane,  a  family  of  Bq  vs  Bi  curves    |§, £3-^     is  r 

Rewriting  Equation  (3-16)  to  separate  the  term  containing  the 
able  coefficient  a.  from  the  terms  of  the  fixed  coefficients. 
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i  =  ^ft(V>^  +ILd^(^a)7f 


(3-17) 

TO  ^ 

Suppose  any  two  members  of  the  family    >  t^i^j}        »   one  for  a.^a,, 
and  the  other  for  a.  =  a   .   Let  the  increment  in  a.  be  ziaor  a   -a.,  , 
then  the  increments  in  the  coordinates  B,  and  Bq  at  any  frequency  asso- 
ciated with  A>a.  are 

4  B0=-da; &-/<$,)"£  (3-i8) 

^  p, 

ition 

aj  the  ratio  of  the  changes  in  the  coordinates  ZIBq/^B^  is  independent 
of  a  at  all  frequencies,  and  that  the  distance  between  any  two  points 
at  the  same  frequencies  one  on  the  curve  for  a   and  the  other  on  the 
curve  for  ai?  is  linearly  varying  with  A a^  . 

For  the  family     I  £,{2^     on  tne  ^i  plane,  consider  a  point  at 
C0n  :  0)  i   on  each  member  of  the  family.   If  a.  assumes  a  continuous  set 
of  values,  then  the  locus  of  the  point  of  0)^   s  0)   ,  on  the  £  ,  plane  is 

obtained  by  eliminating  a.  in  Equation  (3-17)  and  replacing  0)     by  q)       , 
i  n      nl 

namely 


'^fa-^w<*^)-£$tw*k     <3~19> 


B=_  4U*.)0>*  s+t^*(  WW-WW*   )   (3-19a) 

Since  A   is  fixed  to  ^  i  on  tne  ^i  plane,  and  <D  is  fixed  at  Ld    . 
x  y  \  '  1  n  nl 

Equation  (3-19)  describes  a  line  with  a  slope  of  —  -^tiUjl^HL-  pass 

'  '■&&  f&i 

on  the  £>-.  plane.   Thus  the  equation  defining  the  constant  6in  line  for 

Qn   -     Ld   i   on  tne  *?i  plane  with  carying  a„  is  obtained, 
n     nl  J-  i 

The  axis  intersections  of  the  constant  cd     line  are  found  by  letting 

n 

B^  =   0  and  BQ  r  0  in  Equation  (3-19)  or  (3-19a) 


ing  through  the  point  (  |I  a^(S,)  COJf _1  ,  _|I  a^/S,*)  0)^  j 


Bl  ^tersect.  .-  £  ^£C^a£L  J^  |r*/V^'>  ^ 


«K 


s£  a  ^ /A-,f*)ft«)-*^ft..<<.>\  (3.20) 


B0  Ut.r.ee,  _-|aX  (^-^^2^1)       (3-21> 

In  particular  if  i  =  2  then  noting  that  -(&04.,-4k-z  ~  $£ 
and  fa  z   -1,  Equations  (3-19),  (3-20)  and  (3-21)  reduce  to 
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tV  ^;i(Ei+Z.V4,--.-l;  )  0-22) 


p:,  ,,.,,,, ■■...=  -21  a,  ,.  , 


For  a  fourth  order  system  whose  Bg  vs  Bi  carves  are  defined  by  Equa- 
tion (3-4)  which  is  for  the  transformed  system,  the  constant  6)    t   lines 
for  a)nt  r  0)ti  on  the  S,  plane  with  a2t  varying  is  obtained  from  Equa- 
tion (3-19)  or  (3-22)  by  proper  substitutions,  thus 

BOt  =  ^5i)f.Blt-CUt*(l-«($I)cOtli)  (3^25) 

with  axis  intersections  given  by 

Blt  intersection  =  G)t*  (  \   -cp_.it,  ):'-zi)  (3-26) 


u0t 


intersection  =-  2&-   /  j-  />(:,;i<3f;)  (3-27) 


Equation  (3-25)  defines  the  constant  u)nt  lines  for  all  values  of 

(a)  t  on  the  $^  plane  for  fourth  order  systems.   The  slopes  of  the  cc  . 
stant  &)_t  lilies  are  proportional  to  the  frequencies  and  inversely  propor- 
tional to  $m    (  0^\  ,  ~     2*).      For  5  r  0,  the  slopes  3 
all  constant  ®nt   lines  on  the  ^  =  0  plane  are  vertical  lines, 
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In  view  of  Equations  (3-26)  and  (3-27),  if    6)tl  <  ^ 

Blt  *ntersection  >    0 
Bgt.  intersection  <  0    and  dually. 
For  the  maximum  values  of  the  B^t  intersection  an  I 

BQt  intersection,  differentiating  Equations  (3=26)  and  (3 

Max.  Blt  intersect.  =   1/27$*    at    a)tJ  =  I /3  $j 
Min.  B0t  intersect.  =  -27/2048  $*         at  u)ti^3/<?<?± 

Thus  the  constant  0JRt  line  of  the  fourth  order  syst 
right  with  increasing  slope  as  63nt  is  increased.   Then  at  a)  t  =  1/3$ 
the  maximum  3^t  intersect,  is  reached.   As  the  frequency  CO  ^  Ls 
beyond  1/3 5>  then  the  constant  k)nt  line  moves  back  to  the  lefi 
steadily  increasing  slope. 

Fig.  3-23,  3-24  and  3-25  are  families  of  B0t  vs  Blt  cur 
L  ;a  1      plotted  on  the  £  =  0,  .1  and  .2  planes  respectively. 
coefficient  a2t  is  varied  from  «1  to  1.0  with  A&2t.   *  •  *  •  ^''-~   '  instant 
6d  nt  lines  are  shown  some  pertinent  values  of  frequencies s  end  those  can 
be  compared  with  the  constant  &-)nt  curves  in  Fig.  3-3  to  Fig.  3-22 
are  the  families      '"fe^alt  * 

The  property  of  the  constant  co   line  can  be  used  in  constructing 
a  B«t  vs  Bi  t  curve  on  a  constant  £,  plane  for  any  value  of  «,,  whet! 
are  on  the  plane  any  two  curves  each  for  a  known  value  of  a2t«   TM 
illustrated  in  Fig.  3-26. 

Suppose  two   I -^  curves,  one  for  &2t  -    «30  and  the  other  for 
a£t  =  .40,  are  precons  true  ted  on  the  *?  =  .5  plane.  The  line  sej 
joining  the  pairs  of  two  points  at  the  same  frequency  on 
are  the  constant  60nt  lines.   For  example,  the  points  P^  and  I  ?  are  at 
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Fig.  3-24  B   vs  B   curves  oh  the  £=.1  plane. 
(4th  order  systems) 


144 


\:.. ;:~\vi: : 


Fig.  3-?5   B   vr,  3   curves (4th 
on  the   $=.2  plane. 


^Bxia 


WwWF^VMEM 


mmmm^m^M 


.D06- 


,002. 


Fig.  3-26   Illustation  of  constructing 
a  BQt  vs  B   curve  by  using 

constant  CO  .  lines. 
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the  same  frequency  q)       z    .  25  on  the  two  curves,  one  for  ^t  -  •  30  s 
the  other  for  a2t  z    .40  respectively,  end  the  line  segment  joining  the 

points  P^  and  ?2   denoted  by  L  25  *-s  fche  constant  0)nt  line  f or  &)  t  = 

j — ' 
All  members  of  the  family    l-5/^tlr   will  have  the  t0nt  -  .25  points 

or  extension  of  the  line  segment  L^S  • 

Suppose  the   I  5  curve  for  a2(_  =  .32  is  to  be  constructed.   Si 
the  increments  ^B^t  and  <^B2t  associated  with  ^^j-  are  proportional  tc 
•A  a2t  at  every  frequency,  the  point  Q  which  is  located  at  one- fifth  way 
from  the  point  P,  to  the  point  P?  on  the  constant  '0at  line  L  95  is  the 
point  at  CC't  -    .  25  on  the   l#g  curve  for  a2t  =  .32.   Any  other  poir»,  s 
the  curve   I  5   for  a2t  =  .32  are  determined  analogously,  then  connecting 
those  points  with  a  smooth  curve,  the  desired  curve  is  obtained  as  shown 
by  the  dotted  curve. 

With  the  aid  of  the  constant  63  t  line,  one  can  visualize  the  changes 
of  BQt  vs  Bi t  curves  with  changing  a2t  .   For  example  a  point  P^  at 
6)nt  =  .25  moves  along  the  line  L  25  *-n  Fig.  3-26  as  a- t   increases  from 
.3.   If  a  point  is  located  on  a  constant  CO   line  for  Cu  fc  =  6),  on  the 

b  =  ^  plane,  then  the  value  of  a2t  for  which  the  1  5^  curve  passes 
through  the  point  can  be  determined  from  the  linearity  property  of  the 
constant  0)    .    line. 

Now  consider  an  M  point  located  at  (x)nt  =  63^  on  the  -j  curves 
and  one  asks  if  the  M  point  insures  ■--■  ?"  >,  for  the  system,.  For  the 
problem  to  be  specific,  suppose  the  M  point  is  chosen  at  '^j-  =  .5  on 
I  t5  curve  for  a2t  =  1.0.  The  family  I  ,5,-{a«Ji  is  shown  in  Fig.  3-27. 
By  applying  the  basic  rule  to  determine  system  stability  to  the  point 
M  (.375,  .125)  with  the  l<5  curve  for  a2t  :  1.0,  it  is  obvious  that  the 
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system  has  a  pair  of  complex  conjugate  roots  with  <£  less  than  .5. 

It  is  seen  in  Fig.  3-27  that,  for  a£t  less  than  or  equal  to  .  7S  the 
entire  areas  enclosed  by  the  1 . 5  curves  in  the  first  quadrant  are 
R5  region.   For  a2t  =  »8  and  .9  the  areas  enclosed  by  the  curves  are 
divided  into  two  regions,  one  the  R^  region  and  the  other  the  nc 
region.   For  a2t  =1.0  the  curve  gives  almost  no  R  5  region. 

The  curves  in  Fig.  3-27  agree  with  the  curves  in  Fig.  3-13  t 
For  £  =  .5,  the  critical  value  of  a2t  (mentioned  previously^  seems  te 
lie  between  a2t  =  .7  and  a2t  :  .8.   If  this  critical  value  of  su*  is  de- 
termined, then  the  M  point  chosen  at  any  point  on  the  I  ^5  curves  for  a«t 
less  than  the  critical  a2t  will  guarantee  the  system  £  greater  than 
equal  to  .5  whenever  the  M  point  is  in  the  first  quadrant.   In  the  fol- 
lowing paragraphs  the  critical  value  of  a„  denoted  by  a2tc  is  investi- 
gated for  each  "£  . 
The  critical  value  of  a2t. 

Inspection  of  the  figures  shows  that  in  order  for  a  I-5  curve  tc 
form  a  region  R^  ,  the  curve  must  encircle  a  point  in  the  region  R   lis 
the  clockwise  direction.  This  implies  as  one  possible  situation  that 
the  points  of  the  extrema  of  the  coordinates  B.   and  B~  must  occux  in  a 
definite  order;  the  maximum  point  of  BQt  followed  by  the  maximum  point 
of  Bifc  (and  then  minimum  B~  followed  by  minimum  B.  ). 

Let  &\max  and  (^>Qmax     he  the  frequencies  at  which  Blc  and  Bgt  are 

r->  r~7 

maximum  for  a  I -c   curve.   Then  in  order  for  a  L:  curve  to  form  a  region 

R   ,  the  condition 

^Omax  <-        ^lmax 

must  be  fulfilled.   (Note:   even  if  the  condition  60«    <  &1    is  not 
N  Omax    Imax 
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fulfilled,  there  may  exist  R   and  this  is  considered  later.} 
illustrated  in  Fig.  3-28.   In  Fig.  3-28a,  ^Qmax  <   a'l' 
R.  is  formed,  whereas  in  Fig.  3-28b,  Cjl) Omax  >  ^Imax  and  there  is 
region  R  .   In  Fig.  3-28c  the  critical  case  ^Qmax  =  Climax  -  - 
and  the  region  R,  Is  formed. 

Since  uJq^j,  and  <%max  are  functions  of  a2t  and  ^  in  the  fourth 
order  curves,  the  conditions  C^q^x  ~    ^Imax  can  be  exPressed  in  '"- 
of  a2t  and  ^  thus  for  a  given    a^    the  critical  val 
deterained. 

The  values  of  6)j.max  and  ii)oraax  are  f°und  bY  differentiating  the 
equations  defining  the  Bot  vs  Bit  curves  Equation  (3-4)  which  is  rewrit- 
ten below. 

30t=aAt^t  -Acsx^t  -<fccs^t 

Noting  that     &(i)=   0  at  £  =  0,     v§($)  =  0  at  ^  --    .5   &<$,) 
:  0  at   £=  0,  and  .707,  in  determining  the  conditions  &^r.iax  =  &) 
it  is  convenient  to  divide  the  interval  0  £  $    ~  1.0  into  3  sub- 
intervals  0  £  S  *  .5,   .5  <  S^l  .707  and  .707  <  £  <  1.0.  The  spe- 
cial cases  £  =  0,   .5,   .707  and  1.0  are  considered  separately. 

Differentiating  Blt  and  Fqc  with  respect  to  (50 nt  in  Equation  (3-4) 
and  setting  to  zero 

(3-28) 

! 


"ot  rr\<x*. 


B«* 


(b) 


Bit  wax 
cdi-max 


co±  iNC.rea.siN3 


a)eYnax=6)nnax 


Fig.  3-28  An  illustration  of  behavior  of 

B0t  vs  Eit  curves. (The  region  R^  and  extreme  points) 
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If     Cfe&O     and      (prfO     ,    then  by  solving  each  of  Equate-    (3-28) 
for   0)nt 


(3-29) 


where  CO   is  the  frequency  at  the  extremum  of  B^t  and  like- 
wise for  ^Qx* 
For  the  interval  of  $   ,  0  <  %<  .5,  4>  2  >   0,   <A  >  0  and  A<  0. 


Therefore 


01  -^Lj^L-fk* 

lraax  "     3<fr  tj  ^       3  ^ 


(3-30) 


^Omax 


3^       Irtf       22±l: 


Note  that  the  negatives  of  the  double  signs  in  Equation  (3-29)  give 
Chimin  an<^  ^Orain  an<*  tnose  are  negative  numbers,  and  also  in  Equation 
(3-30)   &>imax  and  &0max  are  positive  numbers. 

Next  one  must  find  the  conditions  on  a£t  which  are  equivalent  to  the 
conditions  ^Qmax  ==  ^lmax*   This  is  done  conveniently  with  the  aid  of 
Equation  (3-28)  instead  of  imposing  the  conditions  ^Qmax  ==(0ima 
Equation  (3-30). 

Since  ^lmax  and  ^Omax  are  solutions  of  Equation  (3-28),  the  sub- 
stitutions Gu>nC  =  ^lmax  and  ^nt  =  ^Omax  *-nto  tlie  f*-rst  B^-   fc^e  se 
of  Equation  (3-28)  respectively  must  satisfy  the  Equation  (3-28), 
thus      a2t<fizt-  24>?ti)/7?i<i/:   -t  3<P*(c0fw&*)    =0 
2a2t  -  3  4>a  il\  ,>  a  %-4.<fy  (a)0-m  i<xf  **  C 
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solving  for  a2t  and  completing  the  squares 


(*> ^m 


l^L 


v2.    ,3,25 


a»t-«*  (««..« +  -^)--^ 


Equation  (3-31)  gives  the  relation  between  a2ts  ^imax  ar*d  ^Qmax*  -^t~ 
ing  the  signs  of  Cp2i    Cp^   and  ^4,   and  -$$!*    >2$  for  0  ^  <f, 
Equation  (3-31)  which  ire  parabolas  can  be  sketched  as  in  .   ,    29. 

The  two  equations  of  Equation  (3-31)   are  plotted  on  the  same 
ordinate  system,  a2t  on  the  ordinate  and  ^imx  »  ^Omax  cn  fc^e  a^:' 
Both  parabolas  pass  through  the  origin  and  are  concave  upward  since 

-3^4/^2  >  °»   2  ^3  >  0  for  0^  £<  .5.  The  vertices  of  the  parabolas 

<#- 
are  below  the  frequency  axis  for  both  since     — ,J.        <     0  ai  1 

-2,2b  4>* /S  d>3  <    0  and  they  occur  at  63Lmax  =  ~^/^^>0  and 

^Omax  =-/.?^/^  <  0.   Thus  the  conditions 

^Omax  =  ^lraax 
can  be  expressed  in  terms  of  a2t  and  a2tc  as  eL       ==  3.  tc 
where  a2tc  is  the  value  of  a2t  at  which  ^Omax  ~  ^Itrtax*   3^  letting 
60tc  =  6)1  max  =  ^0max>  Equation  (3-31)  forras  a  set  of  two  equations 
with  two  unknown  to  give 

(3-32) 

ft)  - . — _ 

tC    »$(!  +  $ +  $a) 

By  proceeding  as  above,  it  is  easy  to  see  that  Equation 
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(The  1st  of  Eq.  5-31 )  (  /  (The  2nd  of  Eq.  3-31) 


Fig.  3-29  The  relation  between  a,^  ,  <a.    and  <J 

2tc  ~max      Omax 
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for  the  intervals  .5  <    $  <    .707  and   .707  <g  <  1.0.   This  is 
Equation  (3-31)  holds  for  any  of  the  3  intervals  of  ^  provided  0)^^ 
and  t^Omax  defined  by  Equation  (3-30)  exist.   Therefore  for  the  two  in- 
tervals of  $,  .5  <  $  <  .707  and   .707  <-  $  <  1.0  it  is  necesc . 
check  if  the  value  of  a2tc  determined  by  Equation  (3-32)  gives  positive 
real  Cx\mSLX   and  4>0max. 

For  .5  ^  $  <  .707,  <fy  <  0,   <$,  <  0,   therefore  from  Equation 
^  1  max  is  real  positive  for  all  a2f   In  order  for  &Vt3ia>.  '   '     titive 
real  a2t  is  restricted  by,  noting  that    A <  0S 

The  expressions  (3-32)  and  (3-33)  are  compared  as  follows s 

«        ?     a   *>* 

where  .5  +■   $  ^.  .707 


2  tc    ^  Sa  (  I  -  $*)       "       #  (4  $  2-  '  ) 
Since     |  -  $2  >  o  ,   ^  §*-  I  >  0 

Since     ^-i    <0     ,    (5a-lf^6 

and       a,t,  <  \ 

thus  a2tc  defined  by  Equation  (3-32)  gives  positive  real  ft\nax  and 
&>0max  implying  that  Equation  (3-32)  is  valid  for  the  interval 
.5  <  $<    .707. 
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The  validity  of  Equation  (3-32)  for  the  interval  .707  <   ^  <  1.0 
can  be  proved  likewise. 

For  5  =  0,   .5,  .707  and  1.0,  Equation  (3-32)  is  also  valid.   If 
"5=0  the  l0  curves  for  fourth  order  systems  form  the  region  Rq  for  all 
values  of  a2t  implying  that  a2tc  >  co  which  agrees  with  Equation 

If  £  =  .5   then  Equation  (3-28)  becomes 

to  give   a),  -m  a*  -  Ia*i/3 

Let  ^imax  =  <^0max>   tnen  a2tc  =  • 7^  which  is  equal  to  the  value  of 
a2tc  evaluated  from  Equation  (3-32)  with  ^  =  .5 

If  $  -  .707,  the  solution  of  Equation  (3-28)  obtained  likewise 
gives  a2tc  z    .5  which  is  equal  to  the  value  of  a2tc  evaluated  from  Equa- 
tion (3-32)  with  $  =  .707. 

If  $   =  1.0  then  Equation  (3-32)  is  undefined,  but  lim-(a2tc)  = 
.375  which  agrees  with  the  value  of  a2tc  obtained  by  solving  Equation 
(3-28)  with  $  =  1.0 

Thus  a  general  expression  of  a2tc  as  a  function  of  ^  is  obtained. 
Table  3.1  below  which  well  agrees  with  the  curves  in  Fig.  3-3  to  3-22 
shows  &2tc  and  GJt  for  some  pertinent  values  of  "£  .   The  B«  and  E, 
are  the  coordinates  of  the  point  at  a)nt  =  60tc  on  the  1^*  curve  for 


Table  3.1 

£ 

a2tc 

&)tc 

B0c 

Blc 

.05 

.50250E+02 

.4988lEt01 

.62499Ef03 

.25000Et02 

.10 

.12750Ef02 

.24775E+01 

.39050Ef02 

.62490E+01 

.15 

.58056Et01 

.16347E+01 

.77051E+01 

,27751 

.20 

.33750E+01 

.12097Et01 

.24320E+01 

,15594Et01 

.25 

.22500Ef01 

,95238Et00 

.99187E+00 

.30 

.15389E+01 

.77938E+00 

.47532E+00 

.68930E+00 

.35 

.12704E+01 

.65486E+00 

.25443Er00 

.5Q423E+QQ 

.40 

.10313E+01 

.56090E+00 

.14764Ef00 

.3840lEf00 

.45 

.86728Et00 

.48748E+00 

.91110E-01 

.3Q159E*00 

.50 

. 75000E+00 

.42857Et00 

.59038E-01 

.2427: 

.55 

.66322E+00 

.38032E+00 

.39813E-01 

.19^28E<-00 

.60 

.59722E+00 

.34014Et00 

.27762E-01 

,16642Ef00 

.65 

.54586E*00 

.30621E+00 

.19923E-01 

.141G2EJ-Q0 

.70 

.50510E+00 

.27723E+00 

.14661E-01 

•12107E^00 

.75 

.47222E*00 

.25225E+00 

.11033E-01 

.10516Ef00 

.80 

.44531E+00 

.23053Et00 

.84696E-02 

.92325G-01 

.85 

.4230lE*00 

.21151E+00 

.66209E-02 

.81865E-01 

.90 

.40432E+00 

.19475E+00 

.52617E-02 

. 73264E-01 

.95 

.38850E*00 

.17990E+00 

.42450E-02 

.66135E-01 

1.00 

,37500E<-00 

.16667E+00  ' 

.34722E-02 

.60185E-01 
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The  necessary  condition  on  the  coefficient  a£t  far  existence  of  the 
region  R^  on  the  BQt  vs  Bit  plane  and  the  design  charts. 

As  seen  in  the  derivation  of  a2tc  in  terms  of  £>  s  the  '  ^  curve  of 
a  fourth  order  system  for  a  given  ^  on  the  Bq^  vs  3t^  pi 
points  of  extrema  v/hose  locations  depend  upon  the  value  of  a^t  with  re- 
spect to  a2tc  and  the  encirclements  of  the  curve ,  which  are  cJ- 
lated  to  the  existence  of  the  region  R-  ,  depends  upon  the  location 
the  extreme  points. 

Another  feature  related  to  the  existence  of  the  region  R<»  is  the 
behavior  of  the  curve  at  Large  values  of  &)nf   This  depends  upon  the 
value  of  5  and  will  be  described  with  the  basic  equation  of  the  curve. 
Equation  (3-4)  which  is  again  written  below. 

B,t  =  a.zt&(4 )  t&nt  +  <p?(*>)^i    +  fat  <,, ,  cCht 

v  (3-4) 

For  0<^.5,     <&($)>Qt        <%($)>    0  and  <&(<?)<Q     indicate 
that  the  curve  tends  to  infinity  in  the  third  quadrant.   From  Equation 
(3-29)  it  is  seen  that  there  are  no  60-^  min  and  ^Qmin'   Therefore  if 
a2t  —  a2tc   tnen  the  curve  vents  down  in  the  first  quadrant  and  enters 
into  the  third  quadrant  after  passing  through  the  fourth  quadrant,  thus 
the  entire  area  enclosed  by  the  curve  in  the  first  quadrant  Ls  1 
R<  .   If  a2t  is  increased  beyond  a.2tc      then  the  curve  encircles  in  the 
counter-clockwise  direction,  and  there  exists  the  region  R.  provided  the 
curve  does  not  enter  into  the  second  quadrant.   Referring  to  Fig.  3- 28b 

which  is  the  case  of  a2t  >  a2tc   ^or  ^  **  ^  -  ,5>   -^  a2tc  is  n 
much  larger  than  a£tc  then  the  curve  enters  into  the  third  [u 
through  the  fourth  quadrant  after  forming  a  region  R^  in  1 
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quadrant.   (Note  that  in  this  case  the  curve  forms  a  closed  lo 

first  quadrant  and  the  area  enclosed  by  the  loop  is  non-R^  ,)   If  »,,«.  I? 

greater  than  the  value  of  a2tc  such  that  the  curve  enters  into 

quadrant  by  passing  through  the  origin,  then  there  exists  no  region  R  , 

This  value  of  &2t   denoted  by  a2to  can  be  found  by  setting  B^t  -  B  c 

in  Equation  (3-4)  and  solving  for  &2^   »   thus  one  obtains  a.~       in  tern 

^  as 

where  &).«  is  the  frequency  at  B,t  =  BQ  r  0  when  a«t  ~ 
Thus  for  0  <  £>  ^   .5  the  condition  on  a.2t   for  the  existence  of  the 
region  R^  is  0  <£   ^^ <  a?t  *   ^  a2t  aPProac^es  a2t  t*ie  reg£°B  R^ 
gets  smaller  and  it  is  in  the  region  near  to  the  origin. 

Noting  that  Equation  (3-34)  holds  for  04  5^  1. 0a  apd  by  comparing 
Equations  (3-32)  and  (3-34),  it  is  seen  that 

2-ztc     <  a*to      6-or  0    <  $   C  .00-7 

**tc   >  ^to       ^  W<    $    <    U 

The  relation  between  R*  ,  a?  and  a?  for  the  intervals  of  5  » 
.  5  <  £  4  .707  and  .707  <  ^  4  1.0  can  be  analyzed  in  a  similar  way  as 
for  the  case  of  0  <4  <  .5 

Fig.  3-30a  and  3-30b  are  sketches  of  the  cusrves  for  .  5  <  "S  ^  .707 
and  .  707<£  ^  1.0  respectively  with  various  values  of  a£t  .  The  curves 
are  labeled  by  numbers  from  1  to  5  and  a2t  increases  with  the  Labeling 
number.  The  arrow  on  each  curve  indicates  the  direction  cf  increasing 
60 nf  T^e  dotted  portions  of  the  curves  merely  indicate  the  behavioi  of 

159 


4  * 


Fig.  3-3<»a.  Illustration  of  the  behavior  of  the  Botvs.  Blt  curves 
for  .5  <  ^  =.707 
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Fig.  3-30b. Illustration  of  the  behavior  of  the  Botvs.  B±t  curves 
for  .707^  5=1.0 
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the  curves  after  leaving  the  first  quadrant. 

For   .5<*5^  -707,    <#(€)>  0  ,    <k^  <   0  and   <&££)  <  0 
indicate  the  curves  tend  to  infinity  in  the  second  quadrant..   Fn   Equa- 
tion (3-29),  it  is  seen  that  there  is  no  positive  real  fr--,  r 
^lmin  °ut  there  may  exist  ^Omin  f°r   .  5  ^  ^  £  .707.   For  a2t  - 
which  is  the  case  of  curve  1  or  2  in  Fig.  3-30a  the  curve  bends  dowc 
the  first  quadrant,  enters  into  the  fourth  and  then  intc  tl 
drant,  thus  the  entire  area  enclosed  by  the  curve  in  the  first  quadrant 
is  the  region  R*  .   If  a2t  <  a2t  <  a2to  c^en  as  s^own  by 
Fig.  3-30a,  the  curve  encircles  in  the  counter-clockwise  '  the 

first  quadrant  and  enters  into  the  fourth  and  the  third  the 
second,  thus  a  region  R^  is  formed  below  the  earlier  part  of  the  curve. 
If  a2  >  a2to   (curve  3  or  4)  then  the  curve  does  not  enter  into  the 
fourth  quadrant  and  no  R,  is  formed.   Thus  for  .5  ^  $  4  .707,  the  rela- 
tion between  R,,  ,  a0fc  and  a0i_   is  the  same  as  for  0  <  4  ==  .5. 
5    ^tc      ^to 

For  .  707<$~4  i-Oj  the  curves  are  sketched  in  Fig.  3-30b.   Noting 
that  a2to<  a2tc  ,  it  is  seen  that  if  a2t<  a2to   (curve  *  &n&   2)  the?  the 
entire  area  enclosed  by  the  curve  in  the  first  quadrant  is  the  region 
R^  .   If  a^to^  a2t<a2tc   (curve  3  and  4)  then  the  curve  forms 
loop  in  the  first  or/and  the  fourth  quadrant  encircling  in  the  clock 
direction  to  give  the  region  R,  inside  the  loop  in  the  first  quadrant. 
The  curve  5  does  not  form  the  R^  region,  thus  for  <,7G7<4-  1»0  the  con- 
dition on  a2t  for  the  existence  of  the  region  R^  is   0  <£  a^  <C    ^.  . 

In  particular  for  •£  =  1.0   a2tc  =  *^7^  from  Table  3.1  and  the  region 
R-L  q  is  the  region  for  an  over-damped  system.   Thus  the  necessary  condi- 
tion on  a2t  f°r  a  fourth  order  system  to  have  an  over-damped  r 
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R^g  is  0<  a2t  <  .375  which  agrees  with  the  condition  derived  from 

the  statement  in  References  1  and  2,  that  is, 

The  necessary  condition  of  existence  of  an  over-dairipei  region 
for  the  M  point  is  to  have  n-2  distinct  positive  real  roots  of 
the  polynomial  F"(-S)  =  0  where  F(S)  =  0  is  the  characteristic 
equation  and  the  primes  denote  differentiations  with  respect  to  S. 

This  can  be  shown  as  follows: 

For  a  fourth  order  system,  the  characteristic  equation  is 

and      YM  ~   ''Z^-Zd   -t-  2  &2t   =  0 
The  roots  of  the  second  equation  above  are 


In  order  for  the  roots  to  be  two  distinct  positive  real 

\t  <('*sf(6)   =-375" 
The  analysis  of  the  relation  between  R   and  a2t  made  above  can  be  sum- 
marized as  in  Table  3.2. 

Fig.  3-31  to  3-44  show  the  families    I  $,{2^   on  the   5s  »3,  .4, 
.5,  .6,  .7,  .8  and  1.0  planes  for  a£t  up  to  about  a-2tc    •   (For  ^  3  .3 
and  .4  plane  a2t  is  UP  to  1.0.)  Since  in  a  design  problem^  the  value  of 
^  is  chosen  deliberately  at  some  pertinent  values s  and  since  the  curves 
for  any  intermediate  values  of  a2t  shown  in  the  figures  can  be  constructed 
easily  by  means  of  the  constant  Cdnt   lines  (as  shown  previously  in  Fig. 
3-26)  and  since  a2t  =  a2tc  is  sufficient  for  a  curve  to  form  the  region 
R^  ,  the  figures  may  be  used  as  design  charts  for  fourth  order  systems. 

An  inspection  of  the  curves  reveals  that  the  arguments  developed  for 
a2tc  »  a?r   an(*  t^e  reS*-on  K-<  agree  with  the  curves. 
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Table  3.2 

(a2tc  is  given  in  Table  4.1  or  by  Equation  (3-32) 
a2t0  is  8iven  °y  Equation  (3-34)) 


Interval  of  £ 

a2t 

R^ 

0  <   $  $  .707 
a2tc  <  a2t0 

0  *  a2t  £  a2tc 

Whole  area  in  the  first  quadrant  en- 
closed by  the  \~1    curve. 

a2tc<  a2t<  a2t0 

Area  nearer  to  the  origin  in  the 
first  quadrant  enclosed  by  the  IZ    , 
(The  area  in  the  closed  loop  is 
non-R^  ) 

a2t  ^  a2t0 

No  R 

.707<£  £  1.0 
a2tc>  a2t0 

0^a2t^  a2t0 

Whole  area  in  the  first  quadrant  en- 
closed by  the  ("7  . 

a2t0<^a2t<a2tc 

Area  in  the  first  quadrant  enclosed 
by  the  closed  loop  encircling  in  the 
clockwise  direction. 

a2t>a2tc 

No  R 
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3-4.  Application  of  fourth  order  charts. 

Since  the  curves  constructed  in  the  previous  ■ 
fourth  order  system  and  the  nature  of  the  curves  a? ^ 
versal  to  the  fourth  order  systems,  they  can  he 
and  design  of  fourth  order  systems.   It  Must  be  note^i 
plotted  for  some  pertinent  values  of  the  parage tc 
the  curves  depend  and  there  must  be  employed  soate 
to  apply  to  a  specific  system. 

In  this  section,  the  fourth  order  charts  are  appLie  ' 
cal  examples.   The  procedures  and  methods  of  an.-  . 
system  are  the  same  as  developed  in  References  1  &&£   2. 
Example  3.1 

Given  a  fourth  order  system  with  open  loop  tran»f— 

where  K,  Z  and  P  are  variable  parameters. 
The  problem  is  to  analyze  the  effects  of  varying  K,,   1    -  :  I 
system  stability. 

The  closed  loop  transfer  function  is 

_C/0=  ICS  +  K  Z 

and  transforms  into,  by  using  a  transformation  S  s  (2+*P)s, 


ICg 


c  ,Ay  _ Bid3       (£±E)* 


thus  the  transformed  coefficients  are 
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zprK 


(z+r)3 

^'2? 

(2+Ff- 

K 

2,^  = 


The  transformed  characteristic  equation  has  unity  coeffi 

4      3 
the  terms  s  and  s  and  the  charts  can  be  applied , 

Suppose  K  and  Z  are  fixed.   By  varying  P,   a?  s  a^t  and  a.Jf  arg 

varied  indicating  that  both  the  curves  and  the  M  point  move.   Firsi 

sider  the  change  of  the  curves.   Noting  that  the  shaps  of  the  eui 

!   ' 
(Family   !^        )  depends  only  upon  a2t  ,  the  change  of  the  curve- 
can  be  analyzed  by  investigating  the  relation  between  a«t  and  P. 

Assuming  P  is  any  value  greater  than  zero,  from  Equation  (3=38)  a-t 
is  a  function  of  P.  To  see  the  dependence  of  a£t  on  P<,  &2t  is  di££eren= 
tiated  with  respect  to  P 

implies  that  a«t  is  a  monotonically  decreasing  function  of  P.  a«t  i 
mum  is  obtained  at  P  r  0   thus 

In  view  of  Table  3.1,   a2tc  :  .5  is  between  £  =  .7  aad   >^  .: 
(a2t  =  .5  corresponds  to  ^  =  .707).   Noting  that  &it  ^   a2tc  provides 
no  region  of  R^  for   ^  »  .707,  for  very  small  values  of  Ps  there  is  ao 
region  for  the  M  point  to  guarantee  <    ^    .707. 

As  P  increases  from  zero,  a.^   decreases  from  .5  and   appro.-*.:,  he. , 
when  P  tends  to  infinity.   Therefore  if  ^  ^  .707  is  requite]  ■ 


be  considerably  larger  than  zero.   For  example  if  an  over- -dimp&d   tr&ia 
sient  response  is  required,  then  a2tc  =  .375  indicates  that     F>  1.995. 

The  locus  of  the  M  point  with  respect  to  the  1^  curves  must  be 
analyzed  in  order  to  see  the  system  relative  stability.  The  locus  of 
the  M  points  is  defined  by,  from  Equation  (3-38), 


(2   +  p  >3 

(3-39) 


Bm=  = 


Suppose  P  is  fixed,  then  a2t  is  fixed  to  give  a  fixed  set  ©f  \ -? 
The  locus  of  the  M  point  can  be  analyzed  on  the  BQt  vs  B^t  plana  wlhiere 
the  family   |  te\  a  is  plotted. 

For  the  problem  to  be  specific,  suppose  P  is  fixed  to  a  muaber  very 
close  to  zero.   Then  a2t  is  very  close  to  .5  and  the  locus  of  the  M  p&int 
for  varying  K  and/or  Z  is  analyzed  on  the  chart  for  a~t  £  .5.   For  this 
case,  the  locus  of  the  M  point  is  defined  by 


E 


it— y- 


(3-45) 


which  is  obtained  from  Equation  (3-39)  by  neglecting  P. 

Equation  (3-40)  is  analyzed  as  follows. 

If  K  is  treated  as  a  parameter  then  each  value  of  K  determines  the 
abscissa  of  the  M  point  according  to  the  first  of  Equation  (3=40)  and 
ordinate  of  the  M  point  varies  linearly  with  Z  according  to  the  second 
of  Equation  (3-40),  thus  the  M  point  moves  along  a  vertical  line.  Thus 
for  a  set  of  values  of  K,  there  results  a  family  of  vertical  lines  aloag 
which  the  M  point  moves. 
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If  Z  is  treated  as  a  parameter  then  by  eliminating  K  froa  Eq» ■  .- 
(3-40),  the  locus  of  the  M  point  becomes 

J30t:=.5-2Blt 
which  is  a  family  of  lines  with  slopes  of  ,5Z  passing  through  the 
Thus  for  each  value  of  Z,  the  M  point  moves  along  a  line  defined  by  Equa= 
tion  (3-41)  as  K  varies.   In  Fig.  3-45  which  is  a  duplicate  ©f  the 
IB},<32t     f°r  a2t  ~  '•*     t*ie  ^oc^  °^  tae  M  points  (the  family  of  the  lines 
mentioned  above)  are  shown.   If  K  =  1.0  then  the  M  point  moves  along  the 
vertical  line  for  K  =  1.0  as  Z  varies  and  if  Z  s  .10  then  the  M  point 
moves  along  the  slant  line  for  Z  =  .10  as  K  varies s  thus  for  each  value 
of  K  (or  Z),  there  is  a  limiting  value  of  Z  (or  K)  for  a  specified  syitea 
stability.   For  example  if  K  =  1.0  then  for  the  system  to  have   $  ^  .5 
the  value  of  Z  must  be  Z,  <  Z  <  Z^  where  Z,  is  the  lower  bound  of  z  for 
the  M  point  to  stay  inside  the  region  R  c  and  is  found  at  the  lover  inter 
section  point  of  the  vertical  line  for  K  =  1.0  and  the   1.5  carve  and! 
analogously  for  Z«.   In  like  manner  for  a  fixed  value  of  Zs  say  Z  s   .10),, 
K  must  be  less  than  1.5  in  order  for  the  system  to  have  a  stability 
greater  than   ^=  .3  since  the  intersection  point  of  the  lines  for 
K  s  1.0  and  Z  «  .  10  is  about  on  the  curve.   Also  it  is  seen  free,  Fig. 
3-45  that  if  K  is  greater  than  1.5  then  it  would  be  very  hard  to  make 
the  system  have   ^  3=  .3  and  if  K  >   1.52  then  there  is  no  value  of  Z  to 
make  the  system  have  ^>  ^  .3  because  the  maximum  abscissa  of  the  ■  j, 
curve  is  about  .19  (K  =  (8)  (.19)  =  1.52) 

Thus  for  P  r  0  the  system  stability  for  varying  K  and/or  Z  is  ana- 
lyzed.  For  any  value  of  P  the  analysis  can  proceed  in  the  same  way. 

Now  suppose  the  system  given  by  Equation  (3-36)  is  obtained  from  & 
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third  order  system  by  cascading  a  filter  network,  that  iss  given  &  third 
order  system  open  loop  transfer  function 

Cascade  a  compensator  represented  by 

6rJfc\=       |CciS"!f)  (3=43) 

If  K  =  4  then  the  system  (3-42)  is  at  the  stability  limit  which  cam  be 
readily  seen  from  the  third  order  chart  and  has  the  velocity  constant 
Kv  =  2. 

Consider  the  problem  of  designing  a  lag  compensator    ik   is  repre- 
sented by  Equation  (3-43)  so  that  the  system  stability  is  greater  than  in 
terms  of  ^  ,    £=  .5,  with  minimum  allowable  velocity  constant  Kv  =  ^° 

The  problem  is  to  determine  K,  Z  and  P  in  Equation  (3-36)  so  that 
the  specifications  are  met.   By  following  the  analysis  made  above  and 
from  Fig.  3-45,  K  =  1.0,  Z  =  .1  and  P  ==  0  will  give  the  solution.  The 
value  of  P  can  be  determined  from  the  Kv  restriction  or  in  a  lag  eompem= 
sator  P  =  . 1Z  (which  is  a  reasonable  pole-zero  ratio  as  illustrated  in 
Reference  1)  gives  P  =  .01  and  this  P  provides  Kv  -  5  which  is  acceptable 

in  this  problem,  thus  the  compensated  system  is 

l.o  (S  +.Q 


<j(S)  = 


S  (S^+ZS  rX)(s-+.oi) 


Example  3.2 

Consider  a  unity  feedback  servo  system  with  its  transfer  function 

C  fc-  \  -    |QQ° 

and  it  is  desired  to  design  a  lag-lead  compensator  in  order  for  the  sys- 
tem to  meet  the  following  performance  requirement; 


1)  The  velocity  constant  Ky  is  to  be  maintained  at  100. 

2)  The  settling  time  is  to  be  less  than  .5  sec. 

3)  The  system  is  to  have  a  "$  greater  than  .7. 

By  cascading  a  lag- lead  compensator  to  the  original  system  the  com- 
pensated system  becomes  of  fourth  order  and  the  transfer  function  is 

where  Z^,  P^  and  Z£  ,  Po  are  for  the  lag  and  lead  sections 
respectively  of  the  compensator. 

There  are  five  parameters  K,  Z^,  P, ,  Z£  and  P2  to  be  determined. 
K  can  be  determined  from  the  Kv  specification  by  choosing  a  suitable  at- 
tenuation factor  of  the  compensator  and  P,  the  pole  of  the  lag  section 
can  be  neglected  in  fcitrovic's  method  without  danger,  thus  there  reamin 
three  parameters  Z^,  Z£  and  ?2     to  be  determined.   (This  kind  of  problem 
is  solved  in  Reference  1  and  in  this  example  it  is  to  show  how  the  fourth 
order  charts  can  be  applied.) 

If  the  compensator  is  chosen  such  that  the  zero  to  pole  ratio  of  the 
lag  section  is  to  be  the  reciprocal  of  the  zero  to  pole  ratio  of  the  lead 
section,  then  K  must  be  set  to  1000,  thus  the  transfer  function  which  is 
to  be  considered  in  the  Mitrovic's  method  is 

to  give  the  characteristic  equation 

p(£)  =  £*+0O+B )$3i</°°0+ioPz)$\  1000(2,  +-2a  )£  +  1 000 H,^=  0 
with   the   restriction  on  zeros  and  poles  of   the  compensator  as 

1 '  V  ft  J      x 
thus   there  remain   three  parameters    to  be  determined. 
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There  seems  to  be  two  ways  of  applying  the  fourth  otS&r   charts  to 
solve  the  problems   (1)  by  choosing  a  suitable  P2  »  the  coefficients  t^j, 
a3  and  a^  of  the  characteristic  equation  are  determined  to  give  a  fixed 
set  of  Bq  vs  B^  curves   (for  various  values  of  £  )  and  then  one  consults 
with  the  charts  for  the  families   IfelAt   (Fig.  3-3  to  Pig.  3-22) 3  and 
(2)  without  choosing  P2  first  then  the  coefficients  &2  &n,i   a*  °'£  £fe 
characteristic  equation  are  variable  and  one  consults  with  the  eharts 
for  the  families   I  £,{<32t}    (Fig.  3-31  to  Fig.  3-44).   In  this  eaeample 
the  problem  is  solved  by  the  first  method  and  in  the  next  example  the 
second  method  is  shown. 

In  order  to  apply  the  fourth  order  chart ,  the  system  closed  Icop 
transfer  function  is  transformed  into,  by  using  the  transformation 
S  =  a3s  =  (10  +  P2)s  , 

tfefd  +  ^t^d  +  ao1 


(3-44) 


(10  + S)' 

From  the  design  specification  (3),  the  M  point  must  be  inside  a  R  -? 
region.   If  one  chooses   £  =  .8  (R  g  is  inside  Rj   for  a  fixed  &2t  ^ 
R  3  exists)  then  from  the  design  specification  (2)  G)n   the  frequency  at 
the  M  point  on  the   1.8  curve  must  be  greater  than  10  assuming  the  set- 
tling time  is  approximately  4-/  £  Qjn  ,  and  this  frequency  iss  in  the 
transformed  coordinates  system, 

CiV=aVd3  -  10/00 +  RO 
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p^r  ^+(d3  +  aa.-fc/t+alt^  +  aat 

,.vp,r..       ~       ..      1000+  \0  fi                          1 

I  0O0 

UW£r<_     dgt-  (10+Fk)*                   ^*~ 

lO  +  ft 

|C?00(2,+€.) 

a/5+—       /.  -     ,     r-,     xo. 

d it ions  on  a2t  are 

R-8      aat<.4<*5 
R.*  azt  < .  vov. 

R,.o         ait<.375 
For  the  existence  of  R  g  ,  the  equation  for  a2  in  Equation  (3-44)  and 
the  above  condition  on  a2t  requires  ?2  ^   49. 

If  one  chooses  &2t   =  .3  which  is  a  sufficiently  small  number  so  that 
it  provides  a  good  flexibility  to  locate  the  M  point  to  guarantee  £  £ . 7 
then  P2  =  63.85,  thus  one  obtains  the  following  data 

6^>.I37  or;  a  Cg  WfW 
At  this  point  it  is  checked  if  the  value  of  P2  determined  on  the 
basis  of  the  region  Rg>  is  sufficiently  large  so  that  the  larger  of  the 
newly  introduced  roots  of  the  system  is  much  greater  than  the  real  parts 
of  the  dominant  conjugate  roots. 

p 

Next  one  consults  with  the  family    M^K^t    f°r  a?t  r  '^   wn*-cn 
is  in  Fig.  3-8.   Fig.  3-46  is  a  duplicate  of  Fig.  3-8.   The  curves  for 
"C  z    .7,  .8,  .9  and  1.0  with  frequencies  on  each  curve  are  available. 
It  is  seen  on  the  chart  that  the  line  B^t  =  4  X  10" 2  is  very  close  to  the 
[78  curve  for  .150  <  CO   t  <■    .200  which  is  an  allowable  range  of  C0nt   at 
the  M  point  of  this  problem.   Thus  the  abscissa  of  the  M  point  is  set  to 
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Fig.  3-46  Solution  of  Example  3.2 
with  the  family  ||^,a? 
for  a0,=.3 
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This  a^t  and  P2  substituted  into  Equation  (3-44)  gives 

2,-1-^  =  16 
Let  Z2  =  .1P2  tl»en  Z^  =   6.385  and  Z±   =  9.615  to  give 

From  Equation  (3-44)  a   =  20.9  X  10~4  thus  the  M  point  is  located  at 
(alt  ,  aQt)  =  (4  X  lO-2  ,  20.9  X  10'4)  and  this  point  is  at  C0fl£  =  . 150 
on  the  Hg  curve  which  is  a  satisfactory  M  point  for  the  problem. 

Since  Z~  -    • 1  Po  *-s  chosen,  then  Pi  =  .1  Z,  must  be  chosen  in  view 
of  the  chosen  compensator  pole-zero  ratios,  thus  P^  -  „9615  and  the  com- 
pensated system  transfer  function  is 

1000  (.S+<?.6l5)(S  +  6.3#5-) 


m=- 


S(S+\0)(S+.<Ui&(&  + &Z8S) 


tion  between  a^t  and  (Z^-t-  Z2)  as  seen  from  Equation  (3-44),   (Z-^-t-  Z2) 
associated  with  alt  =  4.25  X  10"2  is 

again  using  the  proportionality  relation  between  a~  and  ZjZ2 

and  so  determined  M  point  denoted  by  M2  in  Fig.  3-46  is  also  satisfactory, 
then  the  compensated  system  transfer  function  is 

&C  ;~  S(S+\0)(S+\.o62XS  +  63.85) 
Thus  if  a  first  try  is  not  satisfactory  then  a  second  try  can  be  made 
without  too  much  calculation  and  the  direction  whether  to  increase  or  de- 
crease in  a^t  in  the  second  try  can  be  determined  by  inspecting  the  re- 
sults of  the  first  try. 
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Example  3.3 

Consider  the  same  problem  of  Example  3.2  with  the  £  specification 
as  t,  z  1.0  with  the  settling  time  specification  omitted.  The  problem 
can  be  solved  on  the  charts  of  the  families    'S.-f^t} 

By  choosing  some  suitable  pole-zero  ratios  of  the  compensator $  Z^ 
and  Z2  can  be  expressed  in  terms  of  P2  which  is  the  only  parameter  to  be 
determined. 

By  choosing  a  compensator  such  that 

p{   =,001  £ 

then  from  Equation  (3-44),   a   ,  a   and  a    which  must  be  determined 
are  expressed  in  terms  of  P2  only,  namely 

3  __  loootioP 


a  =  __i£ 

The  latter  two  of  Equation  (3-45)  define  the  locus  of  the  H   points 
on  the  BQ  vs  B,   plane  if  a  set  of  continuous  values  of  P2  are  supplied, 
and  this  curve  represents  all  locations  of  the  M  points  (among  which  a 
suitable  M  point  must  be  selected)  associated  with  the  compensator  whose 
pole-zero  ratios  are  set  as  above. 

As  P2  varies,  the  location  of  the  M  point  as  well  as  the  \^  curve 
on  the  ^  =  constant  plane  varies.   Since  on  a  ^  =  constant  plane  the 
variation  of  the  I-  curves  can  be  analyzed  by  the  constant  &)nt  lines 
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and  since  the  location  of  the  M  point  moves  along  the  curve  defined  by 

the  latter  two  of  Equation  (3-45),  one  can  find  a  value  of  P2  =  ?2     such 

that  at  P^  ,  the  M  point  is  on  the  1^  curve  for  a2t  determined  by  P2  . 

To  plot  the  locus  of  the  M  points  defined  by  Equation  (3-45),  the 

following  table  is  constructed  by  remembering  that  for  ^  =  1.0  , 

a0  _  =  .25  a~_  =  .375 
2t0      1  2tc 


p2 

70 

80 

90 

100 

110 

120 

130 

140 

a3  =  10+  P2 

80 

90 

100 

110 

120 

130 

140 

150 

a2t 

.265 

.222 

.19 

.165 

.146 

.13 

.117 

.10 

alt 
(10-3) 

15.1 

12.2 

9.9 

8.27 

7.05 

6.03 

5.25 

4.02 

a0t 
(10-5) 

11.9 

9.75 

8.1 

6.85 

5.85 

5.05 

4.4 

3.44 

The  locus  of  the  M  point  is  plotted  in  Fig.  3-47  from  the  table 
above.  Fig.  3-47  is  the  chart  of  the  family    |e,^t\   ^or  a2t  —  -^ 
(It  is  a  duplicate  of  Fig.  3-44).  The  i^g  curves  for  a2t  =  .2,  .18, 
.16,  .  .  .  are  as  shown.   The  frequencies  on  the   1 1.0  curves  for 
a2t  =  *2>  »18  and  .16  are  labeled  whereas  the  frequencies  on  the   Ij.  q 
curves  for  a2t  ^  .14  are  not  shown.   However  the  frequencies  on  any  curve 
can  be  determined  by  using  the  constant  &>nt  lines;  for  example,  if  one 
constructs  a  constant  &>nt  line  for  (O  t   =  .06  by  drawing  a  line  which 
passes  through  the  two  points  each  at  <Dnt  =  .06  on  the   |^  q  curves  for 
a2t  r  •  I6  an(i  a2t  *  *18  tnen  tne  line  crosses  the  I^q  curve  for 
a2t  =  .14  at  the  point  marked  by  X  sign  and  this  point  is  of  6)nt  r  .06 
on  the  IY.0  curve  for  a2t  =  »14« 
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On  the  chart  with  the  locus  of  the  M  points  plotted  as  shown,  the 
relation  between  the  family    I  $,{%&}    an<J  tne  locus  of  the  M  points  is 
analyzed.   On  each  point  on  the  locus  of  the  M  points  there  is  associated 
a  value  of  P2  which  determines  a2t  .   For  example  if  P2  =  80  then  , 
a2t  =  .222.  Noting  that  the   lj  Q   curve  for  a2t  =  .222  will  be  shifted 
to  the  right  of  the  curve  for  a2t  =  .2,  if  P2  =  80  is  chosen  then  the  M 
point  will  be  inside  the  I 1.0  curve  for  a2t  =  .222  implying  that  P£  r  80 
makes  the  system  over-damped.   Consider  the  M  point  for  P2  =  120.   For 
?2   =  120  &2   s  '13  and  the  M  point  will  be  outside  of  the   lj_  q  curve  for 
a2t  =  .13  since  in  Fig.  3-47,  the  'i#o  curve  for  a2t  =  . 14  is  very  close 
to  the  M  point  for  P2  =  120  and  it  can  be  estimated  that  the  li.o  curve 
for  a2t  s  .13  will  be  between  the  \\tQ   curves  for  a2t  =  .12  and  a£t  =  .14, 
thus  P2  r  120  makes  the  system  under-damped  with  a  fairly  high  ^  .  Con- 
sider the  M  point  for  P2  =  100  which  gives  a2t  =  .165.   Let  the  point  Q^ 
and  Q~  be  the  points  at  which  the  locus  of  the  M  points  intersects  with 
the  I  i#Q  curves  for  a2t  =  .16  and  a2t  =  .18  respectively.   a2t  of  the 
\  ■,    q  curves  varies  from  .16  to  .18  between  the  points  Q^  and  Q2.  Noting 
that  the  portion  of  the  locus  of  the  M  points  between  Q^  and  Q2  is  nearly 
a  straight  line,  and  that  the  M  point  for  P2  =  100  is  at  about  1/4  way 
from  Q^  to  Q2,  by  using  a  linear  interpolation,  one  can  predict  that  the 
jl.O  curve  for  a2t  =  .165  will  pass  through  the  point  M  (alt  ,  aQ  )  which 
implies  that  if  P2  =  100  is  chosen  then  the  M  point  will  be  on  the   I^q 
curve,  thus  P2  r  100  gives  the  solution. 

The  frequency  at  the  point  M  (alt  ,  aQ  )  can  be  determined  with  the 

constant  &)nt  line  for  6)nt  z   .09  as  drawn  in  the  figure  which  shows 

that  U)   t  (or  cT.  )  =  .09.   From  the  relation  between  transformed  and 
nt  x     c  ' 


un trans formed  system  frequencies  one  can  obtain 

which  is  the  value  of  the  double  real  roots  of  the  compensated  system, 
and  the  compensated  transfer  function  is 

r  r<w         IOOO(S  +  l.oVS  +  IO') 

which  happens  to  be  a  cancellation  compensation. 

It  must  be  noted  that  one  can  construct  the  L  g  curve  for  a2t  = 
.165  (which  is  for  ?2  -    100)  by  using  the  constant  Q    t   lines  to  notice 
that  the  point  M  (alfc  ,  aQt.)  in  Fig.  3-47  is  slightly  inside  of  the  [^  Q 
curve  for  a^   r  .165.   The  effect  of  increasing  P2  on  the  system  £, 
(against  £>  =  1.0)  can  be  seen  from  the  figure,  that  is,  by  increasing 
p   from  70,  the  system  ^  is  decreased  from  £  >  1  and  it  assumes 
€,  ==  1  at  Pa  :  100.   The  same  analysis  can  be  made  for  ?2  <  70  if  the 
locus  of  the  M  points  is  plotted  on  the  <£  r  1.0  plane  for  that  range  of 
P9. 
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4.    Design. 

In  general,  design  or  compensation  of  a  given  system  is  to  reshape 
the  system  equation  by  using  some  compensators  so  that  the  system  roots 
are  placed  at  suitable  locations  on  the  S-plane  with  due  considerations 
of  the  specified  performance  requirements. 

In  the  usual  design,  the  steady  state  characteristics  are  specified 
by  designating  the  minimum  allowable  steady  state  error  which  can  be 
interpreted  in  terms  of  the  velocity  error  constant  Kv  (or  K  ,  K«  etc. 
depending  upon  the  system  type  number)  and  the  transient  response  char- 
acteristics can  be  interpreted  in  terms  of  the  system  dominant  roots 
locations. 

In  this  section,  assuming  that  the  design  specifications  can  be 
interpreted  in  terms  of  the  steady  state  error  constants  and  the  loca- 
tions of  the  dominant  roots  on  the  S-plane,  an  analytical  as  well  as 
graphical  method  of  designing  cascaded  compensators  is  developed  by  using 
fundamental  properties  of  the  basic  Mitrovic's  equations  and  curves  in 
order  to  minimize  the  amount  of  trial  and  error  which  is  usually  encoun- 
tered in  design  problems  covered  in  the  literature.   Also  it  turns  out 
that  the  labor  can  be  considerably  economized  and  the  difficulty  of 
choosing  some  suitable  values  of  poles  of  the  compensators  (which  is 
essential  in  the  current  Mitrovic's  method)  is  somewhat  removed. 
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4-1.  Design  of  single  section  compensator. 

Consider  an  open  loop  transfer  function 

where  KQ  is  the  minimum  required  forward  gain  and  D(S)  is  a 
polynomial  of  S  with  all  positive  coefficients, 

and  the  design  specifications  require  that  the  system  dominant  roots  be 

located  at 

without  hurting  steady  state  performance  which  is  specified  by  the  velo- 
city constant  Kv  (or  Ka,  K-,  etc.  according  to  the  system  type  number). 

Since  the  gain  is  adjusted  to  Kq  the  Mitrovic's  working  point 
M(a,  ,  aQ)  associated  with  Equation  (4-1)  insures  the  specified  velocity 
constant  Ky  but  in  general  it  does  not  give  the  roots  at  the  specified 
locations.   In  order  for  the  system  to  meet  both  Ky  specification  and 
the  desired  roots  locations  the  point  M(a^  ,  slq)   which  guarantees  the 
specified  velocity  constant  must  be  at  &)n  =  &)ni  on  the   S  ,  curve, 
and  moreover  the  point  at  CO     =  od  ,      on  the   $,  curve  must  be  such  a 
point  that  it  guarantees  the  system  ^  not  less  than  5i  • 

Let  the  Bq  vs  Bi  plane  be  the  one  on  which  the  system  of  interest 
is  analyzed,  and  let  the  Mitrovic's  variables  for  the  original  system  in 
Equation  (4-1)  be  B?  and  BQ  with  the  superscript  "o"  to  denote  the 
original  system.   Then  the  Mitrovic's  curve  for   5=  ^i  °f  tne  original 
system  is 
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by  choosing  the  coordinate  axes  such  that  b£  axis  ===  B,  axis  and  Bq 
axis  ===  Bq  axis,  the  curve  in  Equation  (4-2)  can  be  constructed  on  the 
Bq  vs  B^  plane  and  this  is  denoted  by   I  £  in  Fig.  4-1  with  the  super- 
script "o"  to  indicate  for  the  original  system. 

Let  the  coefficient  of  the  S   term  in  the  polynomial  D(S)  be  d 
assuming  a  type  one  system  in  Equation  (4-1).   Then  the  pair  (d.  ,  K-.) 
determines  the  location  of  the  Mitrovic's  working  point  on  the  BQ  vs  B, 
plane  and  it  is  located  at 

as  shown  in  Fig.  4-1.   The  radius  vector  0M°  has  a  fixed  slope 

dt  ~Kv 

and  this  vector  is  denoted  by  L^ 

The  point  at  63n  =  4)  ■,   on  the   "£   curve  is  located  at  the  point 
denoted  by  Mj  which  in  general  does  not  coincide  with  M°  requiring  a  com- 
pensation.  Thus  the  point  M^  is  the  required  location  of  the  Mitrovic's 
working  point  on  the  \-e,     curve  in  order  to  meet  the  roots  specification 
and  the  coordinates  of  Mj  on  the  Bq  vs  B^  plane  are,  by  letting 
0)n  =  6)nl  in  Equation  (4-2), 

The  radius  vector  OmJ}  has  a  slope  B§  (  ^  ,  0)nl)/  B°  ($i>&n±)° 
The  Yl.  curve,  M°  and  M°|  points  configuration  depends  entirely  upon 
the  original  system  and  the  design  specifications. 

The  ultimate  goal  is  somehow  to  make  the  points  M°  and  M^  coincide 
since  if  M°  =  M°  then  both  K^   and  the  roots  specifications  are  met. 
Since  there  is  no  way  to  move  M°  and  M°j  in  Equation  (4-1),  a  filter 
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Fig,  4-1.  The  locus  of  the  point  Ng  at(^=^na 
on  the  Bo  vs.  B]_  cxirve  for  X=  ¥l 
of  the  compensated  system. 
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network  single  section  compensator  is  cascaded  to  the  original  system  to 
give  a  compensated  system  equation 

6g)»   P(^P)  <4"3> 

where  Z  and  P  are  zero  and  pole  of  the  compensator  and  are 
variable  parameters.  K  is  the  new  gain  which  has  to  be 
adjusted  to  meet  the  steady  state  performance  requirement. 
Let  Ck    be  a  number  such  that  Z  :R4  ,  then  from  the  Kv  requirement 
K  must  be  adjusted  such  that 

Substitution  of  this  K  into  Equation  (4-3)  gives  the  compensated  system 
equation  as 

^;   D(S)(s  +  p)  l  ; 

where  Z  is  expressed  in  terms  of  P  and  o(  . 

The  compensated  system  in  Equation  (4-4)  has  its  own  Mitrovic's 
curves  and  Mitrovic's  working  point.   By  cascading  the  filter  network, 
the  significant  points  and  curves  for  the  original  system  on  the  Bq  vs 
B^  plane  are  changed  into  those  of  the  compensated  system.   Those  changes 
can  be  analyzed  mathematically  as  well  as  graphically. 

Let  Mitrovic's  curves  and  the  points  (corresponding  to  M°  and  Mj  of 
the  original  system)  be  superscripted  by  "c"  for  the  compensated  system. 
The  characteristic  equation  of  the  compensated  system  is 

P(S)(S  +  P)+  k0*(s  +  fM)  =° 

or  SP[S)  +PD(S)+   K**£    +    lr0?=0  (4-5) 
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It  is  seen  from  Equation  (4-5)  that  any  Bq  vs  B^  curve  for  the  com- 
pensated system  is  composed  of  two  component  curves,  one  from  the  poly- 
nomial  SD(S)  of  which  all  of  the  coefficients  are  fixed,  thus  providing 
a  fixed  curve,  and  the  other  from  the  polynomial  PD(S)  which  is  the 
original  system  polynomial,  from  which  the  original  system  curve  is  ob- 
tained, multiplied  by  P.   Let  B^  and  BQ  be  the  Mitrovic's  variables  for 
the  polynomial  SD(S)  and  let  the  coordinate  axes  be  such  that  B^  axis -i 
B^  axis  S  B.  axis  and  Bq  axis  s  Bq  axis  =  B  axis,   then  the  equation 
defining  Bq  vs  B^  curve  for   £  =  5..  of  the  compensated  system  are 

The  interpretation  of  Equation  (4-6)  on  the  Bq  vs  B-^  plane  is  as  fol- 
lows.  First  of  all  the  Bq  vs  B^  curve  for   £  =  ^  of  the  compensated 
system,  namely  the  B0  vs  B?  curve  for  $  =  ^,  on  the  BQ  vs  B,  plane, 
depends  only  upon  P.   Assume  the  Bq  vs  B-I  curve  for  $,  denoted  by   I  £ 
(which  is  obtained  from  the  polynomial  SD(S))  be  constructed  as  shown  in 
Fig.  4-1.   On  the  two  curves    It,       and   !^   ,  the  points  at  any  fre- 
quency CO  =  CO1   can  be  located.   If  the  coordinates  of  the  point  at 

ordinates  of  the  point  at  £0  r  6)^  on  the    $.  curve  componentwise, 

rc 

then  there  locates  the  point  at  0)n  z    CL)^  on  the   '  %     curve.   By  run- 

ning  &)'  from  0  to  60"  ,  the   '  £,     curve  for  the  frequency  range  from  0 

to  CO"  can  be  constructed  which  is  not  shown  in  Fig.  4-1  since  it  is  not 

r-tc 
necessary  to  construct  the   '  £..  curve  at  this  point;  instead  one  focuses 

c  Hc 

attention  to  the  particular  point  M,  which  is  at  a)  z    CO  ^  on  the   <  €, 


curve.  M§  is  a  moving  point  as  P  varies  since  Mj  is  at  6^  =  0Jn|  on  the 

r-7C  pc 

I  £   curve  and  the   »  €        curve  moves  as  P  varies  whereas  the  point  Mj 

rnc 

is  a  fixed  point  at  Od      =  6dn4    on  the   '  <  curve.   Thus  the  point  Mj 
n     ni  j,  « 

is  not  indicated  at  a  particular  location  on  the  BQ  vs  B,  plane  in  Fig. 
4-1  but  the  movement  of  M^  affected  by  P  is  investigated. 

By  fixing  the  frequency  at  0.)     =    GO  *  in  Equation  (4-6)  the  coordin- 
ate of  M^  on  the  Bq  vs  B^  plane  is  determined  by 

B±  =  BjCS,  >  tb&)  +P  B^, ,  ^i) 

which  is  the  equation  defining  the  locus  of  the  M^j  points  on  the  Bq  vs 

B  plane  as  P  varies.  Noting  that  B{  (^-,6J  ,) ,  B'  ("^i»^n£)» 

Bi  (  ^  ,  0)  .)  and  BQ  (  t,    ,  6J  j )  are  fixed  numbers,  Equation  (4-7)  defines 

a  line  on  the  Bq  vs  Bj  plane  such  that  it  passes  through  the  point 

(B[  (  $i>0dni),   BQ*  (  £.,d)  .))   denoted  by  Q  which  is  the  point  at 

^n  =  ^ni  on  the  I*?!  curve  with  a  slope  of  Bo  (  \>  ^ni^l  (  V  ^nl* 
which  is  the  slope  of  the  radius  vector  OM^  .  This  can  be  visualized  by 
eliminating  P  in  Equation  (4-7)  to  obtain 


*-*«•*->  H^ISfa-* 


(?,>*W> 


thus  the  locus  of  M,  points  defined  by  Equation  (4-7)  is  a  line  parallel 
to  OMj  and  passes  through  the  point  Q  (B[  (  -£  ,  a)  ^) ,  BQ  (  ^    ,  0)  ^)) 
on  the  Bq  vs  Bj^  plane.   This  is  denoted  by  L  in  Fig.  4-1.   At  each  point 
on  the  line  L,  there  is  associated  a  value  of  P  according  to  Equation 
(4-7)  and  the  direction  of  increasing  P  is  in  the  direction  of  the  radius 
vector  OMj  .   At  P  -  0,  the  point  at  60 n  r  (jj    ,on  the  compensated  BQ 
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vs  B,  curve  for    =  '.,  (Mj  on  I  l-      )   is  at  Q  which  is 

Equation  (4-7).   As  P  increases  from  zero,  M^  moves  along  the  line  L 

the  direction  determined  by  0M°  .   Thus  the  desired  location  of  the 
d 

rovic's  working  point  (to  insure  the  roots  specification)  which 
inally  at  Md  jumps  to  the  point  Q  upon  cascading  a  single  section  fil 
network  with  the  pole  at  the  origin  and  moves  along  tl         as  the 
pole  of  the  filter  network  increases  such  that  the  distance  from  Q  to  M, 
is  proportional  to  P.   If  a  filter  network  with  fini 
then  it  jumps  from  the  point  M,  to  a  point  on  the  li 

It  is  readily  seen  that  the  extension  of  the  line  L  in  th 
direction  of  the  vector  OM^j  beyond  the  point  Q  .. 
tive  P  which  means  the  pole  of  the  compensate, 
S-plane. 

The  orientation  of  the  line  L  depends  upon  the  locations  of  the 
Mj  and  Q  which  both  may  be  in  any  quadrant  on  the  B^  vs  I 
ing  upon  the  original  system  and  the  design  specifications 

So  far  the  effects  of  cascading  a  single  section  filter  networ 
the  M^  point  are  analyzed.   Remembering  that  in  general  the  Mitrovic's 
working  point  determined  by  the  pair  (a.  ,  a  )  the  coefficients  of  the 
S  and  3   terms  respectively  of  the  characteristic  equation  io  not 
identical  to  the  desired  M  point  located  at  j.     ~  j  n^  on  the 
thus  it  is  necessary  to  adjust  the  variable  parameters  ti   bfing  the  t 
points  together  to  a  suitable  location.   To  do  this  the  behavior  of  Che 
point  determined  by  the  coefficients  a,  and  a^  must  be  investigated. 

The  point  M  is  at  (d,  ,  KQ)  for  the  original  system.   By  cascading 
the  filter,  this  point  is  moved  to  some  other  location  on  the  Bq 
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plane  since  the  coefficients  of  the  characteristic  equation  are  changed. 
Let  the  point  M  be  the  point  for  the  compensated  system  analogous  to 
the  point  M°  of  the  original  system.   Then  from  Equation  (4-5)  the  char- 
acteristic equation  of  the  compensated  system,  the  M°  point  is  defined  by 

MC--=MC('d,P-hfCo^  ,  KoP) 
which  implies  that  the  location  of  the  M°  point  on  the  BQ  vs  B,  plane  de- 
pends upon  both  c<  and  P. 

The  locus  of  the  M  points  on  the  Bq  vs  B^  plane  for  oL    and  P  vari- 
ables is  defined  by 

3±  =d,"P  +  Ke<* 

(4-8) 

By  eliminating  P  in  Equation  (4-8)  and  noting  that  K  =  KLd,  ,  one  obtains 

Bo-KWBi-ICocO  (4-9) 

which  is  an  equation  of  lines  with  the  slopes  fixed  at  K     and  the  B,-axis 
intersections  given  by  K^  which  is  a  linear  function  of  ck   ,  thus  the 
loci  of  the  Mc  joints  are  a  family  of  lines  parallel  to  the  line  lJL  . 
In  Fig.  4-2,  some  members  of  this  family  are  shown  with  the  configuration 
of  the  line  L,  the  vectors  0M°  ,  OM^  and  Q  point  assumed  to  be  the  same 
as  in  Fig.  4-1.  The  line  LJL   (0M°  vector)  itself  is  a  member  of  the 
family  defined  by  Equation  (4-8)  or  Equation  (4-9)  corresponding  to 
ck  z   0.  Any  member  of  the  family  of  the  loci  of  the  Mc  points  can  be  ob- 
tained by  translating  the  line  LA  horizontally  since  the  slope  of  any 
member  is  fixed  at  ICy,  and  for  any  two  members  one  for  d=  q(.  and  the 
other  for  oi  =  o(.2  >  any  two  P°ints  one  on  each  member  read  the  same 
value  of  P  if  the  two  points  are  at  the  same  ordinate,  which  is  readily 
seen  from  the  second  of  Equation  (4-8). 
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Pig.  4-2.  The  family  of  the  loci  of  the  Mc  points 
ar.d  the  locus  of  the  M§  points. 
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Let  the  line  L'  be  a  member  of  the  family.  At  each  point  on  the 
line  L'   there  associates  a  value  of  P  linearly  according  to  Equation 
(4-8).   The  point  for  P  -  0  is  on  the  Bj^  axis.   If  h{     is  divided  into 
two  parts  by  the  B^  axis  above  and  below,  then  the  upper  part  is  for 
P>  0  and  the  lower  part  is  for  P<0  .   Now  let  the  line  L,^  rrscve  horizon- 
tally by  varying  d  .   Then  for   °t  >  0  L^  is  on  the  right  of  L'   and  for 
\   0  it  is  on  the  left  of  L'   ,  thus  the  line  L',  is  well  defined  for 

•     .     -00     -'  ! 

From  the  definition  of  the  point  M  ,  and  the  locus  of  M^  points  Ls 
if  the  Mitrovic's  working  point  M  is  chosen  at  any  point  on  the  I 
then  the  compensated  system  has  roots  at  S  =  -  >i  -<*  m  -I   1  =  '     U)     * 
and  from  the  equations  defining  the  line  L^   ,  it  is  obvious  that  if  the 
Mitrovic's  working  point  M  is  chosen  at  any  point  on  Lj   then  the  com- 
pensated system  satisfies  the  Kv  requirement.  Therefore  in  order  for  the 
compensated  system  to  meet  both  the  Kv  and  roots  specifications ,  the  Mit- 
rovic's working  point  M  must  be  chosen  such  that  the  point  is  on  both 
lines  L  and  L'   and  this  implies  that  the  Mitrovic's  working  point  M 
must  be  at  the  intersection  point  of  the  two  lines.   If  one  is  free  from 
the  parameter  P  appearing  at  each  point  on  both  line  L'   and  L,  and  also 
free  from  any  restriction  on  o(  the  intersection  point  will  be  I 
anywhere  on  the  line  L  if  the  two  lines  are  not  parallel ,  since  a*,  can 

be  varied  to  any  number,  and  each  intersection  point  determines  o<  re- 
sulting infinitely  many  numbers  of  d'S.      However  the  lines  Tj^  and  L  both 
bear  P  at  each  point  on  the  lines  according  to  the  linear  functions  of  E 
defined  in  Equations  (4-7)  and  (4-8).   The  proper  point  must  be  chossn 
among  the  infinitely  many  intersection  points  (as  mentioned  above)  so 
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that  at  the  intersection  point  the  values  of  P  on  Ij.  and  L  must  be  the 
same  values.  This  is  mathematically  speaking  the  two  equations s  Equa- 
tions (4-7)  and  (4-8),  must  be  solved  simultaneously.   From  the  property 
of  the  linearity  of  Equations  (4-7)  and  (4-8),  there  always  exists  one 
and  only  one  solution  point.   Even  in  the  case  that  the  lines  1  "     and  L 
are  parallel  to  each  other,  there  exists  a  solution.  Let  the  solution 
point  be  Ms,  then  the  point  Ms  is  determined  by  solving  Equations  (4-7) 
and  (4-8)  simultaneously,  namely 

G£(SMuV) 


(4-10) 


^saslr[B,,Ctn^)^1?(B'(51)^-^< 


where  Ps  and  o^  are  the  value  of  P  and  o(.  at  the  solution 
point  Ms 

By  substituting  ?s    and  ots  back  into  any  of  Equation  (4-7)  or  (4-8)  the 

location  of  the  point  M  on  the  Bq  vs  B^  plane  is 

Mere     B£  BfCS,  >«V)=B*K  »«Vn)'8?tf,  ■  **»rt 

So  far,  any  restrictions  to  be  imposed  on  P,d  and  the  location  of 
^  in  view  of  practical  situations  have  not  been  considered.,  and  P__  and 
oCs  in  Equation  (4-10)  can  be  any  number  from  negative  infinity  to  posi- 
tive infinity  and  the  location  of  Ms  defined  by  Equation  (4-11)  can  be 
anywhere  on  the  Bq  vs  B^  plane. 
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Now  consider  the  restrictions  on  Pc  and  ds   .   In  view  of  the  fact 
that  the  filter  network  to  be  realized  by  a  passive  RC  network,  the 
first  restrictions  on  ?s    and  o<.s  are 
Ps   >0 
^  >o 
By  considering  a  practical  RC  network  either  too  small  or  too  large  d, 
must  be  avoided.   Thus  one  may  restrict  c*.  such  that 

•I  .<  U.  ■<'  10 
which  is  a  reasonable  restriction. 

In  order  for  the  compensated  system  to  be  a  stable  one,  the  Mitro- 
vic's  working  point  must  be  in  the  first  quadrant,  thus  the  M.  point 
must  be  in  the  first  quadrant.   If  the  solutions  given  by  Equations 
(4-10)  and  (4-11)  are  within  the  restrictions  made  above,  then  the  solu- 
tions may  be  the  proper  ones.   If  the  solutions  violate  any  of  the  re- 
strictions above,  then  this  is  an  indication  that  the  original  system 
cannot  be  compensated  to  meet  the  design  specifications  by  a  practical 
single  section  filter  network.   Thus  one  must  try  to  cascade  another 
section  of  filter  network  or  use  some  other  scheme  for  compensation. 

The  restrictions  on  the  solutions  can  be  interpreted  on  the  Bq  vs 
Bi  plane.   The  restriction  that  requires  the  M   point  to  be  in  the  first 
quadrant  is  fulfilled  if 
Ps  >0 

because  for  that  P  and  o^  ,  any  point  on  Lj  is  in  the  first  quadrant. 
The  restriction  .  1  <^  oC  <  10  indicates  that  the  point  Ms  must  be  on  the 
area  bounded  by  the  Lines  L'   and  l'  n  which  are  the  members  of  the 
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family  of  the  loci  of  the  M°  points  for  <=/,  -  .1  and  ^  =  10  respectively, 
and  those  lines  are  shown  in  Fig.  4-2.  Let  Q*  and  QM  be  the  points  on 
L  at  which  are  crossed  by  l',  and  L,q  q  respectively.  Then  noting  that 
the  M  point  must  be  on  the  line  L,  the  allowed  locations  of  the  Ms 
points  are  between  Q  and  q"  on  L.   The  line  L1,  (or  Liq  q or  both)  may 
cross  L  at  a  point  below  the  B,  axis.   If  this  is  the  case,  then  the 
allowed  locations  of  the  Ms  point  must  be  governed  by  the  restriction 
P  >  0. 

The  existence  of  a  solution  point  Ms  between  Q'  and  Q"  on  the  line 
L  can  be  determined  as  follows.  Let  the  point  Q'  read  P  =  P«  and  p  -  p/ 
on  the  lines  L  and  L',  respectively  and  let  the  point  Q"  read  P  =  Pu  and 
P  =  P*  on  the  lines  L  and  Ljq  q  respectively,  then  there  exists  the  point 
M  between  Q'  and  Q"  if  and  only  if 
?  >  $  &vd-  ft  <  ft' 

or      fj  <p;    and.        pu  >pj 

This  is  an  immediate  consequence  of  the  linear  properties  of  the  equations 
defining  L^   and  L,  thus  one  can  check  whether  there  is  a  point  M5  in 
the  allowable  region  graphically;  however  it  is  much  easier  to  check  it 
by  solving  Equations  (4-7)  and  (4-8)  simultaneously  to  get  the  solutions 
in  Equations  (4-10)  and  (4-11),  since  the  equations  to  be  solved  are 
simple  linear  equations. 

When  the  solutions  turn  out  not  to  be  acceptable,  then  by  choosing 
some  suitable  oC  and  P,  say  cA  =  (7^.  ,  P  z   P.  ,  to  get  an  intermediately 
compensated  system 


Kx  {StSLj) 


^  j       DCsHs+R) 


(4-12) 


where  K^  is  the  new  gain  adjusted  to  give  the  specified  Ky 
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one  can  follow  the  same  procedures  as  developed  above  for  the  system  in 
Equation  (4-12),  that  is  the  system  in  Equation  (4-12)  is  treated  as  the 
system  in  Equation  (4-1)  and  the  polynomial  D(S)(S+Pi)  in  Equation  (4-12) 
plays  the  role  of  tne  polynomial  D(S)  in  Equation  (4-1).   By  cascading 
another  section  of  filter  network  to  the  system  in  Equation(4-12) ,  the 
system  of  interest  is  then 

rfew      k(s+2L)&  +  ^  .,  .,, 

w  }        DCO(s+TjXs+p,>  {        } 

The  system  in  Equation  (4-13)  is  essentially  in  the  same  form  as 

the  system  in  Equation  (4-3)  except  that  K  which  must  be  determined 

2 
enters  into  the  coefficient  of  S   term  (affecting  Bq  vs  B^  curves)  of 

the  characteristic  equation  for  the  system  in  Equation  (4-13).   The  de- 
tails for  this  case  including  how  to  choose  Z±   and  P^  in  the  first  stage 
will  be  discussed  later. 

If  the  solutions  are  acceptable  then  there  still  remains  one  and 
the  final  check,  "Are  the  roots  at  S  s  -  S,cOT,i  ±  /  I  -  -$f  ^n\  dominant?" 

and  this  is  done  by  sketching  the    '  $      curve  for  the  compensated  system 
(final)  to  see  if  the  point  M  is  on  the  region  R^  . 

The  examples  shown  below  illustrate  the  arguments  developed  above 
for  the  single  section  cascaded  compensator. 
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Example  4.1 

Consider  the  transfer  function 


Us)  = 


I  00 


which  has  a  pair  of  complex  conjugate  roots  at  S  =  -2.5     9.2.   It  is 
to  be  compensated,  without  reducing  the  velocity  constant  Kv  ■  20,  so 
that  the  system  has  a  pair  of  complex  conjugate  roots  at 

S=  -  S^  i  ±J7~?*  <4™  =  - 1  o.5  ±  j  1 J 
which  is  equivalent  to   ?,  =  .7  and   GJ   =  15 
D(S)  =  S2+  5S  and  dx  =  5,  Kq  r  100.   The  M°  point  is  at  (5,  100)  and 
the  line  joining  the  origin  and  the  M°  point  has  a  slope  equal  to  K^  =  20 
which  is  the  slope  of  the  family  of  the  lines  L^   (the  loci  of  the  Mc 
point). 

The  compensated  system  equation  is,  noting  that  d^Kv  s  100  s  by  using 
Equation  (4-4) 

^    J~    (Sa+HTSVS  +  p) 
The  Bq  vs  B^  curves   for   the  original   system  are  defined  by 

K  -  < 

to  give  the  location  of  the  M^  point  at 

K  ^K!  (^(^-7,  tir,=  \5)  ,  B0°(  S-.7  ,  0)^,5))  m  f-£  (* »,  Z*$) 
thus  the  M^j  point  is  located  in  the  first  quadrant  and  the  radius  vector 
OMj  gives  the  direction  of  the  line  L  which  is  determined  next.   Note 
that  the  points  M°  and  Mj  are  not  identical  and  a  compensation  is  needed. 
The  characteristic  equation  of  the  compensated  system  is 


in  which  SD(S)  r  S3+  5S2  and  PD(S)  =  P(S2»  5S).   The  Q  point  is  deter- 
mined from  the  Bg  vs  B^  curve  for  the  polynomial  SD(S)  evaluated  at 
£  =  .7  and  a)n  =  15,  that  is 

B,'  =  5-  &  (t )  to-*  +  <%  ( % )  co£ 

to  give  j3('(  <,=.  n ,  U),h  ■:  i  >■■)  -----  J  M 

B0'(*>-7    flJn  -  iz  :■     -3600 

therefore  the  Q  point  is  located  at  (-111,  -3600)  which  is  in  the  third 
quadrant.  The  locus  of  the  Mj  points,  the  line  L  is,  by  using  Equation 
(4-7), 

(4-14) 

which  define  the  line  passing  through  the  point  Q  (-111,  -3600)  with  slope 
of  225/21  =  10.7.   Note  the  slope  10.7  is  the  slope  of  the  radius  vector 

c 
The  locus  of  the  point  M  is,  from  the  characteristic  equation  of 

the  compensated  system 

B,  =  5  P  -f  I  00& 

(4-15) 

I  oo  p 

At  this  point  Equations  (4-14)  and  (^-15)  may  be  solved  simultaneously 
to  get  the  solutions  of  P  and  o(.  without  constructing  any  curves  or 
lines;  however  for  the  illustrative  purpose  some  significant  lines  are 
drawn  in  Fig.  4-3.   Noting  that  the  points  M  ,  Mj  and  Q  are  Deeded  only 
for  determining  the  slopes  and  the  direction  etc.  of  the  lines,  those 
points  are  not  shown  in  Fig.  4-3  but  the  necessary  portions  of  the  lines 
are  drawn. 
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The  line  L  is  constructed  as  follows.   From  Equation  (4-14)  the  B^ 
axis  intersection  is,  by  letting  Bq  =  0  to  obtain  P,  =  16  and  by  substi- 
tuting Pj^  into  the  first  of  Equation  (4-14), 

Bi  intersect  =  225"    at    p=|£ 
thus  the  point  Q-  is  at  (225,  0).   From  Q,  the  line  L  is  directed  into 
the  first  quadrant  with  slope  of  10.7. 

The  locus  of  Mc  point  for  c<  z   2.25  denoted  by  l.       ,.  passes  through 
Q,  with  slope  equal  ti  Ky  =  20  since  from  the  first  of  Equation  (4-15) 
by  letting  BL  =  225  and  P  =  0,  ot  is  found  to  be  2.25.   The  line  Lg 
the  locus  of  Mc  for  cL  =  8  is  drawn  by  noting  in  Equation  (4-15)  that 
the  B,  axis  intersection  of  the  line  for  oC  s  8  is  800  and  L'  ^   and  L' 
are  parallel.   Note  that  on  the  two  lines  L?  95  an<^  ^s  »  ^or  eacn  pair 
of  points  at  the  same  ordinates,  one  on  I*2#25  an<*  tne  otner  on  L«  »  tne 
values  of  P  are  the  same.  At  point  Q,  the  lines  Lg  and  L  intersect  and 
read  P!    =  131.5  on  Lg  and  P,  r  74.5  on  L.   At  the  point  Q^   P'  =  0  and 
Pl   =  16,   thus 

implies  that  a  solution  point  M  exists  between  Q^  and  Q^  on  the  line  L. 
Thus  the  compensation  can  be  accomplished  by  a  single  section  lead  filter 
network  with  2.25  <•  c<  <    8  which  is  an  acceptable  one. 

By  solving  Equations  (4-14)  and  (4-15)  simultaneously  or  by  using 
Equation  (4-10),  one  obtains 

?s  =28.8 
0^  =  3.5 
to  give  the  locations  of  the  M^  and  M°  points  at  Mg  (494, 2880) 
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and  gs=^/^s  =8.2.3 

(C=  \00ds=  350 
then  the  compensated  system  is 


(j(^- 


Finally  the  compensated  system  is  checked  on  the  third  order  chart 
or  by  constructing  Bq  vs  B,  curves  for  the  system  to  see  that  it  is  ade- 
quately compensated.   Note  that  the  characteristic  equation  of  the  com- 
pensated system  factors  as 

to  give  the  roots  at 

S  =  -|0.£±j  11 
S=-I2.g 
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Examp le  4 . 2 

Given  an  open   loop   transfer   function 

& 


6b(s)  = 


S(SZ+  2. OS  t  1-2) 

which  is  analyzed  on  the  third  order  charts  to  be  an  unstable  system  with 
a  pair  of  complex  conjugate  roots  in  the  right  half  plane,  design  a  com- 
pensator to  place  a  pair  of  complex  conjugate  roots  at  ^  ;  .3  and 
C0n   =  .6  with  the  velocity  constant  maintained  at  Ky  :  5. 
For  the  original  system  the  significant  points  are 

=  r£  (.<??,.  5T<?)  (4-16) 

The  radius  vector  OMJj  has  the  slope  =  .59/. 95  =  .63 

By  cascading  a  filter  network  the  compensated  system  is 

to  give  the  characteristic  equation 

from  which  the  BQ  vs  B,  curves  for   ^ =  .3  and  the  Q  point  are  obtained 

as 

B,'=  •  7^  CO^  ■+  ) .  2  ?  atf  -  .f  f  <^  ^ 

(4-17) 

thus  the  point  Q  is  at  (.681,  .09)  and  the  locus  of  m£  point  L  is 
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defined  by,  using  the  coordinates  of  Md  and  Q  in  Equations  (4-16)  and 
(4-17), 

B,-.6ffJ  +  P<S95) 

From  the  characteristic  equation  of  the  compensated  system  the  coeffi- 
cients of  the  S  and  S  terms  are  seen  to  be  1.2P  4  6  X.  and  6p  respec- 
tively, thus  the  locus  of  the  M  points  is 

(4-19) 

By  solving  Equations  (4-18)  and  (4-19)  simultaneously  or  by  using  Equa- 
tions (4-10)  and  (4-11)  with  the  proper  substitutions  of  the  coordinates 

o 
of  the  points  M,  and  Q  from  Equations  (4-16)  and  (4-17)  respectively s  one 

obtains 

F|  = .  O  I  £  G 

ots  =  .112  £s  s  f|  /^  ~.m.$  fC=  &oL  =  .  6  73 

and  the  point  M  at  (.705,  .0996)  which  are  the  acceptable  solutions. 
The  compensated  system  is  then 

G(S)  =  -   -^^-'»g) 

V"  '     S(^2-t-2^^-!.2)(o+.0/^') 

and  can  be  checked  on  the  fourth  order  chart  for  the  dominance  of  the 
roots  at   ^  r  .3  and  Cdn   =  .6 

The  closed  loop  transfer  function  is 
c      .673(* +  •'«?) 

from  which  a2t  r  -%-  s  1'^23-  =  .25.   Noting  that  a2t  is  less  than  &ltc 
for  £  z    .3  (see  Table  3.1  in  Section  3-3)  and  the  point  Mc  is  in  the 
first  quadrant,  the  roots  at   £=  .3  and  0)n  =    .6  are  the  dominant  pair. 
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4-2.  Design  of  double  section  compensators. 

In  a  design  problem,  one  may  find  that  a  given  system  cannot  be 
compensated  to  meet  the  specified  performance  requirements  by  a  single 
section  compensator,  and  this  is  the  case  that  the  set  of  solutions  ob- 
tained by  cascading  a  single  section  filter  network  as  in  Section  4-1 
turns  out  to  be  unacceptable  and  a  multi-section  compensator  is  needed. 

Let  the  set  of  solutions  obtained  by  cascading  a  single  section 
compensator  be 

and  suppose  the  solutions  are  not  acceptable.   This  implies  that  the  orig- 
inal system  must  be  cascaded  with  a  compensator  with  more  than  one  sec 
tion,  thus  the  designer  may  try  to  cascade  a  double  section  compensator 
and  proceed  with  exactly  the  same  approach  as  in  the  case  of  a  single 
section  compensation. 

The  unacceptable  solutions  from  the  single  section  compensation  are 
first  analyzed  before  cascading  a  double  section  compensator.   The  unac- 
ceptable solutions  fit  into  one  of  the  following  cases: 
1)   with    |oU  >] 

la)   {^  >0,     dj  >0 


lc)  ?f<o} 


Id)  ?f     <  Di 


2) 

with 

rf\  <i 

2a) 

pf    >o 

(Xdj.  <  .  i 

2b) 

?f    >o 

^  <0 

2c) 

Pf    <o 

eL  >o 

2d) 

Pf    <o 

°V    <0 

The  cases  b,  c,  d  are  not  acceptable  because  of  the  inreallzability  of 
the  compensator  with  a  passive  network.   The  case  la  is  for  the  noise 
problem  and  the  restriction   o< <  10  is  not  a  theoretical  one.   In  case 
2a  which  is  a  lag  network,  if  oC  is  too  small  then  the  sizes  of  the  com- 
ponents of  the  filter  network  may  become  impractical 4  yet  some  oL    con- 
siderably smaller  than  .1  may  be  acceptable. 

Let    <^x(s)=  ^ be  a  fictitious  compensator 

•J  S  "*"  ]  r 

where  Kc  is  the  attenuation  factor  and  Zf    z   P^  lol,      >  P.  and  oL,      fit 
into  one  of  the  above  cases  which  are  the  solutions  obtained  by  cascading 
a  single  section  compensator  and  G     (s)  is  not  acceptable. 

If  Gc+  (S)  is  cascaded  to  the  original  system  G0(S)  in  Equation 
(4-1),  then  the  transfer  function  G0(S)Gcf (S)   satisfies  the  well-known 
roots  locus  criteria  thus 

at   £— ^cOm^/Fs**^  =ri 

Noting  that  K  >  0,  c/.   can  be  either  positive  or  negative   the 
angle  criterion  in  Equation  (4-20)  for  the  designated  root  r^  is 
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Z6b(s) 


4^)s=f->* 


(4-21) 


where  n  is  an  integer. 
Let  *f    be  the  angle  that  must  be  added  to  the  angle  Z.6-0Cs)        by 
some  compensator  in  order  for  the  compensated  system  to  have  roots  at 
the  designated  locations,  then 


From  Equations  (4-21)  and  (4-22) 

1     "^       S~KPf 


(4-22) 


?■■ 


$>+** 


S+Px 


s-r, 


s+iL 


S+Px 


(4-23) 


s=n 


•i     ^.  <° 


Since  0° <  f  <  180°  indicates  that  Gn(S)  needs  a  lead  angle  in  order  for 
the  compensated  system  to  have  the  roots  at  the  designated  locations „ 
and  similarly  since  -180°  <  <f   <  0°  (or  180°  <  f  <   360°)  indicates  a  lag 
angle,  the  unacceptable  solutions  for  P  and  ck     tell  the  designer  the 
type  of  compensator  needed. 

With  the  aid  of  Fig.  4-4  one  can  obtain  the  following  table  which 
shows  the  proper  type  of  compensator  needed  in  view  of  the  angle  cri- 
terion for  the  possible  combinations  of  the  unacceptable  solutions  r/e 
and  Pf  . 
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i.k^i 


.2.K|<1 


i 

, 

r 

lb.  PjX) 

X 

-?1 

\ 

lc.  P-^0 
di>0 


d<o(1  <  .1 


r 

Id.   P1<cO 

~zl 

*i 

r 

2d. 

\ 

pi<°. 

"Zl 

-pi 

Fig.  4-4.  The  possible  pole  zero  conf igurat j ons 
of  the  unacceptable  single  section 
compensator,  (S+Pi/c*!  )/(s+P1) . 
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Pf   and 

*f 

Zf=Pf/oCf 

Indicated 

compensator 

Pf>0, 

oC  >  10 

zf  >o 

Positive        0°^  <f -^  180° 

lead 

Mfl>i 

Pf  >  0, 

^jL<         0 

Zf<  0 

Positive     180°^  y^360° 

lag 

Pf<0, 

^>  o 

Zf  <  0 

Negative   -180°<f  <  0° 

lag 

Pf<  o, 

*f<  o 

zf^o 

Negative  j      0°  <  <f<  180° 

lead 

Pf  >  0, 

X«£.l 

zf  >o 

Negative  :-180°<  y  <  0° 

lag 

Ki<J 

Pf  >  0, 

*<<<    0 

Zf<0 

Positive  |    180°<  ^<360° 

lag 

Pf<  0, 

<%>    0 

zf<o 

Positive        0°  <  *f<  180° 

lead 

; 

Pf  <  0, 

<^:<    0 

Zf>0 

Negative        0    <  <-? <  180° 

lead 

Fig.  4-4  shows  the  pole  zero  configurations  on  the  S-plane  of  the 
unacceptable  compensator  for  all  possible  combinations  of  the  unaccept- 
able solutions  Pf  and  c^-  obtained  from  the  single  section  compensation. 
Thus  the  results  of  cascading  a  single  section  compensator  to  the  original 
system  gives  the  proper  type  of  compensator  to  be  cascaded.   Upon  realiz- 
ing that  the  system  cannot  be  compensated  by  a  single  section  compensator, 
one  may  try  as  a  next  step  to  cascade  a  double  section  (identical)  com- 
pensator in  which  it  is  not  necessary  to  know  the  proper  type  of  compen- 
sator needed,  since  if  the  system  can  be  compensated  by  such  a  compen- 
sator, then  the  solutions  will  now  be  acceptable  and  the  type  of  compen- 
sator will  be  indicated  by  the  solutions. 

Consider  the  original  system  in  Equation  (4-1)  cannot  be  compensated 
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by  a  single  section  compensator  and  an  identical  double  section  compen- 
sator is  cascaded  to  give  the  compensated  transfer  function 

(*&)  =      ,-  ,  N  , V-  (4-24) 

where  D(S),  P  and  ck    are  defined  as  in  Equation  (4-1)  with 
the  same  design  specifications  as  in  Section  4-1.   From  the  steady  state 
performance  specification  the  gain  K  now  must  be  adjusted  to  give 

or        K~  d\  (Cv  cA  =  fc"o<?( 

The  characteristic  equation  of   the   system  in  Equation   (4-24)   with  K  re- 
placed by  Kr  <  can  be  written  as 

RS )  =  MO(S+  pf  +  IW*( S 4  F>;, f -  o 

or  pco«  rPC )  +aps,l>fc)+pbGH  te&*+  2Mps  +  fop*=  0        (4-25) 

An  inspection  of  Equation  (4-25)  shows  that  the  Bq  vs  B^  curve  for 

t?  r  £   associated  with  the  polynomial  F(S)  consist  of  the  curves  for 

2  2 

the  component  polynomials  S  D(S),  SD(S),  D(S)  and  S  .   The  curves  asso- 
ciated with  the  polynomial  SD(S)  and  D(S)  are  computed  already  in  the 
single  section  compensation  as  the  Bq  vs  B^  and  Bq  vs  B^  curves  respec- 


F(S)  r  S  +  a  S+  aQ  ,   for  which  Mitrovic's  equations  are  h-^   ?  $2^°^    » 

?  2 

Bq  =  COy[   •   If  the  curve  associated  with  the  polynomial  S  D(S)  is  de- 
noted by  the  B"  vs  b"  curve  then  the  Bq  vs  B^  curve  for   $  =  ^i  asso- 
ciated with  the  compensated  system  characteristic  equation,  Equation 
(4-25)  is 
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(4-26) 

which  is  a  curve  that  depends  upon  P  and  o(  .   (Note  that  in  the  single 

section  compensation,  the  compensated  system  curve  depends  only  on  P.) 

The  point  M,  which  is  at  03   r  CO   ■■  on  the  compensated  system  curve  in 

Equation  (4-26)  moves  if  either  P  or  ^  is  varied,  and  the  equation  of 
the  movement  or  the  locus  of  the  Md  points  is  obtained  by  fixing  the 

frequency  at  60  n  =  6)  ,  in  Equation  (4-26),  that  is,  the  locus  of  the 

Mj  point  on  the  B-.  vs  B.  plane  is 

Note  that  in  setting  up  Equation  (4-27)  only  the  first  and  last  terms 
need  new  calculations. 

The  coefficients  a  and  a  of  the  terms  S  and  S   in  the  compen- 
sated system  characteristic  equation  F(S)  in  Equation  (4-25)  are 

2  2 

d^P  +  2K..*  P   (The  d-,P  is  obtained  from  the  expansion  of  the  term 

7  2  c 

P  D(S))  and  KqP  respectively,  and  the  locus  of  the  M  point  for  the 

compensated  system  in  Equation  (4-24)  is  defined  by 

3,^diP2>2Ko^p 

a  (4-28) 

Equations  (4-27)  and  (4-28)  can  be  solved  simultaneously  for  c< 
and  P  which  are  the  solutions  for  the  double  section  compensator.   Since 
Equations  (4-27)  and  (4-28)  are  of  second  order  in  P  and  %    ,  the  simul- 
taneous solution  of  the  two  equations  requires  the  solution  of  a  fourth 
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order  polynomial  which  is  in  general  not  simple  to  solve ;  however  there 
are  some  convenient  methods  to  solve  the  equation  such  as  graphical 
methods,  approximation  technique  etc. 

If  the  solutions  from  double  section  compensation  are  not  accept- 
able, then  a  triple  section  compensation  may  be  tried  with  the  same 
procedure  as  for  the  double  section  compensation;  however  the  results 
require  solution  of  a  sixth  order  polynomial  which  may  be  solved  by  a 
digital  computer. 

The  following  example  illustrates  the  double  section  compensation. 

Example  4.3 

Given  an  open  loop  transfer  function 

Compensate  so  that  roots  are  at 

=  .5  ty  =  15 

without  reducing  the  velocity  constant  1^  ;  28 
The  given  system  is  unstable  with  the  complex  conjugate  roots  in 
the  right  half  of  the  S -plane  which  can  be  seen  from  the  third  order 
chart.   If  a  single  section  compensator  is  cascaded,  then  the  system 
transfer  function  becomes 

rrsx ^ZOo([5fP,M) 

b^;      scs+Q(s-riir)(s+p") 

from  which  the  characteristic  equation  is  written  as 

($*+/6S3-i  I  5  S2  )+  p (s3-H6  S2+  i  XS)  -i-4-2  Od  £  4-  ¥  2  o  p  -  0 

The  polynomial  SD(S)  =  S^  +  16S  -+■  15S  gives  the  Bq  vs  B^  curve  for 
£  =  .5  as 
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from  which  [.-.  ,'  ( S  = •- - 5" ,  CO*,  «I5}  «  -  3  I  5"0 

BD; c  t  = •  5";  ca» =i  &) = -5-0  625 

are  obtiined. 

L 

to  calculate   B^OS—^  Hj='5"")  =  ^^^ 


(4-29) 


(4-30) 


Using  Equation  (4-10)  with  the  values  of  b'  b',  B°  and  B°  in 
Equations  (4-29)  and  (4-30)  and  K^  =  Kq  =  420,  one  obtains  the  set 
of  solutions  for  the  single  section  compensation,  namely 

which  are  not  acceptable,  thus  the  original  system  cannot  be  compensated 
by  a  single  section  compensator. 

The  original  system  is  then  cascaded  with  a  double  section  (identi- 
cal) compensator  as 

The  characteristic  equation  is  then 
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(4-31) 


in  which  S2D(S)  =  S5  +•  16S4  +  15S3 

The  B0  vs  Bi   curve  for  t,  =    .5  associated  with  the  polynomial  S2D(S)  is3 

noting  that  the  coefficient  of  the  S2  term  in  S2D(S)  is  zero, 

f3,"=  —I  4  co^f  +  a)„f 

from  which 

B,;Ys-5\g.V,,-I5)=-3373- 

(4-32) 
0^(5-5,  6a„el5)=  VOttSO 

are  obtained. 

Substituting  the  values  in  Equations  (4-32),  (4-30)  and  (4-29)  into  Equa- 
tion (4-27)  the  locus  of  the  Mj  point  for  the  compensated  system  with  the 
double  section  compensator  is  defined  by 

B0=  702750  +2p(r50l*5)  +  ft***)  +  ¥-20^(^5) 

Using  KQ  =  420  and  d^   r   15  in  Equation  (4-18),  the  locus  of  the  M 
points  is  defined  by 

(4-34) 

By  solving  Equations  (4-33)  and  (4-34)  one  obtains  the  set  of  solutions 

dL^v.o9  2s=?r.3Z 

which  are  the  acceptable  set.   The  compensated  system  transfer  function 
is  then  ^ 


S(S  +  \)(S  +  \5)(S+2l.75)z 
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The  dominance  of  the  roots  at  ^  =  .5  and  Cdn  z   15  can  be  checked 

on  the  third  order  chart  by  factoring  the  characteristic  equation  of 

2 
the  compensated  system  into  two  polynomials  (S  +•  15S1-  225)  and  F,  (S) 

where  Fi(S)  is  a  third  order  polynomial,  or  by  sketching  the   I  .5 

curve  of  the  compensated  system  as  shown  in  Fig.  4-5.   Note  that  the 

roots  at  "$=  .5  and  (.0n   =  15  are  not  the  dominant  pair  but  all  other 

roots  have  ^  greater  than  .5. 
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the  compensated  system. (Example  4.3) 
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4-3.  Successive  application  of  single  section  compensator  design 

procedures. 

When  the  original  system  Gq(S)  requires  a  compensator  with  more 
than  one  section,  the  system  can  be  compensated  by  a  non-identical  multi- 
section compensator,  by  simply  repeating  the  single  section  compensator 
design  technique  (Section  4-1)  with  judicious  choice  of  an  intermediate 
root  location. 

Returning  to  the  equations  giving  the  solutions  of  single  section 
compensation 

(4-10) 

where  the  Mitrovic  s  variables  and  the  constants  are  defined 
as  in  Section  4-1.   Pf  and  o(.£  are  unacceptable  solutions 
associated  with  the  designated  root  r.,  =  -S,o)^\  ^/i-S^^m 

It  is  obvious  in  Equation  (4-10)  that  if  $,  and/or  (0    .    are  changed  to 

J-         nl 

t  i 

some  other  values  £,      ,  and/or  CO    then  Equation  (4-10)  may  give  some 

1  nl 

acceptable  solutions  P,  and  o(   and  this  is  equivalent  to  compensating 
the  original  system  to  locate  the  root  at  r^  =  - %](x^\  +J\  ~S1Z  ^7)1 
which  is  different  from  the  designated  root  location  K  —  - 5  (jl>^ t  -h/h^2  Cd^ \ 
thus  the  system  is  partially  compensated  by  so  determined  single  section 
compensator  of  P.  and  o(,  .   To  this  partially  compensated  system,  the 
procedure  of  designing  a  single  section  compensator  with  the  designated 
root  r«   is  applied  and  the  system  may  be  compensated  to  the  desired 
performance  requirements;  if  not,  the  procedure  is  repeated. 

The  choice  of  the  intermediate  root  r^  or  the  choice  of  C,-^   and 
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u:hi     is  based  on  Equation  (4-10).   Since  the  only  restriction  on  P  is 
P>0,  the  choice  of  '£,  and  (x)n\     must  be  in  such  a  direction  that  the 
numerator  and  denominator  of  first  of  Equation  (4-10)  have  the  same 
signs,  and  yet  that  P  with  the  other  numbers  in  the  second  of  Equation 
(4-10)  must  give  an  d.     such  that  it  is  positive  and  not  too  large  and 
not  too  small.   Thus  one  must  consider  too  many  numbers  simultaneously 
which  is  hard  to  define  by  a  general  rule.   However,  in  view  of  the 
fact  that  for  a  smaller  ^  specification  the  compensation  is  in  general 
easier,  by  choosing  a  <,      smaller  than   2^  an^  analyzing  the  l  <» 
and   I  £>    curves  one  can  choose  a  set  of  suitable  d?l     and  d)./n 

Suppose  <~>j       and  0),',,    are  chosen  and  a  suitable  set  of  P.  and  o( 
is  determined,  then  the  partially  compensated  system  transfer  function 
can  be  written  as 

e>'s)-    PCs^s+fO  <435) 

where  KT0/p(S)  =  &a<:~ )   is  the  original  system  defined  as  in 

Section  4-1, 
which  has  the  velocity  constant  at  the  specified  value  Ky  and  has  a  root 
r'  .   To  compensate  the  system  in  Equation  (4-35)  so  that  it  has  a  root 
fr  =  —  $  CUfli  -r-/i— £**  ^7)i         which  is  the  specified  root,  a  single 
section  compensator  is  cascaded  and  the  transfer  function  becomes 

The  characteristic  equation  of  the  system  in  Equation  (4-36)  is  then 

or        F(^)=5aPC5)+f?SP^^p(5PC>)+fi>P(0>)  f  feM^+Osfi^-^P^    (4-37) 

+  I^RP-0 
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With  the  Mitrovic's  variables  for  the  component  polynomials  in  Equation 
(4-37)  defined  as  in  Sections  4-1  and  4-2,  the  Mitrovic's  curve  for 
5  =  ^  of  the  compensated  system  in  Equation  (4-36)  can  be  written  in 
terms  of  the  component  curves,  namely 

By  fixing  M  z    6)  -i  i*1  Equation  (4-38),  the  locus  of  the  M,  points 
(which  is  at  od     =   CO    ■,    on  the   I  e,   curve  defined  by  Equation  (4-38)) 
on  the  Bq  vs  B-,  plane  is  defined  by 

8,  =  &?+  p,  8/  +  P  (  Bi'  +  F5  B ,° )  +  (C0  * ,  o<  ^  ^ )  60,, , 

where  the  Mitrovic's  variables  on  the  right  sides  of  Equation 

(4-39)  are  the  values  of  corresponding  Mitrovic's  equations 

evaluated  at  <f  =  ■c..      and  u)      z   a)    i  • 
'    '  1        n     nl 

Note  that  in  Equation  (4-39)  the  last  terms  are  from  the  component 
polynomial  K0o/,<a  S  in  Equation  (4-37)  as  in  the  case  of  the  identical 
double  section  compensation  in  Section  4-2,  and  if  U-.    =  dL.      and  P,  i  P 
then  Equation  (4-39)  becomes  identical  to  Equation  (4-27) .   Also  note 
that  in  Equation  (4-39),  the  values  b{,  Bq,  B°,  and  Bq  are  computed  at 
the  first  stag*  to  reveal  P_  and  d^    are  unacceptable,  B,  and  Bq  have  to 
be  calculated. 

By  collect 

tion  (4-37)  and  by  doing  likewise  for  the  S   terms  the  locus  of  the  M 
point  on  the  B.  vs  B1  plane  for  the  compensated  system  in  Equation  (4-36) 
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Is   defined  by 

B.^dipi+KoolOP  +  KofioC 

(4-40) 
%=KoP,P 

Equations  (4-39)  and  (4-40)  are  linear  first  order  equations  which  can 
be  solved  simultaneously  without  difficulty  to  give  a  set  of  solutions 
in  P  and  -\    .   If  the  solutions  are  acceptable  then  the  compensation  is 
accomplished,  and  if  not  then  the  procedures  can  be  repeated.   The  follow- 
ing example  shows  the  non- identical  double  section  compensation  which  is 
seen  in  Section  4-2  with  identical  double  section  compensation. 

Example  4.4 

Consider  the  problem  in  Example  4.3  with  the  same  design  specifica- 
tions.  The  original  transfer  function  is 

a*C  J  ~  SCS+IXS+I5) 
which  cannot  be  compensated  by  a  single  section  compensator.   The  poly- 
nomials to  be  considered  in  the  design  problem  are 

p(S)=  S3-M6S2-*-/5*£ 

^0$)=^  +  l&S*  +  '5S3 
to   give   the   associated  Mitrovic's  curves  as   follows; 

B^  (  $,  uO  «  >X4kte)0h  ■+  I  &  4>3(S  )  ti£  -4  (paCf)o)^ 
bJ  (^  U^)^  i  5  tin-  /£&*)*)»    ~  #3  (*)<<>>? 
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b'i  '     -    ■  i*)a).»,a  +  '6^,c  d-V  -t  d\  (^)a^ 

By  cascading  a  single  section  compensator  to  place  the  root  at  t,  ~    „5 
and   60,-  15  resulted  in  an  unacceptable  set  of  solutions 

p  =  -259.8  U  z   -208 

which  indicates  that  the  system  requires  a  lead  compensation.   At  the 
stage  of  getting  the  unacceptable  solutions  the  following  are  computed; 

BQ'  («;  =.  5,  0)^,-/5)  =  -50te5 

Since  the  system  cannot  be  compensated  to  have  the  root  at   £ =  . 5 
and  Cdn  z   15  by  using  a  single  section  compensator,  an  intermediate  root 
is  chosen  so  that  the  system  is  partially  compensated  by  an  acceptable 
single  section  compensator. 

Let  the  intermediate  root  be  at  ^  =  .3  and  6Jn  z   8  then  for  this 

o   o   '      ' 

root  3^,  Bq,  B^  and  Bq  must  be  calculated.  By  using  the  corresponding 

equations,  the  following  are  calculated: 

%l(%<  ■■',  ton*  ^  \=?u?.£> 
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By  substituting  those  values  into  Equation  (4-10)  one  obtains 

which  are  acceptable  solutions,  thus  the  partially  compensated  system 

has  the  root  at   £  r  .3  and   &)n  =  8. 

To  set  the  root  at   £r  .5  and   £0n  =  15  which  is  required  fron 
the  design  specification,  another  single  section  compensator  is  cascaded 
with  the  partially  compensated  system  to  give  the  system  equation 

r  (  <\  =     Z497^(S+3.73)(S+P/o<) 

^  -'}        S  <TS+Q(  5  +  15)  (  S-tZ2./6)(b+p') 

The  new  values  to  be  computed  are 

B^-tT,  0^=150  =  -3375- 

O     O     •     i     ,,         it 

By  substituting  B^,  Bq,  B^,  Bq,  B-^  and  Bq  for   5  =  .5  and  co      =  15  and 

the  values  of  P^  and  oC   into  Equation  (4-39),  one  obtains  the  locus  of 

c 
the  M^  points 

E^  =  -33  75  +  fc2./6;  (■  3/^)  +  p  (-3i50+2Z/6(ZV-o))  +  2^97(JS)d. 

or  &l=-73/'73+2/68  P+37#55<A 

(4-41) 
-  ,      -^4  3  O  /  7  -^  '  •  39 p  +  561225  JL 

The  locus  of  the  Mc  points  of  the  doubly  compensated  system  is,  by 
using  Equation  (4-40), 

v^--2?2°i  p+i3oe* 

(4-42) 
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Equations  (A-41)  and  (4-42)  are  solved  to  give  the  solutions 

*=S.iq  z  =  7-  %& 

which  is  an  acceptable  set. 

The  compensated  system  is 

and  that  the  original  system  is  compensated  to  the  design  specifications 
is  confirmed  by  the  compensated  system  curve   1^5  and  the  Mitrovic's 
working  point  as  shown  in  Fig.  4-6. 
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Fig.   4-6     The    1T5  aiid     L8  curves  for  the 

compensated  system. (Example  4.4) 
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5.   Conclusion. 

1)  Frequency  response. 

It  is  seen  that  the  frequency  response  of  both  closed  and  open  loop 
systems  can  be  obtained  on  the  Mitrovic's  Bq  vs  B^  plane  for  any  order 
system  and  particularly  for  a  system  with  no  zeros,  the  method  is  simple. 

By  using  the  concept  of  the  constant  bandwidth  curves  on  the  Bq  vs 
Bi  plane  a  system  can  be  designed  so  that  it  meets  the  bandwidth  specifi- 
cation closely  when  the  bandwidth  specification  is  a  critical  one.   For 
the  third  order  system,  a  universal  bandwidth  chart  is  constructed  which 
can  be  used  in  both  analysis  and  design  of  a  third  order  system. 

2)  The  fourth  order  charts. 

Universal  fourth  order  charts  on  the  Bq  vs  B^  plane  are  constructed 
which  can  be  applied  in  both  analysis  and  design  of  any  fourth  order 
system.   With  the  preconstructed  fourth  order  charts,  the  labor  of  con- 
structing the  necessary  curves  which  is  essential  in  the  Mitrovic's 
method  is  removed  or  minimized. 

3)  Design. 

An  analytic  method  has  been  developed  designing  the  required  cas- 
caded compensator  by  using  the  basic  properties  of  the  Mitrovic's  equa- 
tions and  the  Bq  vs  Bn  curves.   The  method  gives  accurate  solutions  of 
the  required  compensator  in  order  for  a  system  to  have  its 
roots  at  arbitrarily  designated  locations  on  the  S-plane,  with  arbi- 
trarily designated  gain. 

If  the  given  system  and  the  design  specifications  are  such  that 
the  system  can  be  compensated  by  a  single  section  compensator,  the  prob- 
lem is  to  solve  two  linear  simultaneous  equations  which  provide  solutions 
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quickly.   For  identical  multi- section  compensation  the  problem  is  to 
solve  a  high  order  polynomial  in  which  the  order  of  the  polynomial  de- 
pends upon  the  sections  needed;  however,  the  problem  can  be  guided  into 
the  ones  of  solving  linear  (first  order)  equations  by  applying  single 
section  compensator  design  technique  successively. 
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APPENDIX  I 
FUNCTIONS  Jk(  <S  )    (GENERAL  FORMULA  AND  TABLE) 

0k<$)   =  -<2^k.1(^  )+  0k.2(^» 

0o<^>  -"  ° 

(For  the  proof  of  the  formula  see  Cll  and  [2]) 


l>\ 

l                            1 
0             0.1         0.2         0.3  !      0.4   j      0.5 

0.6 

0.7  j      1.0 

A 

-1          -1          -1         1-1          -1 

-1 

-1 

j  — . 

-1         ]  ~l 

i 

02 

•               ;              i    ' 

0             0.2      j    0.4         0.6     ;    0.8 

1.0 

1.2     I    1.4         2.0 

\  <Pz 

1.0         0.96       0.84       0.64       0.36       0 

-0.44  |  -0.96     -3.0 

i  $4  ; 

0           -0.392  -0.736|-0. 984  -1.088! -1.000 

-0.672   -0.056   +4.000 

\05  , 

-1.000  -0.882  -0.546  -0.050  +0.510 

+1.000 

+1.246 

+1.038!  -5.000 

1  <Pe   \ 

0 

+0.568; +0.954 

+1.014  +0.680 

0 

-0.824 

-1.398;  +6.000 
+  0.918'  -7.000 

k 

+1.000 

+0.768*0.164 

-0.559   -1.054 

-1.000 

-0.258 

* 

0           -0.722-1.020 

-0.679   +0.164+1.000 

+1.133 

+0.112J  +8.000 

:    ^ 

-1.000  -0.624  40.244 

+0.966  +0.923i    0 

-1.102 

-LOTS!  -9.000 

!    ^10 

0            +0.847  +0.922   +0.099   -0.902   -1.000 

+0.189 

+  1.393  +10.000 
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